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‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGGTGTGTGGCTGTCAAC  >  NC_000913/4053422‑4053718
                                                                                                                                                        |                                                                                                                                                    
gTGAATGGCGCTTCGTTTACAAAGTAGAAAACGAAAGCTCCATCAATGTGCTGAGAACAAGATTGCCTGAAGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTA‑CGACGCATCAGGCCTC‑TGAGc                                                                                                                                                     >  2:739548/1‑150 (MQ=255)
  gAATGGCGCTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCat                                                                                                                                                   >  1:1004449/1‑150 (MQ=255)
  gAATGGCGCTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCat                                                                                                                                                   >  1:309059/1‑150 (MQ=255)
      ggCGCTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCatgatg                                                                                                                                               <  2:663136/150‑1 (MQ=255)
         gcTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCt                                                                                                                                            >  1:67179/1‑150 (MQ=255)
         gcTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCt                                                                                                                                            <  1:312177/150‑1 (MQ=255)
           ttCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCTACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGa                                                                                                                                          <  2:583096/150‑1 (MQ=255)
             cGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATg                                                                                                                                        >  2:657157/1‑150 (MQ=255)
               tttACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCTGATAAGACGCGACCAGCGCCGCATCAGGCATCATGCTCA‑GATGCCTGAt                                                                                                                                         <  2:1470658/150‑1 (MQ=21)
               tttACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTCGTAGGCCTGATAAGACGCGACCAGCGCCGCATCAGGCATCATGCTCA‑GATGCCTGAt                                                                                                                                         <  1:1131914/150‑1 (MQ=18)
                ttACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATCATGCTCA‑GATGTCTGATgcgg                                                                                                                                     <  1:71709/150‑1 (MQ=38)
                taaCAAAGTAGAAAAAGAAAGCGCCATCAATGTGCTGGGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATCATGCTCA‑GATGCCTGATgcgg                                                                                                                                     <  2:1721254/148‑1 (MQ=38)
                 tACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCTGATAAGACGCGACCAGCGCCGCATCAGGCATCATGCTCA‑GATGCCTGATgc                                                                                                                                       <  1:993269/150‑1 (MQ=21)
                 tACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATgcggc                                                                                                                                    <  1:594619/150‑1 (MQ=255)
                  aCAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATgcggcg                                                                                                                                   <  1:135497/150‑1 (MQ=255)
                   cAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGgcgc                                                                                                                                  >  2:54243/1‑150 (MQ=255)
                    aaaGTAGAAAACGAAAGCGACATCAATGAGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATAATTGCTTTTGTAGGCAGGATAAG‑‑TCGTTTA‑CGCCGCATCAGGCAGC‑TCAGCATGATGCCATATGCCGCGCt                                                                                                                                 >  1:832529/1‑150 (MQ=255)
                     aaGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTg                                                                                                                                >  1:1105558/1‑150 (MQ=255)
                          gAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTcgc                                                                                                                           >  2:1446355/1‑150 (MQ=255)
                            aaaCGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGt                                                                                                                         >  2:1046824/1‑150 (MQ=255)
                            aaaCGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGt                                                                                                                         >  2:227879/1‑150 (MQ=255)
                            aaaCGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTA‑CCCCGAAACAGGAATA‑TGAGAATTATGAATGATGCGGCGCTGTACGCGt                                                                                                                         >  1:1396620/1‑150 (MQ=255)
                              aCGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCt                                                                                                                       <  1:187254/150‑1 (MQ=255)
                               cGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCtt                                                                                                                      <  2:1253102/150‑1 (MQ=255)
                                 aaaGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGAGAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTAt                                                                                                                    >  2:901498/1‑150 (MQ=255)
                                  aaGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACGTGTTTAAATTTTTGCTTTTGTTGGCCGGATACG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATc                                                                                                                   <  1:603445/150‑1 (MQ=255)
                                      gcCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTCGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCCTGCTGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGc                                                                                                               <  2:521485/150‑1 (MQ=255)
                                      gcCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGc                                                                                                               <  2:1015741/150‑1 (MQ=255)
                                        cATCAATGTGCTGATAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCt                                                                                                             <  2:1654137/150‑1 (MQ=255)
                                        cATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACTAAACCTGTTGATATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCTTGATGCCTGATGCGGCGCTGGTCGAGTCTTATCAGGCCt                                                                                                             >  2:1388997/1‑150 (MQ=255)
                                        cATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACTAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATTATGCCTGATGCTGCGATGTTCTCGTCTTATCAgtact                                                                                                             >  2:942600/1‑146 (MQ=38)
                                        cATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGTATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑AGAGCATGATGCCTGATGAGTCGCTGGTCGCGTCTTATCAGGCCt                                                                                                             >  2:77467/1‑150 (MQ=255)
                                        cATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCt                                                                                                             >  2:1562017/1‑150 (MQ=255)
                                        cATCAATGTGCTGAGAACAAGATTGCCTGAGGCGCTGGGCAACAAAACCTGTTGAAATTAATGCTTTTGTAGGCCGTATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCt                                                                                                             <  2:37359/150‑1 (MQ=255)
                                            aaTGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGa                                                                                                         >  1:766461/1‑150 (MQ=255)
                                              tgtgCTGAGAAAACGATTGCCTGATGCGCTGGGCAACAAAAGCTGTGTAAATTATTGCTTTTGTAGGCCGGATACG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACt                                                                                                       <  2:1340604/150‑1 (MQ=37)
                                               gtgCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACtt                                                                                                      <  1:1150602/150‑1 (MQ=255)
                                                tgCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCTGATAAGACGCGACCAGCGCCGCATCAGGCATCATGCTCA‑GATGCCTGATg                                                                                                                                        >  1:1268596/1‑118 (MQ=9)
                                                tgCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCTGATAAGACGCGACCAGCGCCGCATCAGGCATCATGCTCA‑GATGCCTGATg                                                                                                                                        <  2:1268596/118‑1 (MQ=9)
                                                  cTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTtct                                                                                                   >  2:1179717/1‑150 (MQ=255)
                                                  cTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCACC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTtct                                                                                                   >  2:312321/1‑150 (MQ=255)
                                                   tGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTtctc                                                                                                  <  1:502027/150‑1 (MQ=255)
                                                      gaACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGAGAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGAATGATGCCTGATGCGGCGCTGGTCGCGGCTTTTCAGGCCTACGACTTTCTcccc                                                                                               >  2:1062767/1‑148 (MQ=255)
                                                        aCAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGATCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCTGGCCTACGACTTTCGCCGCAt                                                                                             >  2:265912/1‑150 (MQ=255)
                                                         cAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCAACGACTTTCTCCGCATc                                                                                            >  1:1606137/1‑150 (MQ=255)
                                                           aGATTGCCTGATGCGCTGGGCAACAAAACCTGTTTAAATTATTGCTTTTGTCGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGAATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCg                                                                                          <  2:1419592/150‑1 (MQ=255)
                                                           aGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCg                                                                                          >  1:1155620/1‑150 (MQ=255)
                                                             tttGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTA‑CGCGGCATCAGGAATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGc                                                                                        <  2:1313145/149‑1 (MQ=255)
                                                             aTTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGc                                                                                        <  2:1263096/150‑1 (MQ=255)
                                                                    gATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaaca                                                                                 >  1:88091/1‑150 (MQ=255)
                                                                    gATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAaaca                                                                                 <  1:400952/150‑1 (MQ=255)
                                                                          cTGGGCAACAAAACCTTTTTAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTAA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGc                                                                           <  2:792630/150‑1 (MQ=255)
                                                                             ggCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCTGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTc                                                                        >  1:1063506/1‑150 (MQ=255)
                                                                               caacaaAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAg                                                                      <  2:1105558/150‑1 (MQ=255)
                                                                               caacaaAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAg                                                                      >  2:769196/1‑150 (MQ=255)
                                                                               caacaaAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAg                                                                      >  2:1489786/1‑150 (MQ=255)
                                                                                aacaaAACCTGATTAAATTATTGCTTTTGTCGGCCGGATAAG‑‑GCGTTTC‑CCCCGCATCAGGCATC‑TGAGCATGATGCCGGTTGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGc                                                                     <  2:335981/150‑1 (MQ=255)
                                                                                 acatAACCTGTTGAATTTATTGCGTTCGTAGTCCTGATTAG‑‑GCGTTTA‑CGCAGCATCAGGCATA‑TGAGCTTGTTGCCTGATGCGGCGCGGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCt                                                                    <  2:830953/150‑1 (MQ=255)
                                                                                 acaaAACCTGTTGAAATTTTTGCTTTCGTAGGCCGGAAAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCTTGATGCCGGATGCGGCGGTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCt                                                                    <  2:825382/150‑1 (MQ=255)
                                                                                 acaaAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCt                                                                    <  2:240265/150‑1 (MQ=255)
                                                                                  caaAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCTCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTa                                                                   <  2:761611/150‑1 (MQ=255)
                                                                                        cTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGc                                                             <  2:1006962/150‑1 (MQ=255)
                                                                                         tGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGcc                                                            <  1:661722/150‑1 (MQ=255)
                                                                                          gTTTAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGccc                                                           <  1:1374917/150‑1 (MQ=255)
                                                                                             taaaTTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAAt                                                        <  1:1001927/149‑1 (MQ=255)
                                                                                              aaaTTTTTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCGTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAACCATCCCAATc                                                       >  2:831485/1‑150 (MQ=255)
                                                                                                    ttGCTTTTGTAGGCCGGATTAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCg                                                 <  1:1037874/150‑1 (MQ=255)
                                                                                                    ttGCTTTTGTAGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCg                                                 >  2:1242975/1‑150 (MQ=255)
                                                                                                     tGCTTTTGTAGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCgg                                                >  1:1582039/1‑150 (MQ=255)
                                                                                                        ttttGTAGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGcac                                             <  2:1004449/150‑1 (MQ=255)
                                                                                                            gTAGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCCCGTCTTATCATGCATACGACTTTCTCTGCATCCGGCATAAACAAAGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAg                                         >  2:918790/1‑150 (MQ=38)
                                                                                                              aGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACCAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGcc                                       >  2:831690/1‑149 (MQ=255)
                                                                                                               ggCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCAc                                      >  1:1695998/1‑150 (MQ=255)
                                                                                                               ggCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCAc                                      >  2:509783/1‑150 (MQ=255)
                                                                                                               ggCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATAATAGCGGCACACAGCAc                                      >  1:616106/1‑150 (MQ=255)
                                                                                                               ggCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTACTCCGCATCCGGCATAAACAACGAGCTTCAGCTCAACAGCCCAATCATTGCGGCACACCGCAc                                      >  1:259916/1‑150 (MQ=255)
                                                                                                                gCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACg                                     >  1:497603/1‑150 (MQ=255)
                                                                                                                 ccGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGCCTTTCTCCGCCTCCGGCATAAACGACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACCCAGCACgg                                    >  1:253454/1‑150 (MQ=255)
                                                                                                                 ccGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACgg                                    >  1:317589/1‑150 (MQ=255)
                                                                                                                  cGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGgc                                   <  2:1371653/150‑1 (MQ=255)
                                                                                                                    gATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGgcag                                 <  2:1325342/150‑1 (MQ=255)
                                                                                                                             cGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGc                          <  2:1129374/150‑1 (MQ=255)
                                                                                                                               tttA‑CGCCGCATCAGGCAGC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAg                        <  2:1119170/150‑1 (MQ=255)
                                                                                                                                ttA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGt                       <  1:712622/150‑1 (MQ=255)
                                                                                                                                     gccgcATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGc                   >  2:858962/1‑150 (MQ=255)
                                                                                                                                     gccgcATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGc                   <  1:224314/150‑1 (MQ=255)
                                                                                                                                         cATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGgtg               >  2:726840/1‑150 (MQ=255)
                                                                                                                                            cAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTCTCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAACCAGCCCAATCATTGCGGCACCCCGCACGGCAGCAACGGCAGTTTGCGgtgtgt            >  1:121273/1‑150 (MQ=255)
                                                                                                                                            cAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGgtgtgt            >  2:954427/1‑150 (MQ=255)
                                                                                                                                             aGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGGtgtgtg           >  1:405386/1‑150 (MQ=255)
                                                                                                                                              ggCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGGTGTGTgg          <  1:198834/150‑1 (MQ=255)
                                                                                                                                               gCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGGTGTGTGGc         <  1:769196/150‑1 (MQ=255)
                                                                                                                                               gCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGGTGTGTGGc         >  2:1673679/1‑150 (MQ=255)
                                                                                                                                               gCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGGTGTGTGGc         <  2:405386/150‑1 (MQ=255)
                                                                                                                                               gCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATAGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGGTGTGTGGc         >  1:1373607/1‑150 (MQ=255)
                                                                                                                                                cATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGGTGTGTGGCt        <  1:433782/150‑1 (MQ=255)
                                                                                                                                                 aTC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGGTGTGTGGCTg       <  2:1490000/150‑1 (MQ=255)
                                                                                                                                                   c‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGGTGTGTGGCTGTc     >  2:732360/1‑150 (MQ=255)
                                                                                                                                                     tGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGTAACGGCAGTTTGCGGTGTGTGGCTGTCa    >  2:202409/1‑150 (MQ=255)
                                                                                                                                                      gAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGGTGTGTGGCTGTCaa   <  1:54243/150‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                                    
GTGAATGGCGCTTCGTTTACAAAGTAGAAAACGAAAGCGCCATCAATGTGCTGAGAACAAGATTGCCTGATGCGCTGGGCAACAAAACCTGTTGAAATTATTGCTTTTGTAGGCCGGATAAG‑‑GCGTTTA‑CGCCGCATCAGGCATC‑TGAGCATGATGCCTGATGCGGCGCTGGTCGCGTCTTATCAGGCCTACGACTTTCTCCGCATCCGGCATAAACAACGAGCTTCAGCTAAACAGCCCAATCATTGCGGCACACAGCACGGCAGCAACGGCAGTTTGCGGTGTGTGGCTGTCAAC  >  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: ‑