Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,504,873 T→A 9.6% intergenic (+90/+62) frlR → / ← yhfS putative DNA‑binding transcriptional regulator/FNR‑regulated pyridoxal phosphate‑dependent aminotransferase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,504,8730TA 9.6% 443.9 / 11.2 161intergenic (+90/+62)frlR/yhfSputative DNA‑binding transcriptional regulator/FNR‑regulated pyridoxal phosphate‑dependent aminotransferase family protein
Reads supporting (aligned to +/- strand):  ref base T (33/104);  new base A (16/2);  total (54/107)
Fisher's exact test for biased strand distribution p-value = 1.30e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 174 total aligned reads displayed.

TCCGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCT  >  NC_000913/3504724‑3505022
                                                                                                                                                     |                                                                                                                                                      
  cGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCctct                                                                                                                                                      <  1:793593/150‑1 (MQ=255)
  cGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCGCCCCCGAACCccct                                                                                                                                                      >  1:1975187/1‑150 (MQ=255)
     aCTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTccc                                                                                                                                                   <  1:2873557/150‑1 (MQ=255)
       tcttctGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCa                                                                                                                                                 >  1:3060423/1‑150 (MQ=255)
        cttctGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTACCCTCACCCAAAACCTCAccccca                                                                                                                                                >  1:3744828/1‑148 (MQ=255)
         ttctGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCaaa                                                                                                                                               <  1:3440298/150‑1 (MQ=255)
         ttctGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCaaa                                                                                                                                               <  2:366688/150‑1 (MQ=255)
            tGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAA‑ggg                                                                                                                                           <  2:555015/150‑1 (MQ=255)
             gCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCAAAACCCAACCCTATCCCCAAA‑gtgg                                                                                                                                          >  1:3378247/1‑150 (MQ=39)
                gcgcCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCCCCCCCAAACCCACACCCCAAA‑GGGGCgc                                                                                                                                       >  1:3501112/1‑149 (MQ=39)
                 cgcCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAg                                                                                                                                      <  2:2763196/150‑1 (MQ=255)
                 cgcCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAg                                                                                                                                      <  1:380798/150‑1 (MQ=255)
                 cgcCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAg                                                                                                                                      <  1:2657179/150‑1 (MQ=255)
                  gcCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCGCCCCTTT‑GGGGAGAgg                                                                                                                                     <  1:3923508/150‑1 (MQ=39)
                   ccAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGTTGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAggg                                                                                                                                    <  1:993209/150‑1 (MQ=255)
                   ccAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAggg                                                                                                                                    <  2:489091/150‑1 (MQ=255)
                   ccAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAggg                                                                                                                                    <  1:2590247/150‑1 (MQ=255)
                   ccAATCGCATCAACTTAAATATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCTTT‑GGGGAGAggg                                                                                                                                    <  1:2062640/150‑1 (MQ=39)
                   ccAAGCGCATCACCTTAACTAGTGATAATAAAAGACCTTCACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTGTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAggg                                                                                                                                    <  2:4096001/150‑1 (MQ=39)
                     aaTCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGa                                                                                                                                  <  1:3516989/150‑1 (MQ=255)
                      ctCGCATCACCCTCCCTCGCGATAATAACAGACATTACCCGGAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGAGGCTGGCGCGTAGACTTTTACATCCCGTACGTTCCCCTCACCCTACCCCTCTCCCCAAA‑GGGGCGAGGGGAc                                                                                                                                 <  2:2886828/149‑1 (MQ=37)
                      aTCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGAc                                                                                                                                 <  1:2111022/150‑1 (MQ=255)
                           aTCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGACCGTTc                                                                                                                            <  1:1450185/150‑1 (MQ=255)
                           aTCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGACCGTTc                                                                                                                            <  1:1450293/150‑1 (MQ=255)
                              aCCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCGCACCCAAAACCACACCCCAAA‑GGGGGGGGGGGCCCggggagt                                                                                                                         >  1:1279518/1‑144 (MQ=38)
                               ccTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACACCCTAACCCACACCCCAAA‑GGGGCGAGGGGACCGGTCAGTa                                                                                                                        >  1:3710833/1‑150 (MQ=255)
                               ccTTAACTATTGATAATAAAAGACATAAACCGTAGGCCGGATAAGATGCGCCAGCATCGCAACCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTACCCACACCCAAAACCACTACCCAAACGGGGCGGGGGGCCCGtgagaa                                                                                                                         >  1:1901392/1‑145 (MQ=37)
                                 ttACCTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACg                                                                                                                      <  2:162582/150‑1 (MQ=255)
                                 ttAACTATTGATAATAAAAGACCTTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACg                                                                                                                      <  2:702419/150‑1 (MQ=255)
                                 ttAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCTCACCCTAACCCTCTCCCCAAA‑GGGGGGTGGGGTCCGTTCAGTTCg                                                                                                                      >  1:3351158/1‑150 (MQ=39)
                                 ttAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACg                                                                                                                      <  2:1558005/150‑1 (MQ=255)
                                     cTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAAt                                                                                                                  <  2:3111750/150‑1 (MQ=255)
                                     cTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAAt                                                                                                                  <  1:2503149/150‑1 (MQ=255)
                                       aTTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAAttt                                                                                                                <  2:3430398/150‑1 (MQ=255)
                                          gATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGc                                                                                                             <  1:3210666/150‑1 (MQ=255)
                                          gATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACACCCAAACCCACACCCCCAA‑CGGGCGAGGGGGCCGGGCCGTACCGAATTTAGc                                                                                                             >  1:1684019/1‑150 (MQ=37)
                                           aaaataaAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTACCCACACCCGAACCCTCTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGCAGTTAGCt                                                                                                            >  2:3392662/3‑150 (MQ=39)
                                             aataaAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGGTGAATTTTACATCCCGTACGTTGCCCACACCCAAACCCACACCCCAAA‑TGGGCGCGGGGGCCCGTCAGTACGGAAATTTTCTTg                                                                                                          >  1:1667780/1‑150 (MQ=25)
                                                 aaaGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCGCCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGCCgg                                                                                                      >  1:2098169/1‑150 (MQ=255)
                                                 aaaGACATTAACCGTAGGCCGGAGAAGATGCGCCAGCATCGCATCCGGCGATGCGGGCGCGTTGAATTTTACATCCCGTACGTTCCCCGCCCCCTAACCCGCACCCCAAA‑TGGGCGAGGGGACCGTTCAGTACGGAATGTTGCTTGACgg                                                                                                      >  2:638668/1‑150 (MQ=38)
                                                  aaGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGt                                                                                                     <  1:3294410/150‑1 (MQ=255)
                                                        ttAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTACCCACACCCTAACCCTCTACCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATc                                                                                               >  2:4066472/1‑150 (MQ=255)
                                                           aCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAg                                                                                            <  1:436836/150‑1 (MQ=255)
                                                           aCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTAAAACCCGTACGGTCCCCACACCCAAAACCACAACCCAAA‑GGGGGCAGGGGGCCCCTCACTAACGAAATTATGATGAAGGGGCTATcaac                                                                                            >  1:1023609/1‑147 (MQ=16)
                                                             cGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTGGCTTGACGGTGCTATCCAGCa                                                                                          <  1:2551114/150‑1 (MQ=255)
                                                             cGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCa                                                                                          <  2:3704190/150‑1 (MQ=255)
                                                               tAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCAtt                                                                                        <  2:3060423/150‑1 (MQ=255)
                                                                 ggCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTaa                                                                                      <  2:265812/150‑1 (MQ=255)
                                                                 ggCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTaa                                                                                      <  2:265823/150‑1 (MQ=255)
                                                                     ggATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATAc                                                                                  <  1:1570673/150‑1 (MQ=255)
                                                                      gATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGACCGGTCCGTACGGGATTTAGCGTGACGGTGCGAGCCAGCCTTAactact                                                                                 >  2:3300136/1‑150 (MQ=38)
                                                                       aTAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTg                                                                                <  2:786254/150‑1 (MQ=255)
                                                                         aaGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGc                                                                              <  2:2007568/150‑1 (MQ=255)
                                                                               cgcCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACAACCAAACACGCACCCCCAA‑GGGGCGAGGGGGCCCGTCAGTAACGAATTTATATTTTAGGTGCTATACACCATTAAATACTTGCAATAca                                                                        >  2:151194/1‑149 (MQ=255)
                                                                               cgcCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTACCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGACCCTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCAGTAAATACTGGCAATAct                                                                        >  2:2700561/1‑150 (MQ=255)
                                                                                  cAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACtctc                                                                     <  2:459280/150‑1 (MQ=255)
                                                                                    gcatcgcatcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCAc                                                                   <  1:2944890/150‑1 (MQ=255)
                                                                                    gcatcgcatcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCAc                                                                   <  2:4115479/150‑1 (MQ=255)
                                                                                    gcatcgcatcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCAc                                                                   <  1:2003294/150‑1 (MQ=255)
                                                                                     catcgcatcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCGCACCCTAACCCGCTCCCCAAA‑GGGGCGAGGGGCCCGGTCAGTACGGAATTTAGCGTGACGGTGCGCTCCCGCATTAAATACTGGCAAAACGCTCAcc                                                                  >  2:2154919/1‑149 (MQ=255)
                                                                                         gcatcCGGCGATGCTGGCGAGTTTAATTTTACTTCCCGTACGTTCCCCTCACCCTACCCCTCTCCCCTAT‑GGGGCGAGGGGACCTTTCAGTACGGAATTTAGCTTGTCGGTGCTATCCAGCATTATATACTGGCAATTCTCTCACGCaaa                                                              <  2:1984283/150‑1 (MQ=255)
                                                                                            tcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATc                                                           <  2:1360732/150‑1 (MQ=255)
                                                                                            tcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCACTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATc                                                           <  1:1009129/150‑1 (MQ=255)
                                                                                             ccGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGCTCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATcc                                                          <  1:1536197/150‑1 (MQ=255)
                                                                                              cGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCg                                                         <  1:2858976/150‑1 (MQ=255)
                                                                                               ggCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCgc                                                        <  2:1912824/150‑1 (MQ=255)
                                                                                                           gcgTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTctct                                            <  2:4042671/150‑1 (MQ=255)
                                                                                                              ttGAATTTTACATCCCGTACGGTCCCCTCACCCTCACCCGCTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCa                                         <  1:2419051/150‑1 (MQ=255)
                                                                                                               tGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCAc                                        <  1:3736538/150‑1 (MQ=255)
                                                                                                               tGAATTTTACATCCCGTACGGTCCCCGCACCCGAACCCTCAACCCAAA‑GGGGGGAGGGGTCCGGTCAGTACGGAATTTAGCTTGACGGTTCTCTCCACCATAAAATACAGGCAAAACTCTCAACCAAAAAACGCAGCAACGtctcgtaac                                        >  2:911044/1‑145 (MQ=255)
                                                                                                                gAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCAcc                                       <  2:2358487/150‑1 (MQ=255)
                                                                                                                gAATTTTACATCCCGTACGGTCCCCTCACCCTAACCCTCTCCCCCAA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCAcc                                       <  2:2358532/150‑1 (MQ=255)
                                                                                                                   ttttACATCCCGTACGTTCCCCTCACCCTAACCCTCCCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCt                                    >  1:3819572/1‑150 (MQ=255)
                                                                                                                   ttttACATCCCGTACGTTCCCCACACCCAAAACCACAACCCCAA‑CGGGCCAGGGGACCCGTCAGTACCGAAATTAGCGTGACGGTGCGAGCCACCCATAAATACAGGCAAAACGCGCACCCAAAAACCCCACCACCGGCACATCAccccc                                    >  2:144158/1‑147 (MQ=255)
                                                                                                                   atttACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCt                                    <  2:3351158/149‑1 (MQ=255)
                                                                                                                    tttACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTg                                   <  2:1294009/150‑1 (MQ=255)
                                                                                                                    tttACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTg                                   <  1:3022531/150‑1 (MQ=255)
                                                                                                                    tttACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTg                                   <  2:3378580/150‑1 (MQ=255)
                                                                                                                    tttACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTg                                   <  2:2087031/150‑1 (MQ=255)
                                                                                                                    cgtACATCCCGTACGTGCCCCTCACCCTAACCCACTCCCCACA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCCAAATCCGCAGCACCGTCTCTTCACCGCTg                                   <  2:3616028/148‑1 (MQ=255)
                                                                                                                          tCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTc                             <  1:2228136/150‑1 (MQ=255)
                                                                                                                           cccGTACGTTCCCCTCACCCAAACCCTCTACCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAAATTAGCGTGTCGGTGCTAAACAGCATTAAATACAGGCAATACTCACACGCAAAAACCGCAGCACCGTCTCTTCACCGCGGCGATTCg                            >  2:2215237/1‑150 (MQ=37)
                                                                                                                            ccGTACGTTCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCgg                           <  1:3788324/150‑1 (MQ=255)
                                                                                                                            ccGTACGGTCCCCTCCCCCTAACCCTCTCCCCCAC‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTCTCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCgg                           <  1:3726746/150‑1 (MQ=255)
                                                                                                                              gTACGTTCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGt                         >  2:1950161/1‑150 (MQ=255)
                                                                                                                                 cGTTCCCCTCACCCTAACCCTCTCCCCACA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTaa                      <  2:3763607/150‑1 (MQ=255)
                                                                                                                                 cGTTCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTaa                      <  1:869188/150‑1 (MQ=255)
                                                                                                                                   ttCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATg                    >  2:2309455/1‑150 (MQ=255)
                                                                                                                                    tCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGc                   <  2:3160675/150‑1 (MQ=255)
                                                                                                                                    tCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGc                   <  2:1264744/150‑1 (MQ=255)
                                                                                                                                        cTCACCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAAt               >  1:2474378/1‑150 (MQ=255)
                                                                                                                                           aCCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGc            >  2:434352/1‑150 (MQ=255)
                                                                                                                                           aCCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGc            >  1:2760597/1‑150 (MQ=255)
                                                                                                                                            cccTAACCCTCTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGca           <  2:659301/150‑1 (MQ=255)
                                                                                                                                                aaCCCTCTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAAt       >  2:2864108/1‑150 (MQ=255)
                                                                                                                                                  cccTCTCCCCAAC‑GGGGGGTGGGGACCGTTCAGTACGTATTTTAGCTTGACGGTGCTATCCATAATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGTTTCGAATTAACCCGAATCGctgcggtg     <  1:4066472/150‑8 (MQ=255)
                                                                                                                                                   ccTCTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCACTACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGtt    >  1:3703280/1‑150 (MQ=255)
                                                                                                                                                    ctctCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTc   >  2:1508235/1‑150 (MQ=255)
                                                                                                                                                    ctctCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTc   >  2:2102152/1‑150 (MQ=255)
                                                                                                                                                     tctcCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCt  >  2:2437464/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TCCGAACTCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAA‑GGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCT  >  NC_000913/3504724‑3505022

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