Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,460,477 A→C 7.0% intergenic (‑113/+45) nrdG ← / ← nrdD anaerobic ribonucleoside‑triphosphate reductase activating protein/anaerobic ribonucleoside‑triphosphate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,460,4770AC 7.0% 584.2 / 27.9 183intergenic (‑113/+45)nrdG/nrdDanaerobic ribonucleoside‑triphosphate reductase activating protein/anaerobic ribonucleoside‑triphosphate reductase
Reads supporting (aligned to +/- strand):  ref base A (92/78);  new base C (4/9);  total (96/87)
Fisher's exact test for biased strand distribution p-value = 1.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

Only 100 of 186 total aligned reads displayed.

GACGATGTCGACAGGATAGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAATATCCGCACAC  >  NC_000913/4460329‑4460626
                                                                                                                                                    |                                                                                                                                                     
  cGATGTCGACAGGATAGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGc                                                                                                                                                    <  1:1396730/150‑1 (MQ=255)
   gATGTCGACAGGATAGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCg                                                                                                                                                   <  2:515882/150‑1 (MQ=255)
   gATGTCGACAGGATAGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCg                                                                                                                                                   >  1:146308/1‑150 (MQ=255)
       tCGACAGGATAGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTa                                                                                                                                               >  2:664127/1‑150 (MQ=255)
        cGACAGGATAGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTAc                                                                                                                                              >  1:512203/1‑150 (MQ=255)
          aCAGGATAGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTAcgc                                                                                                                                            <  1:191849/150‑1 (MQ=255)
            aGGATAGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTAcgccg                                                                                                                                          <  1:857568/150‑1 (MQ=255)
                 aGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATGCGGCGTAAACGCCTTATcc                                                                                                                                     <  2:101484/150‑1 (MQ=38)
                  gTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATGCGGCGTAAAcgccttatc                                                                                                                                      >  2:334970/1‑143 (MQ=38)
                  gTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATGCGGCGTAAAcgccttatc                                                                                                                                      <  1:334970/148‑6 (MQ=38)
                   tACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                                                                                   <  1:790356/150‑1 (MQ=38)
                   tACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCgg                                                                                                                                   <  1:604638/150‑1 (MQ=255)
                    aCTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc                                                                                                                                  <  1:258121/150‑1 (MQ=255)
                     cTGATGATAATTCATTATTTCGCCTGGCTTATCAGTGCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGCCCTGGAGCAATTGCCAAATGTAAGCCGGATAAGGCGTTTACGCCGCATCCGGCa                                                                                                                                 >  2:781609/1‑150 (MQ=255)
                     cTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCa                                                                                                                                 >  2:429784/1‑150 (MQ=255)
                        atgatAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAttt                                                                                                                              <  2:945748/150‑1 (MQ=255)
                          gatAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGc                                                                                                                            <  1:1238021/150‑1 (MQ=255)
                           atAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCt                                                                                                                           <  1:722869/150‑1 (MQ=255)
                           atAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCt                                                                                                                           <  1:558649/150‑1 (MQ=255)
                           atAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCt                                                                                                                           >  2:1211338/1‑150 (MQ=255)
                           atAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCt                                                                                                                           <  1:289813/150‑1 (MQ=255)
                              aTTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAg                                                                                                                        <  2:1159235/150‑1 (MQ=255)
                              aTTCATTATTTCGCCTGGCGTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTCGCTTAg                                                                                                                        >  2:402757/1‑150 (MQ=255)
                               ttCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGc                                                                                                                       >  1:876181/1‑150 (MQ=255)
                                 cATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTc                                                                                                                     >  1:1265966/1‑150 (MQ=255)
                                  attattTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCt                                                                                                                           >  2:403423/1‑143 (MQ=255)
                                    tattTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCaga                                                                                                                  <  2:1283273/150‑1 (MQ=255)
                                          gCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAg                                                                                                            <  2:174546/150‑1 (MQ=255)
                                           ccTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGa                                                                                                           >  2:1268715/1‑150 (MQ=255)
                                              ggCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGc                                                                                                                            <  1:1170828/130‑1 (MQ=255)
                                              ggCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGc                                                                                                                            >  2:1170828/1‑130 (MQ=255)
                                                cTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAAc                                                                                                      <  2:222626/150‑1 (MQ=255)
                                                 ttATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAAcc                                                                                                     <  2:263114/150‑1 (MQ=255)
                                                 ttATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAAcc                                                                                                     <  2:114147/150‑1 (MQ=255)
                                                   aTCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTa                                                                                                   >  1:675442/1‑150 (MQ=255)
                                                          cgcgTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGccc                                                                                            >  1:216143/1‑150 (MQ=255)
                                                             gtgtAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCAtt                                                                                         >  2:563377/1‑150 (MQ=255)
                                                                  ggCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCa                                                                                    <  2:594159/150‑1 (MQ=255)
                                                                   gCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCaa                                                                                   >  1:516976/1‑150 (MQ=255)
                                                                     tGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAAt                                                                                 >  1:740254/1‑150 (MQ=255)
                                                                      gCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATg                                                                                <  1:374214/150‑1 (MQ=255)
                                                                      gCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATg                                                                                <  1:1085638/150‑1 (MQ=255)
                                                                        aTCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGtt                                                                              <  1:314323/150‑1 (MQ=255)
                                                                         tCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGttt                                                                             <  2:369267/150‑1 (MQ=255)
                                                                          ccGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTa                                                                            >  1:46232/1‑150 (MQ=255)
                                                                          ccGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCAAAATGTTTa                                                                            >  2:684422/1‑150 (MQ=255)
                                                                            ggTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAAc                                                                          >  2:750826/1‑150 (MQ=255)
                                                                            ggTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGAAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAAc                                                                          >  1:414184/1‑150 (MQ=255)
                                                                              tAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAAcgc                                                                        >  1:716835/1‑150 (MQ=255)
                                                                               aaTCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAAcgcg                                                                       >  1:610767/1‑150 (MQ=255)
                                                                                aTCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACgcgg                                                                      >  2:831665/1‑150 (MQ=255)
                                                                                 tCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCTGATAAGGCGTTTACGGCGCATCAGGCATTTTCTTAGCTCAGAGTGAAGATTAATCTATCTGGCCATTCCCTAAATGTTTAAggcggc                                                                     <  2:1155269/150‑1 (MQ=38)
                                                                                  cAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACgcggcg                                                                    >  2:134087/1‑150 (MQ=255)
                                                                                      gCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTa                                                                >  1:1404078/1‑150 (MQ=255)
                                                                                             gCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAActtctt                                                         >  2:294151/1‑150 (MQ=255)
                                                                                                   cTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCt                                                   >  1:1143603/1‑150 (MQ=255)
                                                                                                     gTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTa                                                 <  1:694143/150‑1 (MQ=255)
                                                                                                     gTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAAAGTTTAACGCGGCGCTTAACTTCTTCCTGCTTa                                                 >  2:305625/1‑150 (MQ=255)
                                                                                                      tCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTAc                                                >  2:1426584/1‑150 (MQ=255)
                                                                                                       cgcgTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTAcc                                               >  2:728851/1‑150 (MQ=255)
                                                                                                         cgTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAg                                             >  2:843948/1‑150 (MQ=255)
                                                                                                           tCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTa                                                                                                        >  2:327937/1‑89 (MQ=255)
                                                                                                           tCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCg                                           <  2:876181/150‑1 (MQ=255)
                                                                                                              tATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTa                                        <  1:1426584/150‑1 (MQ=255)
                                                                                                               aTCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTaa                                       >  2:499172/1‑150 (MQ=255)
                                                                                                               aTCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTaa                                       >  2:1034677/1‑150 (MQ=255)
                                                                                                               aTCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTaa                                       <  1:1140274/150‑1 (MQ=255)
                                                                                                               aTCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTaa                                       >  1:1029951/1‑150 (MQ=255)
                                                                                                                tCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTaaa                                      >  1:1380477/1‑150 (MQ=255)
                                                                                                                 cAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAAc                                     <  1:774959/150‑1 (MQ=255)
                                                                                                                 cAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAAc                                     <  1:16564/150‑1 (MQ=255)
                                                                                                                 cAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAAc                                     >  1:423944/1‑150 (MQ=255)
                                                                                                                 cAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAAc                                     >  2:508347/1‑150 (MQ=255)
                                                                                                                 cAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAAc                                     <  1:637852/150‑1 (MQ=255)
                                                                                                                  aGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACg                                    <  1:1237576/150‑1 (MQ=255)
                                                                                                                  aGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGGCTTAACTTCTTCCTGCTTACCAGCGTTAAACgg                                   <  1:429784/150‑1 (MQ=255)
                                                                                                                     ccTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGAc                                 >  2:958464/1‑150 (MQ=255)
                                                                                                                      cTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACg                                <  1:1312730/150‑1 (MQ=255)
                                                                                                                        ggAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTg                              >  1:48258/1‑150 (MQ=255)
                                                                                                                        ggAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTg                              <  1:1268715/150‑1 (MQ=255)
                                                                                                                             aaTTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATcc                         >  1:66329/1‑150 (MQ=255)
                                                                                                                               ttGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCgg                       <  1:843948/150‑1 (MQ=255)
                                                                                                                                tGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCggg                      <  1:750826/150‑1 (MQ=255)
                                                                                                                                tGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTCGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCggg                      <  1:402757/150‑1 (MQ=255)
                                                                                                                                 gCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGc                     >  2:1297823/1‑150 (MQ=255)
                                                                                                                                 cccAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGc                     >  1:325887/2‑150 (MQ=255)
                                                                                                                                  ccAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCt                    <  2:610767/150‑1 (MQ=255)
                                                                                                                                   cAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTa                   >  2:1024822/1‑150 (MQ=255)
                                                                                                                                      aTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCt                <  1:151029/150‑1 (MQ=255)
                                                                                                                                       tGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTa               <  1:563377/150‑1 (MQ=255)
                                                                                                                                        gTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTaa              >  1:1167627/1‑150 (MQ=255)
                                                                                                                                         tAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTaaa             >  2:31298/1‑150 (MQ=255)
                                                                                                                                         aatgCCGGATGCGGCGTAAACGCCTTATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTaaa             >  1:682212/4‑150 (MQ=255)
                                                                                                                                          aGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAat            >  1:1005572/1‑150 (MQ=255)
                                                                                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAata           >  2:502866/1‑150 (MQ=255)
                                                                                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAata           >  2:853627/1‑150 (MQ=255)
                                                                                                                                            gCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAatat          >  2:778728/1‑150 (MQ=255)
                                                                                                                                                    aGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAATATCCGCacac  >  2:1105865/1‑150 (MQ=255)
                                                                                                                                                    aGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAATATCCGCacac  >  2:1234086/1‑150 (MQ=255)
                                                                                                                                                    aGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAATATCCGCacac  <  2:148323/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GACGATGTCGACAGGATAGTACTGATGATAATTCATTATTTCGCCTGGCTTATCAGTCCGCGTGTAGGCTGCATCCGGTAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGACCTGGAGCAATTGCCAAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTGCTTAGCTCAGAGTGAAGATTAACCTATCTGCCCATTCCCCAAATGTTTAACGCGGCGCTTAACTTCTTCCTGCTTACCAGCGTTAAACGGACGTGCATCCGGGCTACCTAAATATCCGCACAC  >  NC_000913/4460329‑4460626

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: