Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,053,570 G→T 6.9% intergenic (+152/+77) hemN → / ← yshB coproporphyrinogen III oxidase, SAM and NAD(P)H dependent, oxygen‑independent/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,053,5700GT 6.9% 320.1 / 9.7 112intergenic (+152/+77)hemN/yshBcoproporphyrinogen III oxidase, SAM and NAD(P)H dependent, oxygen‑independent/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (53/51);  new base T (2/6);  total (55/57)
Fisher's exact test for biased strand distribution p-value = 2.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.15e-01

Only 100 of 112 total aligned reads displayed.

GTGAATGGCGCTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGGTGTGTGGCTGTCAAC  >  NC_000913/4053422‑4053718
                                                                                                                                                        |                                                                                                                                                    
gTGAATGGCGCTTCGTTTACAAAGTAGAAAACGAAAGCTCCATCAATGTGCTGAGAACAAGATTGCCTGAAGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACGACGCATCAGGCCTCTGAGc                                                                                                                                                     >  2:739548/1‑150 (MQ=255)
  gAATGGCGCTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCat                                                                                                                                                   >  1:1004449/1‑150 (MQ=255)
  gAATGGCGCTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCat                                                                                                                                                   >  1:309059/1‑150 (MQ=255)
      ggCGCTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCatgatg                                                                                                                                               <  2:663136/150‑1 (MQ=255)
         gcTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCt                                                                                                                                            >  1:67179/1‑150 (MQ=255)
         gcTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCt                                                                                                                                            <  1:312177/150‑1 (MQ=255)
           ttCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCTACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGa                                                                                                                                          <  2:583096/150‑1 (MQ=255)
             cGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATg                                                                                                                                        >  2:657157/1‑150 (MQ=255)
               tttACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCTGATAAGACGCGACCAGCGCCGCATCAGGCATCATGCTCAGATGCCTGAt                                                                                                                                         <  2:1470658/150‑1 (MQ=21)
               tttACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTCGTAGGCCTGATAAGACGCGACCAGCGCCGCATCAGGCATCATGCTCAGATGCCTGAt                                                                                                                                         <  1:1131914/150‑1 (MQ=18)
                ttACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCATGCTCAGATGTCTGATgcgg                                                                                                                                     <  1:71709/150‑1 (MQ=38)
                taaCAAAGTAGAAAAAGAAAGCGCCATCAATGTGCTGGGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCATGCTCAGATGCCTGATgcgg                                                                                                                                     <  2:1721254/148‑1 (MQ=38)
                 tACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCTGATAAGACGCGACCAGCGCCGCATCAGGCATCATGCTCAGATGCCTGATgc                                                                                                                                       <  1:993269/150‑1 (MQ=21)
                 tACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATgcggc                                                                                                                                    <  1:594619/150‑1 (MQ=255)
                  aCAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATgcggcg                                                                                                                                   <  1:135497/150‑1 (MQ=255)
                   cAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGgcgc                                                                                                                                  >  2:54243/1‑150 (MQ=255)
                    aaaGTAGAAAACGAAAGCGACATCAATGAGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATAATTGCTTTTGTAGGCAGGATAAG‑‑TCGTTTACGCCGCATCAGGCAGCTCAGCATGATGCCATATGCCGCGCt                                                                                                                                 >  1:832529/1‑150 (MQ=255)
                     aaGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTg                                                                                                                                >  1:1105558/1‑150 (MQ=255)
                          gAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTcgc                                                                                                                           >  2:1446355/1‑150 (MQ=255)
                            aaaCGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGt                                                                                                                         >  2:1046824/1‑150 (MQ=255)
                            aaaCGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGt                                                                                                                         >  2:227879/1‑150 (MQ=255)
                            aaaCGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACCCCGAAACAGGAATATGAGAATTATGAATGATGCGGCGCTGTACGCGt                                                                                                                         >  1:1396620/1‑150 (MQ=255)
                              aCGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCt                                                                                                                       <  1:187254/150‑1 (MQ=255)
                               cGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCtt                                                                                                                      <  2:1253102/150‑1 (MQ=255)
                                 aaaGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGAGAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTAt                                                                                                                    >  2:901498/1‑150 (MQ=255)
                                  aaGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACGTGTTTAAATTTTTGCTTTTGTTGGCCGGATACG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATc                                                                                                                   <  1:603445/150‑1 (MQ=255)
                                      gcCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTCGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCCTGCTGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGc                                                                                                               <  2:521485/150‑1 (MQ=255)
                                      gcCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGc                                                                                                               <  2:1015741/150‑1 (MQ=255)
                                        cATCAATGTGCTGATAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCt                                                                                                             <  2:1654137/150‑1 (MQ=255)
                                        cATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACTAAACCTGTTGATATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCTTGATGCCTGATGCGGCGCTGGTCGAGTCTTATCAGGCCt                                                                                                             >  2:1388997/1‑150 (MQ=255)
                                        cATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACTAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATTATGCCTGATGCTGCGATGTTCTCGTCTTATCAgtact                                                                                                             >  2:942600/1‑146 (MQ=38)
                                        cATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGTATAAG‑‑GCGTTTACGCCGCATCAGGCATCAGAGCATGATGCCTGATGAGTCGCTGGTCGCGTCTTATCAGGCCt                                                                                                             >  2:77467/1‑150 (MQ=255)
                                        cATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCt                                                                                                             >  2:1562017/1‑150 (MQ=255)
                                        cATCAATGTGCTGAGAACAAGATTGCCTGAGGCGCTGGGCAACAAAACCTGTTGAAATTAATGCTTTTGTAGGCCGTATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCt                                                                                                             <  2:37359/150‑1 (MQ=255)
                                            aaTGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGa                                                                                                         >  1:766461/1‑150 (MQ=255)
                                              tgtgCTGAGAAAACGATTGCCTGATGCGCTGGGCAACAAAAGCTGTGTAAATTATTGCTTTTGTAGGCCGGATACG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACt                                                                                                       <  2:1340604/150‑1 (MQ=37)
                                               gtgCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACtt                                                                                                      <  1:1150602/150‑1 (MQ=255)
                                                tgCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCTGATAAGACGCGACCAGCGCCGCATCAGGCATCATGCTCAGATGCCTGATg                                                                                                                                        >  1:1268596/1‑118 (MQ=9)
                                                tgCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCTGATAAGACGCGACCAGCGCCGCATCAGGCATCATGCTCAGATGCCTGATg                                                                                                                                        <  2:1268596/118‑1 (MQ=9)
                                                  cTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTtct                                                                                                   >  2:1179717/1‑150 (MQ=255)
                                                  cTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCACCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTtct                                                                                                   >  2:312321/1‑150 (MQ=255)
                                                   tGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTtctc                                                                                                  <  1:502027/150‑1 (MQ=255)
                                                      gaACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGAGAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGAATGATGCCTGATGCGGCGCTGGTCGCGGCTTTTCAGGCCTACGACTTTCTcccc                                                                                               >  2:1062767/1‑148 (MQ=255)
                                                        aCAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGATCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCTGGCCTACGACTTTCGCCGCAt                                                                                             >  2:265912/1‑150 (MQ=255)
                                                         cAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCAACGACTTTCTCCGCATc                                                                                            >  1:1606137/1‑150 (MQ=255)
                                                           aGATTGCCTGATGCGCTGGGCAACAAAACCTGTTTAAATTATTGCTTTTGTCGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGAATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCg                                                                                          <  2:1419592/150‑1 (MQ=255)
                                                           aGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCg                                                                                          >  1:1155620/1‑150 (MQ=255)
                                                             tttGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACGCGGCATCAGGAATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGc                                                                                        <  2:1313145/149‑1 (MQ=255)
                                                             aTTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGc                                                                                        <  2:1263096/150‑1 (MQ=255)
                                                                    gATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaaca                                                                                 >  1:88091/1‑150 (MQ=255)
                                                                    gATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaaca                                                                                 <  1:400952/150‑1 (MQ=255)
                                                                          cTGGGCAACAAAACCTTTTTAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTAACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGc                                                                           <  2:792630/150‑1 (MQ=255)
                                                                             ggCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCTGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTc                                                                        >  1:1063506/1‑150 (MQ=255)
                                                                               caacaaAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAg                                                                      <  2:1105558/150‑1 (MQ=255)
                                                                               caacaaAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAg                                                                      >  2:769196/1‑150 (MQ=255)
                                                                               caacaaAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAg                                                                      >  2:1489786/1‑150 (MQ=255)
                                                                                aacaaAACCTGATTAAATTATTGCTTTTGTCGGCCGGATAAG‑‑GCGTTTCCCCCGCATCAGGCATCTGAGCATGATGCCGGTTGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGc                                                                     <  2:335981/150‑1 (MQ=255)
                                                                                 acatAACCTGTTGAATTTATTGCGTTCGTAGTCCTGATTAG‑‑GCGTTTACGCAGCATCAGGCATATGAGCTTGTTGCCTGATGCGGCGCGGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCt                                                                    <  2:830953/150‑1 (MQ=255)
                                                                                 acaaAACCTGTTGAAATTTTTGCTTTCGTAGGCCGGAAAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCTTGATGCCGGATGCGGCGGTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCt                                                                    <  2:825382/150‑1 (MQ=255)
                                                                                 acaaAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCt                                                                    <  2:240265/150‑1 (MQ=255)
                                                                                  caaAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCTCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTa                                                                   <  2:761611/150‑1 (MQ=255)
                                                                                        cTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGc                                                             <  2:1006962/150‑1 (MQ=255)
                                                                                         tGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGcc                                                            <  1:661722/150‑1 (MQ=255)
                                                                                          gTTTAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGccc                                                           <  1:1374917/150‑1 (MQ=255)
                                                                                             taaaTTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAAt                                                        <  1:1001927/149‑1 (MQ=255)
                                                                                              aaaTTTTTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCGTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAACCATCCCAATc                                                       >  2:831485/1‑150 (MQ=255)
                                                                                                    ttGCTTTTGTAGGCCGGATTAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCg                                                 <  1:1037874/150‑1 (MQ=255)
                                                                                                    ttGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCg                                                 >  2:1242975/1‑150 (MQ=255)
                                                                                                     tGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCgg                                                >  1:1582039/1‑150 (MQ=255)
                                                                                                        ttttGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGcac                                             <  2:1004449/150‑1 (MQ=255)
                                                                                                            gTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCCCGTCTTATCATGCATACGACTTTCTCTGCATCCGGCATAAACAAAGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAg                                         >  2:918790/1‑150 (MQ=38)
                                                                                                              aGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACCAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGcc                                       >  2:831690/1‑149 (MQ=255)
                                                                                                               ggCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCAc                                      >  1:1695998/1‑150 (MQ=255)
                                                                                                               ggCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCAc                                      >  2:509783/1‑150 (MQ=255)
                                                                                                               ggCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATAATAGCGGCACACAGCAc                                      >  1:616106/1‑150 (MQ=255)
                                                                                                               ggCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTACTCCGCATCCGGCATAAACAACGAGCTTCAGCTCAACAGCCCAATCATTGCGGCACACCGCAc                                      >  1:259916/1‑150 (MQ=255)
                                                                                                                gCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACg                                     >  1:497603/1‑150 (MQ=255)
                                                                                                                 ccGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGCCTTTCTCCGCCTCCGGCATAAACGACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACCCAGCACgg                                    >  1:253454/1‑150 (MQ=255)
                                                                                                                 ccGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACgg                                    >  1:317589/1‑150 (MQ=255)
                                                                                                                  cGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGgc                                   <  2:1371653/150‑1 (MQ=255)
                                                                                                                    gATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGgcag                                 <  2:1325342/150‑1 (MQ=255)
                                                                                                                             cGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGc                          <  2:1129374/150‑1 (MQ=255)
                                                                                                                               tttACGCCGCATCAGGCAGCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAg                        <  2:1119170/150‑1 (MQ=255)
                                                                                                                                ttACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGt                       <  1:712622/150‑1 (MQ=255)
                                                                                                                                     gccgcATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGc                   >  2:858962/1‑150 (MQ=255)
                                                                                                                                     gccgcATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGc                   <  1:224314/150‑1 (MQ=255)
                                                                                                                                         cATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGgtg               >  2:726840/1‑150 (MQ=255)
                                                                                                                                            cAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTCTCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAACCAGCCCAATCATTGCGGCACCCCGCACGGCAGCAACGGCAGTTTGCGgtgtgt            >  1:121273/1‑150 (MQ=255)
                                                                                                                                            cAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGgtgtgt            >  2:954427/1‑150 (MQ=255)
                                                                                                                                             aGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGGtgtgtg           >  1:405386/1‑150 (MQ=255)
                                                                                                                                              ggCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGGTGTGTgg          <  1:198834/150‑1 (MQ=255)
                                                                                                                                               gCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGGTGTGTGGc         <  1:769196/150‑1 (MQ=255)
                                                                                                                                               gCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGGTGTGTGGc         >  2:1673679/1‑150 (MQ=255)
                                                                                                                                               gCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGGTGTGTGGc         <  2:405386/150‑1 (MQ=255)
                                                                                                                                               gCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATAGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGGTGTGTGGc         >  1:1373607/1‑150 (MQ=255)
                                                                                                                                                cATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGGTGTGTGGCt        <  1:433782/150‑1 (MQ=255)
                                                                                                                                                 aTCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGGTGTGTGGCTg       <  2:1490000/150‑1 (MQ=255)
                                                                                                                                                   cTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGGTGTGTGGCTGTc     >  2:732360/1‑150 (MQ=255)
                                                                                                                                                     tGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGTAACGGCAGTTTGCGGTGTGTGGCTGTCa    >  2:202409/1‑150 (MQ=255)
                                                                                                                                                      gAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGGTGTGTGGCTGTCaa   <  1:54243/150‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                                    
GTGAATGGCGCTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCAGGCATCTGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGGTGTGTGGCTGTCAAC  >  NC_000913/4053422‑4053718

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