Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,070,413 T→G 6.8% intergenic (‑12/+102) yihR ← / ← yihS putative sulphoquinovose mutarotase/sulphoquinovose isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,070,4130TG 6.8% 226.8 / 4.8 91intergenic (‑12/+102)yihR/yihSputative sulphoquinovose mutarotase/sulphoquinovose isomerase
Reads supporting (aligned to +/- strand):  ref base T (47/37);  new base G (3/4);  total (50/41)
Fisher's exact test for biased strand distribution p-value = 6.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01

ACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGC  >  NC_000913/4070264‑4070559
                                                                                                                                                     |                                                                                                                                                   
aCCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCCTACAAGGCAGCTTAATTAACGACATACAGCCTCCCGt                                                                                                                                                     >  1:259345/1‑150 (MQ=255)
 cccTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTg                                                                                                                                                    <  2:1283724/150‑1 (MQ=255)
   cTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAg                                                                                                                                                  <  1:991724/150‑1 (MQ=255)
    tGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAgg                                                                                                                                                 <  2:1323216/150‑1 (MQ=255)
     gAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGc                                                                                                                                                <  2:1593417/150‑1 (MQ=255)
        aGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCACGGGAGGCTgt                                                                                                                                             <  1:1199336/150‑2 (MQ=255)
        aGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCACGGGAGGCTgt                                                                                                                                             <  1:1744868/150‑2 (MQ=255)
         ggTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCCCCTTAATTAACGACATACAGCCTCCCGTGAGGCTGAt                                                                                                                                           >  2:695625/1‑150 (MQ=255)
         ggTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGAt                                                                                                                                           <  2:1063834/150‑1 (MQ=255)
         ggTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCACCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGATGCTGAt                                                                                                                                           >  1:322943/1‑150 (MQ=255)
          gTTTATTCAGCCAGACCAGCACCTACGGTCACAATGGTCGGGTGATAATCGCCCGTGGAAAGTGAAACGGTATGTCCTCTGAAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTTACGACATACAGCCTCCCGTGAGGCTGATg                                                                                                                                          <  1:889378/150‑1 (MQ=255)
             aaTTCAGCCAGACCAGCACCGACGTTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCa                                                                                                                                       <  2:1195779/150‑1 (MQ=255)
                tCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAAc                                                                                                                                    <  2:1123414/150‑1 (MQ=255)
                 cAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAAcc                                                                                                                                   <  2:1016031/150‑1 (MQ=255)
                 cAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAAcc                                                                                                                                   <  1:1560042/150‑1 (MQ=255)
                  aGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAAccc                                                                                                                                  <  1:373260/150‑1 (MQ=255)
                     cAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCaaa                                                                                                                               >  1:966565/1‑150 (MQ=255)
                       gACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCaaaaa                                                                                                                             <  2:256934/150‑1 (MQ=255)
                          cAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCCCCTTAATTAACGACATACAGCCTCCCGTGAGGCTGCTGTCAAACCCAAAAACTg                                                                                                                          >  2:100484/1‑150 (MQ=255)
                           aGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGc                                                                                                                         <  1:1193664/150‑1 (MQ=255)
                             cACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATTTTGGTGATTTGCATACAAGCCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCAc                                                                                                                       >  2:1234348/1‑150 (MQ=255)
                              aCCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCAcc                                                                                                                      <  2:585195/150‑1 (MQ=255)
                                 gACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCggc                                                                                                                   >  1:1511262/1‑149 (MQ=255)
                                  aCGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGAc                                                                                                                  <  1:158219/150‑1 (MQ=255)
                                  aCGGTCACAATGGTCGCGTGATAATCCCCCGCGGAAAGTGAAACTGTCTGTCCACTGCAATGCAAATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGAc                                                                                                                  <  2:650904/150‑1 (MQ=255)
                                    ggTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACgg                                                                                                                >  2:1251906/1‑150 (MQ=255)
                                    ggTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACgg                                                                                                                >  1:48359/1‑150 (MQ=255)
                                       cacaATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTcc                                                                                                             >  2:1034472/1‑150 (MQ=255)
                                        acaATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACACAAAAACTGCACCGGACGGTccc                                                                                                            >  2:1671741/1‑150 (MQ=255)
                                         caATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCAACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTccca                                                                                                           >  1:1726286/1‑149 (MQ=255)
                                          aaTGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGACTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCt                                                                                                          >  1:37247/1‑150 (MQ=255)
                                             ggTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCCCCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAACACTGCACCGGACGGTCCCCTCGc                                                                                                       >  1:1741712/1‑150 (MQ=255)
                                             ggTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCATCGc                                                                                                       <  1:982715/150‑1 (MQ=255)
                                                cgcgTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGcccc                                                                                                    >  2:275964/1‑150 (MQ=255)
                                                     gATAATCCCCCGCGGCAAGTGAAACGGGCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTgg                                                                                               <  1:1812617/150‑1 (MQ=255)
                                                       tAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTgggg                                                                                             >  2:376543/1‑150 (MQ=255)
                                                       tAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGAGAATTGCATACAAGGCACCTTAATTAACGACATACCGCCTCCCGTGAGGCTGCTGTCAAACCCAAAAACTCCACCGGACGGTCCCCTCGCCCCTTTgggg                                                                                             >  2:944664/1‑150 (MQ=255)
                                                           cccccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAGTTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGagag                                                                                         <  2:1754421/150‑1 (MQ=255)
                                                                  ggCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAgg                                                                                  <  2:1572393/150‑1 (MQ=255)
                                                                  ggCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAgg                                                                                  <  2:1407088/150‑1 (MQ=255)
                                                                  ggCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAgg                                                                                  <  1:1470676/150‑1 (MQ=255)
                                                                   gCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAggg                                                                                 <  2:419712/150‑1 (MQ=255)
                                                                     aaGTGAAACGGGCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGgag                                                                               <  2:1793486/150‑1 (MQ=255)
                                                                       gTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATTCAAGGCACCTTAATTAACGACATACAGCCTCCCGGGCGGCTGATGTCAAACCCATAAACTGCTCCGGACGGTCCCCTCGCCCCTTTTGGGTGTGGGTTAGGGagag                                                                             >  1:1722802/1‑150 (MQ=255)
                                                                         taaaCGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGCGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAggg                                                                           <  1:1158980/149‑1 (MQ=255)
                                                                            aCGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGaa                                                                        >  2:1458515/1‑150 (MQ=255)
                                                                              ggTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGaaaa                                                                      >  2:480706/1‑150 (MQ=255)
                                                                                 cTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAAccc                                                                   >  1:1337243/1‑150 (MQ=255)
                                                                                     ccACTGCAATGCATATTGGTGATTTGCATACAATGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGAAGTCAAACCCAAAAACTGCACAGGACGGTCCCCTCGCCCCGGTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCAc                                                               >  2:800651/1‑150 (MQ=255)
                                                                                     ccACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAACACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGTAGAGGGGAAACACCGCAc                                                               >  2:807556/1‑150 (MQ=255)
                                                                                     ccACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAACACCCGCAc                                                               >  2:677905/1‑150 (MQ=255)
                                                                                     ccACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCAc                                                               >  2:671485/1‑150 (MQ=255)
                                                                                     ccACTGCAATCCATATTGGTCATTTGATTACAAGGCACCTAAATTAACCTCATACAGCCTCCAGTGAGGCTTATGTCAAACACAAAAACTGCACCGGACGGGCCCCTCGCCCATTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCAc                                                               <  1:115400/150‑1 (MQ=255)
                                                                                        cTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTAGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCCCCCt                                                            >  2:473646/1‑150 (MQ=255)
                                                                                        cTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCGGATGTCAATCCCAAAAACTGCACCGGACGGTCCCCTCGCCCCATTGGGGAGAGGGTTAGGGAGCGGGGAAAACCCGCAcccg                                                            >  2:1295931/1‑149 (MQ=255)
                                                                                            aaTGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGATGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCACTTTGGGGAGAGGGTTAGGGAGAGTGGAAAACCCGCACCATAttt                                                        >  1:1242754/1‑150 (MQ=255)
                                                                                               gCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCa                                                     >  2:1029250/1‑150 (MQ=255)
                                                                                                cATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTTCGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAg                                                    <  2:655860/150‑1 (MQ=255)
                                                                                                 atatTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGa                                                   <  1:1163989/150‑1 (MQ=255)
                                                                                                  tatTGGTGATTTTCATACAAGGCACCTAAATAAACGACATACACCCTCCCGGGAGGCTGATGTCAAACCCAAAATCTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCTCCCTATTTTCAGAt                                                  <  1:35548/150‑1 (MQ=38)
                                                                                                   atTTGTGTTTTGAAGACAAGGCACATTAATTAACGACATACAGCCTCCTTGGAGGCTGTTGTCAAACCAAAAAACTGCCCGGGACGGTCCCCTCGCCCCTTTGGGGAGAGTGTTAGGGATAGGGGAAAACCCGCACCCTATTTTCAGATc                                                 <  2:699830/150‑1 (MQ=255)
                                                                                                     ttgggATTTTCATACAAGCCAACTTAATTAAATACATACAGCCTCCCGTGAGGCTGATGGCAAACCCAAAAACTCCACCGTACGGTCCCCTCGCCCATTTGGGGAGAGGGTTTGGGATAGGGGAAAACCCGCACCCTATTTTCAGATCAg                                               <  2:125179/146‑1 (MQ=255)
                                                                                                     tGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAg                                               >  2:51498/1‑150 (MQ=255)
                                                                                                     tGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAAACAGCACCCTATTTTCAGATCAg                                               >  2:1047244/1‑150 (MQ=255)
                                                                                                     tGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCATTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAg                                               >  2:1564341/1‑150 (MQ=255)
                                                                                                      ggTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGt                                              >  1:838687/1‑150 (MQ=255)
                                                                                                      ggTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGt                                              >  2:1209774/1‑150 (MQ=255)
                                                                                                           tttGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGACAACCCGCACCCTATTTTCAGATCAGTTAttt                                         >  1:1113009/1‑150 (MQ=255)
                                                                                                           tttGCATACAAGGCACCATAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCTGATCAGTTAttt                                         >  2:693758/1‑150 (MQ=255)
                                                                                                             tGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGAACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCg                                       >  2:637017/1‑150 (MQ=255)
                                                                                                               cATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCa                                     >  1:432146/1‑150 (MQ=255)
                                                                                                                aTACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCAt                                    <  1:1034472/150‑1 (MQ=255)
                                                                                                                aTACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCAt                                    >  1:523453/1‑150 (MQ=255)
                                                                                                                 tACAAGTCACATTAATTAACGACAAACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCAtt                                   <  2:1028636/150‑1 (MQ=255)
                                                                                                                             taattaaCGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGCAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATCTCcagcag                       >  2:492463/1‑150 (MQ=255)
                                                                                                                             taataacCGACATACCGCCCCCCGTGTGGCTGATGTAAAACCCAAAAACTCCCCCGGACGGCCCCCTCCCCCTTTTGGGGAGAGGGTTAGTGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCcagcag                       <  1:636114/150‑1 (MQ=255)
                                                                                                                               attaaCGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCGTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGAc                     >  1:1417239/1‑150 (MQ=255)
                                                                                                                                  aaCGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCg                  >  2:19811/1‑150 (MQ=255)
                                                                                                                                  aaCGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCg                  >  2:694873/1‑150 (MQ=255)
                                                                                                                                  aaCGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCg                  >  2:158042/1‑150 (MQ=255)
                                                                                                                                   aCGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCACGCAGACACGc                 >  1:370916/1‑150 (MQ=255)
                                                                                                                                      gacatcAGCCTCACGGGAGGCTGTATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTg               >  2:1418320/6‑150 (MQ=39)
                                                                                                                                      aCATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGATAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGc              <  1:275964/150‑1 (MQ=255)
                                                                                                                                        catcAGCCTCACGGGAGGCTGTATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTGTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGCTATTTCGCATTAATATCCAGCAGACCCGCTGCa             >  2:1729906/4‑150 (MQ=255)
                                                                                                                                          aCAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAAcc          >  2:347690/1‑150 (MQ=255)
                                                                                                                                            aGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGc        <  2:48359/150‑1 (MQ=255)
                                                                                                                                             gCCTCCCGTGTGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCt       <  2:1232372/150‑1 (MQ=255)
                                                                                                                                             gCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACACGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCt       >  1:239321/1‑150 (MQ=255)
                                                                                                                                             gCCTCCCGTGAGGCTGATGTCAAACACAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGcg       >  2:689708/1‑149 (MQ=255)
                                                                                                                                              ccTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTg      >  1:1058518/1‑150 (MQ=255)
                                                                                                                                                 cccGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAg   <  2:57375/150‑1 (MQ=255)
                                                                                                                                                  ccGTGTGGCTGTTGGCAAACCCAAAAACTCCCCCGGAGGGCCCCCACGCGCCTTGGGGGCGAGGGTTAGGGAGCGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGc  <  1:170667/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
ACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGC  >  NC_000913/4070264‑4070559

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: