Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 431,107 A→T 6.1% intergenic (+155/+22) yajD → / ← tsx HNH nuclease family protein/nucleoside channel, receptor of phage T6 and colicin K

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913431,1070AT 6.1% 360.0 / 6.9 123intergenic (+155/+22)yajD/tsxHNH nuclease family protein/nucleoside channel, receptor of phage T6 and colicin K
Reads supporting (aligned to +/- strand):  ref base A (65/50);  new base T (2/6);  total (67/56)
Fisher's exact test for biased strand distribution p-value = 1.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 124 total aligned reads displayed.

CGTAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAGTAACGAG  >  NC_000913/430959‑431250
                                                                                                                                                    |                                                                                                                                               
cGTAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGaa                                                                                                                                                 >  2:429251/1‑149 (MQ=255)
 gTAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAAt                                                                                                                                                >  1:1009838/1‑149 (MQ=255)
  tAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATa                                                                                                                                               >  1:1048461/1‑149 (MQ=255)
  tAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATa                                                                                                                                               >  1:1392808/1‑149 (MQ=255)
       ttGCCTGATGCGCTACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATcc                                                                                                                                          <  2:907302/149‑1 (MQ=255)
       ttGCCTGATGCGCTACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATcc                                                                                                                                          >  1:940995/1‑149 (MQ=255)
         gCCTGATGCGCTACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCg                                                                                                                                        <  2:1150229/149‑1 (MQ=255)
         gCCTGATGCGCTACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCg                                                                                                                                        <  2:164077/149‑1 (MQ=255)
            tGATGCGCTACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGATTTATTCCCGGCa                                                                                                                                     <  1:1011700/149‑1 (MQ=255)
            tGATGCGCTACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGATTTATTCCCGGCa                                                                                                                                     <  1:665016/149‑1 (MQ=255)
             gATGCGCTACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTAGGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGATTTATTCCCGGCAt                                                                                                                                    <  1:875255/149‑1 (MQ=255)
             gATGCGCTACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCAt                                                                                                                                    >  1:565738/1‑149 (MQ=255)
              aTGCGCTACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCAtt                                                                                                                                   >  2:467321/1‑149 (MQ=255)
                gcgcTACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGATTTATTCCCGGCAtttc                                                                                                                                 <  2:1222515/149‑2 (MQ=255)
                gcgcTACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCAtttt                                                                                                                                 <  2:213398/149‑1 (MQ=255)
                 cgcTACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGATTTATTCCCGGCAtttct                                                                                                                                <  1:748054/149‑3 (MQ=255)
                  gcTACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCa                                                                                                                               <  2:270221/149‑1 (MQ=255)
                      cGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAAt                                                                                                                           >  1:1116704/1‑149 (MQ=255)
                           aTCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTATCAATCAGAAATGCCGGGAATAAGTCCCGGAATTTACATAATCAGaa                                                                                                                      >  2:121270/1‑149 (MQ=255)
                           aTCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGaa                                                                                                                      <  1:287338/149‑1 (MQ=255)
                            tCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAg                                                                                                                     <  1:658503/149‑1 (MQ=255)
                              aGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGtt                                                                                                                   >  1:541087/1‑149 (MQ=255)
                              aGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGtt                                                                                                                   >  1:1579702/1‑149 (MQ=255)
                                  cTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTaa                                                                                                               >  2:1426954/1‑149 (MQ=255)
                                  cTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTaa                                                                                                               >  1:1056955/1‑149 (MQ=255)
                                   tACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAAc                                                                                                              <  1:238950/149‑1 (MQ=255)
                                   tACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAAc                                                                                                              >  2:1214508/1‑149 (MQ=255)
                                      gTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGATTAAATCCCGGGATTTTCATAATCAGAAGTTGTAACcta                                                                                                           <  1:414050/149‑1 (MQ=255)
                                      gTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACcta                                                                                                           >  1:344716/1‑149 (MQ=255)
                                       ttATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACctac                                                                                                          <  1:936383/149‑1 (MQ=255)
                                          tttCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCtactac                                                                                                       <  1:911767/149‑1 (MQ=255)
                                             cAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAg                                                                                                    <  2:1579702/149‑1 (MQ=255)
                                             cAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAg                                                                                                    <  2:868362/149‑1 (MQ=255)
                                             cAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAg                                                                                                    <  1:224940/149‑1 (MQ=255)
                                             aagaaaTATATTGAATTTACGTGCTTTTGTTGGCGGGAGAAAGTGTTTACGCCTCATCATGCATGAACAATGCACACGTTGTTAACATTCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACAAg                                                                                                    <  2:376741/145‑1 (MQ=255)
                                              aGCAATATATTGATTTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACATAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCAGTTTGATAATAAGAAGTTGTAACCTACTACCAgg                                                                                                   >  1:118664/1‑149 (MQ=38)
                                              aGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAgg                                                                                                   >  2:1254286/1‑149 (MQ=255)
                                                cAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTa                                                                                                 >  1:545130/1‑149 (MQ=255)
                                                cAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTa                                                                                                 >  2:100006/1‑149 (MQ=255)
                                                 aaTATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGTAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTaa                                                                                                <  1:624870/149‑1 (MQ=255)
                                                     tatTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAaccac                                                                                            <  2:1195873/149‑1 (MQ=255)
                                                     tatTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAaccac                                                                                            <  2:1392808/149‑1 (MQ=255)
                                                      atTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAaccacc                                                                                           >  2:1461812/1‑149 (MQ=255)
                                                      atTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAaccacc                                                                                           <  1:494501/149‑1 (MQ=255)
                                                       ttGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCAccc                                                                                          >  2:1311771/1‑149 (MQ=255)
                                                          aaTTTACGTGCTTTTGTTGGCCGGACTAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCACAg                                                                                       >  1:988038/1‑149 (MQ=255)
                                                          aaTTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAg                                                                                       >  2:1429844/1‑149 (MQ=255)
                                                          aaTTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAAAAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAg                                                                                       >  1:1186785/1‑149 (MQ=255)
                                                               aCGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGt                                                                                  >  2:439011/1‑149 (MQ=255)
                                                                 gTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTag                                                                                >  2:489688/1‑149 (MQ=255)
                                                                  tGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTaga                                                                               <  2:376472/149‑1 (MQ=255)
                                                                     ttttGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCTCACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTGCCAGGTAACCACCCCAGCCGGTAGAGCg                                                                            >  1:209297/1‑149 (MQ=255)
                                                                         gttgGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACg                                                                        >  1:1267706/1‑149 (MQ=255)
                                                                              ccGGACAAATCGTTTAAGGCGGATCCGGCATGTACAAATCAAAAGTTGTTAAAAATAATAAATGGCGGGTATAAATCCCGGCATTTTAATAATCCCCATTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGaa                                                                   <  2:1020699/149‑1 (MQ=255)
                                                                              ccGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGaa                                                                   >  2:1340769/1‑149 (MQ=255)
                                                                              ccGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGAGAACGTTGaa                                                                   >  1:750888/1‑149 (MQ=255)
                                                                              ccGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAAACTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGaa                                                                   >  1:1079418/1‑149 (MQ=255)
                                                                                 gACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAgtt                                                                >  2:654205/1‑149 (MQ=255)
                                                                                    aaaGCGTTTACGCCGCATCCGTCATGAACAAAGCACACTTTGTTAACAATCAGAAATGCCGTGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCTGTAGAGCGAACGTTGAAgttgcc                                                             <  2:601505/149‑1 (MQ=255)
                                                                                    aaaGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGTTAGATCGAACGTTGAAgttgcc                                                             >  2:334459/1‑149 (MQ=255)
                                                                                    aaaGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGATCGTTGAAgttgcc                                                             >  1:333846/1‑149 (MQ=255)
                                                                                    aaaGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAgttgcc                                                             >  2:579868/1‑149 (MQ=255)
                                                                                     aaGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAgttgccg                                                            >  1:1196324/1‑149 (MQ=255)
                                                                                        cGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAgttgccgttg                                                         <  2:1048461/149‑1 (MQ=255)
                                                                                          tttACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAgttgccgttgcc                                                       >  1:1102410/1‑149 (MQ=255)
                                                                                                ccgcATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTATCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAagtt                                                 >  2:470763/1‑149 (MQ=255)
                                                                                                    aTCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAagttcagt                                             <  2:1176307/149‑1 (MQ=255)
                                                                                                     tCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAagttcagtt                                            >  2:4789/1‑149 (MQ=255)
                                                                                                     tCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCAGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCATCCGGAAGAGCGAAAGATGACGTTGCCGTTGCCGAagtacagtt                                            >  1:1150225/1‑149 (MQ=255)
                                                                                                       cGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCt                                          <  2:202408/149‑1 (MQ=255)
                                                                                                         gCATGAACAAAGCACAGGTTGTTAACAATCAGAAATGGCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGc                                        <  2:757425/149‑1 (MQ=255)
                                                                                                                 aaaGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATcgtcgt                                >  1:1584859/1‑149 (MQ=255)
                                                                                                                   aGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTc                              <  1:1311771/149‑1 (MQ=255)
                                                                                                                   aGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTc                              >  1:1342359/1‑149 (MQ=255)
                                                                                                                    gCACACGTTGTTAACTATCATAAATGCTGGGAATAAATGCAGGCATTTTCATAATGAGAAGTTGTAACCTACTTCCTTGGAACCTCCCCATCCGGTTGAGCGAACGTTGAAGTTGCCGTTGCGGATGGTCTTTTTTGTATGGTCGTTcc                             <  2:1079418/149‑1 (MQ=255)
                                                                                                                      acacGTTGTTAAGAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGCTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCAc                           <  1:262762/149‑1 (MQ=255)
                                                                                                                      acacGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCAc                           <  1:429251/149‑1 (MQ=255)
                                                                                                                      acacGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCAc                           <  1:1099108/149‑1 (MQ=255)
                                                                                                                       cacGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCTGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCAGCATCGTCGTTCCACt                          >  2:697815/1‑149 (MQ=255)
                                                                                                                       cacGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACt                          >  1:684431/1‑149 (MQ=255)
                                                                                                                        acGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTg                         <  2:100225/149‑1 (MQ=255)
                                                                                                                        acGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTg                         >  1:564503/1‑149 (MQ=255)
                                                                                                                        acGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTg                         >  1:298396/1‑149 (MQ=255)
                                                                                                                        acGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCAGAAGTTCAGTTCTCCATCGTCGTTCCACTg                         >  1:810831/1‑149 (MQ=255)
                                                                                                                         cGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGGTGTAACCTACTACCAGGTAACCACCCCAGCCGGGAGAGCGAACGTTGAAGTCGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGa                        <  1:1037819/149‑1 (MQ=255)
                                                                                                                               tAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACcg                  >  1:757492/1‑149 (MQ=255)
                                                                                                                                aacaaacaGAAATGCCGGGAATAAAACCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGGACAGTTCTGCATCCTCGATCCACTGACCACcgt                 >  1:174495/1‑149 (MQ=255)
                                                                                                                                    aTCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCGCCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTg             >  1:1395144/1‑149 (MQ=255)
                                                                                                                                      cAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGcc           <  1:1214508/149‑1 (MQ=255)
                                                                                                                                       aGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCa          <  1:489688/149‑1 (MQ=255)
                                                                                                                                        gAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAg         <  2:1056955/149‑1 (MQ=255)
                                                                                                                                        gAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAg         >  1:891231/1‑149 (MQ=255)
                                                                                                                                        gAAATGCCGCGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAg         >  2:1288221/1‑149 (MQ=255)
                                                                                                                                         aaaTGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCCGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAGt        >  1:704017/1‑149 (MQ=255)
                                                                                                                                           aTGCCGTGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACAACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGGGCCAGTaa      >  2:1146222/1‑149 (MQ=255)
                                                                                                                                            tGCCGGGAATAATTCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGGTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAGTAAc     <  2:940995/149‑1 (MQ=255)
                                                                                                                                            tGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCa          <  2:920147/144‑1 (MQ=255)
                                                                                                                                              ccGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTGCCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAGTAACGa   <  2:209297/149‑1 (MQ=255)
                                                                                                                                              ccGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAGTAACGa   >  2:708356/1‑149 (MQ=255)
                                                                                                                                               cGGGATTTATTCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAGTAACGAg  >  2:1180518/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
CGTAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACGTTATTTCAGCAATATATTGAATTTACGTGCTTTTGTTGGCCGGACAAAGCGTTTACGCCGCATCCGGCATGAACAAAGCACACGTTGTTAACAATCAGAAATGCCGGGAATAAATCCCGGCATTTTCATAATCAGAAGTTGTAACCTACTACCAGGTAACCACCCCAGCCGGTAGAGCGAACGTTGAAGTTGCCGTTGCCGAAGTTCAGTTCTGCATCGTCGTTCCACTGACCACCGTCGTGCCAGTAACGAG  >  NC_000913/430959‑431250

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