Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 1,069,320 A→G S463G (AGT→GGT)  purH → bifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/inosine monophosphate cyclohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077931,069,3200AG99.6% 3787.8 / ‑4.0 1390S463G (AGT→GGT) purHbifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/inosine monophosphate cyclohydrolase
Reads supporting (aligned to +/- strand):  ref base A (6/0);  new base G (1366/18);  total (1372/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 1399 total aligned reads displayed.

CAATGATCATGTAGCGTTAGTATCTGATGGATTTTTCCCTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGAGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACATTTTAAACACTAATTTTAAA  >  NC_007793/1069223‑1069419
                                                                                                 |                                                                                                   
cAATGATCATGTAGCGTTAGTATCTGATGGATTTTTCCCTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGt                                                                                                   >  1:4261966/1‑100 (MQ=255)
cAATGATCATGTAGCGTTAGTATCTGATGGATTTTTCCCTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGt                                                                                                   >  1:88399/1‑100 (MQ=255)
cAATGATCATGTAGCGTTAGTATCTGATGGATTTTTCCCTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGt                                                                                                   >  1:2023862/1‑100 (MQ=255)
cAATGATCATGTAGCGTTAGTATCTGATGGATTTTTCCCTATGGGAGATACAGTTGAACTTGCAGCACAACACGGTATAAAGGCAATTATCCAACCGGGt                                                                                                   >  1:1083980/1‑100 (MQ=255)
   tGATCATGTAGCGTTAGTATCTGATGGATTTTTCCCTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGggtggt                                                                                                >  1:1623901/1‑100 (MQ=255)
    gATCATGTAGCGTTAGTATCTGATGGATTTTTCCCTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGtt                                                                                               >  1:2812237/1‑100 (MQ=255)
    gATCATGTAGCGTTAGTATCTGATGGATTTTTCCCTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGtt                                                                                               >  1:3850721/1‑100 (MQ=255)
       cATGTAGCGTTAGTATCTGATGGATTTTTCCCTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGa                                                                                            >  1:4856526/1‑100 (MQ=255)
        aTGTAGCGTTAGTATCTGATGGATTTTTCCCTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCATTTATCCAACCGGGTGGTTCGAt                                                                                           >  1:5455274/1‑100 (MQ=255)
         tGTAGCGTTAGTATCTGATGGATTTTTCCCTATGGGGGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGAtt                                                                                          >  1:1585309/1‑100 (MQ=255)
         agTAGCGTTAGTATCTGATGGATTTTTCCCTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGAtt                                                                                          >  1:3974882/2‑100 (MQ=255)
          gTAGCGTTAGTATCTGATGGATTTTTCCCTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTa                                                                                         >  1:634360/1‑100 (MQ=255)
          gTAGCGTTAGTATCTGATGGATTTTTCCCTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTa                                                                                         >  1:4978315/1‑100 (MQ=255)
          gTAGCGTTAGTATCTGATGGATTTTTCCCTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTa                                                                                         >  1:1000497/1‑100 (MQ=255)
                 aaGTATCTGATGGATTTTTCCCTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCa                                                                                  >  1:4434725/2‑100 (MQ=255)
                     aTCTGATGGATTTTTCCCTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGAt                                                                              >  1:4051791/1‑100 (MQ=255)
                       cTGATGGATTTTTCCCTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTc                                                                            >  1:1827796/1‑100 (MQ=255)
                       cTGATGGATTTTTCCCTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTc                                                                            >  1:3590410/1‑100 (MQ=255)
                       cTGATGGATTTTTCCCTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTc                                                                            >  1:2158993/1‑100 (MQ=255)
                         gATGGATTTTTCCCTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCaa                                                                          >  1:1733768/1‑100 (MQ=255)
                         gATGGATTTTTCCCTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCaa                                                                          >  1:5516612/1‑100 (MQ=255)
                         gATGGATTTTTCCCTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCaa                                                                          >  1:3636209/1‑100 (MQ=255)
                         gATGGATTTTTCCCTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCaa                                                                          >  1:2105617/1‑100 (MQ=255)
                           tGGATTTTTCCCTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAAtt                                                                        >  1:5569718/1‑100 (MQ=255)
                                  ttCCCTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATgg                                                                 >  1:682875/1‑100 (MQ=255)
                                  ttCCCTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGAGTGGTTCGATTAAAGATCAAGATTCAATTGATATgg                                                                 >  1:5360995/1‑100 (MQ=255)
                                   tCCCTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGc                                                                >  1:945309/1‑100 (MQ=255)
                                   gcccTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGc                                                                >  1:2571887/2‑100 (MQ=255)
                                    cccTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACTGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGcc                                                               <  1:4711125/100‑2 (MQ=255)
                                      cTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCtaa                                                             >  1:2689582/1‑100 (MQ=255)
                                        aTGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCtaata                                                           >  1:1396945/1‑100 (MQ=255)
                                        aTGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCtaata                                                           >  1:3292662/1‑100 (MQ=255)
                                                          ctggcagcaCAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATgg                                         >  1:2747065/4‑100 (MQ=255)
                                                           ttGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATGTGGCTAATAAACATGGTATTGCAATGGt                                        >  1:2330656/1‑100 (MQ=255)
                                                           ttGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGt                                        >  1:5028557/1‑100 (MQ=255)
                                                            tGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTGAAGATCAAGGTTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTa                                       >  1:1444800/1‑100 (MQ=255)
                                                            tGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCATTGGTa                                       >  1:1991405/1‑100 (MQ=255)
                                                            tGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTa                                       >  1:4316247/1‑100 (MQ=255)
                                                            tGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTa                                       >  1:370651/1‑100 (MQ=255)
                                                            tGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTa                                       >  1:2628519/1‑100 (MQ=255)
                                                            tGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTa                                       >  1:3247059/1‑100 (MQ=255)
                                                            tGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTa                                       >  1:3901017/1‑100 (MQ=255)
                                                            tGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTa                                       >  1:752918/1‑100 (MQ=255)
                                                            tGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTa                                       >  1:4941481/1‑100 (MQ=255)
                                                            tGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTa                                       >  1:4360952/1‑100 (MQ=255)
                                                            tGCAGCACAACATGGTATAAAGGCAATTATCCAACCGGGCGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTa                                       >  1:5322564/1‑100 (MQ=255)
                                                            cgcagcaCAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTa                                       >  1:5113508/2‑100 (MQ=255)
                                                             gcagcaCAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAg                                      >  1:5153878/1‑100 (MQ=255)
                                                             gcagcaCAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAg                                      >  1:4659019/1‑100 (MQ=255)
                                                              cagcaCAACATGGTATAAAGGCAATTATCCAACCGGGGGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGt                                     >  1:2897975/1‑100 (MQ=255)
                                                               agcaCAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTc                                    >  1:1170571/1‑100 (MQ=255)
                                                               agcaCAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTc                                    >  1:3759555/1‑100 (MQ=255)
                                                               agcaCAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTc                                    >  1:2368649/1‑100 (MQ=255)
                                                                gcaCAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTca                                   >  1:3384353/1‑100 (MQ=255)
                                                                gcaCAACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTca                                   >  1:121190/1‑100 (MQ=255)
                                                                acacaACATGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTca                                   >  1:310419/2‑100 (MQ=255)
                                                                 aacaacaTGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTcac                                  >  1:2975739/2‑100 (MQ=255)
                                                                  acaacaTGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTcaca                                 >  1:1871283/1‑100 (MQ=255)
                                                                  acaacaTGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTcaca                                 >  1:3479837/1‑100 (MQ=255)
                                                                   gaacaTGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAgatcgga                                >  1:5744483/2‑94 (MQ=255)
                                                                   caacaTGGTATAAAGGCAATTATCCAACCGGGGGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAgatcgga                                >  1:5407998/1‑94 (MQ=255)
                                                                      tctGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCa                             >  1:3183454/3‑100 (MQ=255)
                                                                      agtGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCa                             >  1:2854058/3‑100 (MQ=255)
                                                                       gtGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCAc                            >  1:3534156/2‑100 (MQ=255)
                                                                       aTGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCAc                            >  1:2519956/1‑100 (MQ=255)
                                                                        tGGTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACt                           >  1:2433145/1‑100 (MQ=255)
                                                                        aggTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACt                           >  1:1255887/2‑100 (MQ=255)
                                                                         ggTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTc                          >  1:4482122/1‑100 (MQ=255)
                                                                         ggTATAAAGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTc                          >  1:4904274/1‑100 (MQ=255)
                                                                             ctaaGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACa                      >  1:2215179/3‑100 (MQ=255)
                                                                             aaaaGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACa                      >  1:3800125/2‑100 (MQ=255)
                                                                              aaaGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACAt                     >  1:1120381/1‑100 (MQ=255)
                                                                              aaaGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACAt                     >  1:3427655/1‑100 (MQ=255)
                                                                               taGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACAtt                    >  1:1926090/2‑100 (MQ=255)
                                                                               aaGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACAtt                    >  1:1298291/1‑100 (MQ=255)
                                                                               aaGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACAtt                    >  1:3342011/1‑100 (MQ=255)
                                                                               aaGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACAtt                    >  1:1514611/1‑100 (MQ=255)
                                                                                tggCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACAttt                   >  1:5257557/2‑100 (MQ=255)
                                                                                aGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACAttt                   >  1:4153648/1‑100 (MQ=255)
                                                                                aGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACAttt                   >  1:1695300/1‑100 (MQ=255)
                                                                                aGGCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACAttt                   >  1:3006682/1‑100 (MQ=255)
                                                                                 ggCAATTATCCAACCGGGTGGTTCGATTAAAGATTAATCTTCCCTTGCTCTGGCTAATAAACCTGGTATTGCACTGTTAGTCACAGGCACTCGCCCtttt                  >  1:2461131/1‑100 (MQ=255)
                                                                                 agCAATTATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACAtttt                  >  1:2070162/2‑100 (MQ=255)
                                                                                     catATCCAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACATTTTAAac              >  1:3935996/3‑100 (MQ=255)
                                                                                         cccAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACATTTTAAACACTa          >  1:815781/2‑100 (MQ=255)
                                                                                          ccAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACATTTTAAACACTaa         >  1:4614901/1‑100 (MQ=255)
                                                                                          ccAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACATTTTAAACACTaa         >  1:5621711/1‑100 (MQ=255)
                                                                                           gtaCCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACATTTTAAACACTAAt        >  1:4522616/3‑100 (MQ=255)
                                                                                           ctaCCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACATTTTAAACACTAAt        >  1:2287702/3‑100 (MQ=255)
                                                                                           ctaCCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACATTTTAAACACTAAt        >  1:1784040/3‑100 (MQ=255)
                                                                                           cAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACATTTTAAACACTAAt        >  1:5073006/1‑100 (MQ=255)
                                                                                           cAACCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACATTTTAAACACTAAt        >  1:569978/1‑100 (MQ=255)
                                                                                            taCCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACATTTTAAACACTAAtt       >  1:5785161/2‑100 (MQ=255)
                                                                                            gaCCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACATTTTAAACACTAAtt       >  1:1038688/2‑100 (MQ=255)
                                                                                            gaCCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACATTTTAAACACTAAtt       >  1:4783824/2‑100 (MQ=255)
                                                                                            ccccGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACATTTTAAACACTAAtt       >  1:1344037/3‑100 (MQ=255)
                                                                                            caCCGGGTGGTTCGATTAAAGATCAAGATCCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACATTTTAAACACTAAtt       >  1:3046200/2‑100 (MQ=255)
                                                                                             aCCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACATTTTAAACACTAAttt      >  1:5201434/1‑100 (MQ=255)
                                                                                             aCCGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACATTTTAAACACTAAttt      >  1:4094688/1‑100 (MQ=255)
                                                                                              ccGGGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACATTTTAAACACTAAtttt     >  1:4201981/1‑100 (MQ=255)
                                                                                                 |                                                                                                   
CAATGATCATGTAGCGTTAGTATCTGATGGATTTTTCCCTATGGGAGATACAGTTGAACTTGCAGCACAACATGGTATAAAGGCAATTATCCAACCGAGTGGTTCGATTAAAGATCAAGATTCAATTGATATGGCTAATAAACATGGTATTGCAATGGTAGTCACAGGCACTCGACATTTTAAACACTAATTTTAAA  >  NC_007793/1069223‑1069419

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑