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

ACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGC  >  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)
         ggTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCCCCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑At                                                                                                                                           >  2:695625/1‑150 (MQ=255)
         ggTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑At                                                                                                                                           <  2:1063834/150‑1 (MQ=255)
         ggTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCACCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGATGCTG‑At                                                                                                                                           >  1:322943/1‑150 (MQ=255)
          gTTTATTCAGCCAGACCAGCACCTACGGTCACAATGGTCGGGTGATAATCGCCCGTGGAAAGTGAAACGGTATGTCCTCTGAAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTTACGACATACAGCCTCCCGTGAGGCTG‑ATg                                                                                                                                          <  1:889378/150‑1 (MQ=255)
             aaTTCAGCCAGACCAGCACCGACGTTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCa                                                                                                                                       <  2:1195779/150‑1 (MQ=255)
                tCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAAc                                                                                                                                    <  2:1123414/150‑1 (MQ=255)
                 cAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAAcc                                                                                                                                   <  2:1016031/150‑1 (MQ=255)
                 cAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAAcc                                                                                                                                   <  1:1560042/150‑1 (MQ=255)
                  aGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAAccc                                                                                                                                  <  1:373260/150‑1 (MQ=255)
                     cAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCaaa                                                                                                                               >  1:966565/1‑150 (MQ=255)
                       gACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCaaaaa                                                                                                                             <  2:256934/150‑1 (MQ=255)
                          cAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCCCCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑CTGTCAAACCCAAAAACTg                                                                                                                          >  2:100484/1‑150 (MQ=255)
                           aGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGc                                                                                                                         <  1:1193664/150‑1 (MQ=255)
                             cACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATTTTGGTGATTTGCATACAAGCCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCAc                                                                                                                       >  2:1234348/1‑150 (MQ=255)
                              aCCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCAcc                                                                                                                      <  2:585195/150‑1 (MQ=255)
                                 gACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCggc                                                                                                                   >  1:1511262/1‑149 (MQ=255)
                                  aCGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGAc                                                                                                                  <  1:158219/150‑1 (MQ=255)
                                  aCGGTCACAATGGTCGCGTGATAATCCCCCGCGGAAAGTGAAACTGTCTGTCCACTGCAATGCAAATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGAc                                                                                                                  <  2:650904/150‑1 (MQ=255)
                                    ggTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACgg                                                                                                                >  2:1251906/1‑150 (MQ=255)
                                    ggTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACgg                                                                                                                >  1:48359/1‑150 (MQ=255)
                                       cacaATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTcc                                                                                                             >  2:1034472/1‑150 (MQ=255)
                                        acaATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACACAAAAACTGCACCGGACGGTccc                                                                                                            >  2:1671741/1‑150 (MQ=255)
                                         caATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCAACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTccca                                                                                                           >  1:1726286/1‑149 (MQ=255)
                                          aaTGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGACTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCt                                                                                                          >  1:37247/1‑150 (MQ=255)
                                             ggTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCCCCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAACACTGCACCGGACGGTCCCCTCGc                                                                                                       >  1:1741712/1‑150 (MQ=255)
                                             ggTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCATCGc                                                                                                       <  1:982715/150‑1 (MQ=255)
                                                cgcgTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGcccc                                                                                                    >  2:275964/1‑150 (MQ=255)
                                                     gATAATCCCCCGCGGCAAGTGAAACGGGCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTgg                                                                                               <  1:1812617/150‑1 (MQ=255)
                                                       tAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTgggg                                                                                             >  2:376543/1‑150 (MQ=255)
                                                       tAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGAGAATTGCATACAAGGCACCTTAATTAACGACATACCGCCTCCCGTGAGGCTG‑CTGTCAAACCCAAAAACTCCACCGGACGGTCCCCTCGCCCCTTTgggg                                                                                             >  2:944664/1‑150 (MQ=255)
                                                           cccccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAGTTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGagag                                                                                         <  2:1754421/150‑1 (MQ=255)
                                                                  ggCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAgg                                                                                  <  2:1572393/150‑1 (MQ=255)
                                                                  ggCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAgg                                                                                  <  2:1407088/150‑1 (MQ=255)
                                                                  ggCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAgg                                                                                  <  1:1470676/150‑1 (MQ=255)
                                                                   gCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAggg                                                                                 <  2:419712/150‑1 (MQ=255)
                                                                     aaGTGAAACGGGCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGgag                                                                               <  2:1793486/150‑1 (MQ=255)
                                                                       gTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATTCAAGGCACCTTAATTAACGACATACAGCCTCCCGGGCGGCTG‑ATGTCAAACCCATAAACTGCTCCGGACGGTCCCCTCGCCCCTTTTGGGTGTGGGTTAGGGagag                                                                             >  1:1722802/1‑150 (MQ=255)
                                                                         taaaCGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGCGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAggg                                                                           <  1:1158980/149‑1 (MQ=255)
                                                                            aCGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGaa                                                                        >  2:1458515/1‑150 (MQ=255)
                                                                              ggTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGaaaa                                                                      >  2:480706/1‑150 (MQ=255)
                                                                                 cTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAAccc                                                                   >  1:1337243/1‑150 (MQ=255)
                                                                                     ccACTGCAATGCATATTGGTGATTTGCATACAATGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑AAGTCAAACCCAAAAACTGCACAGGACGGTCCCCTCGCCCCGGTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCAc                                                               >  2:800651/1‑150 (MQ=255)
                                                                                     ccACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAACACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGTAGAGGGGAAACACCGCAc                                                               >  2:807556/1‑150 (MQ=255)
                                                                                     ccACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAACACCCGCAc                                                               >  2:677905/1‑150 (MQ=255)
                                                                                     ccACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCAc                                                               >  2:671485/1‑150 (MQ=255)
                                                                                     ccACTGCAATCCATATTGGTCATTTGATTACAAGGCACCTAAATTAACCTCATACAGCCTCCAGTGAGGCTT‑ATGTCAAACACAAAAACTGCACCGGACGGGCCCCTCGCCCATTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCAc                                                               <  1:115400/150‑1 (MQ=255)
                                                                                        cTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTAGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCCCCCt                                                            >  2:473646/1‑150 (MQ=255)
                                                                                        cTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCGG‑ATGTCAATCCCAAAAACTGCACCGGACGGTCCCCTCGCCCCATTGGGGAGAGGGTTAGGGAGCGGGGAAAACCCGCAcccg                                                            >  2:1295931/1‑149 (MQ=255)
                                                                                            aaTGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGATG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCACTTTGGGGAGAGGGTTAGGGAGAGTGGAAAACCCGCACCATAttt                                                        >  1:1242754/1‑150 (MQ=255)
                                                                                               gCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCa                                                     >  2:1029250/1‑150 (MQ=255)
                                                                                                cATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTTCGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAg                                                    <  2:655860/150‑1 (MQ=255)
                                                                                                 atatTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGa                                                   <  1:1163989/150‑1 (MQ=255)
                                                                                                  tatTGGTGATTTTCATACAAGGCACCTAAATAAACGACATACACCCTCCCGGGAGGCTG‑ATGTCAAACCCAAAATCTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCTCCCTATTTTCAGAt                                                  <  1:35548/150‑1 (MQ=38)
                                                                                                   atTTGTGTTTTGAAGACAAGGCACATTAATTAACGACATACAGCCTCCTTGGAGGCTG‑TTGTCAAACCAAAAAACTGCCCGGGACGGTCCCCTCGCCCCTTTGGGGAGAGTGTTAGGGATAGGGGAAAACCCGCACCCTATTTTCAGATc                                                 <  2:699830/150‑1 (MQ=255)
                                                                                                     ttgggATTTTCATACAAGCCAACTTAATTAAATACATACAGCCTCCCGTGAGGCTG‑ATGGCAAACCCAAAAACTCCACCGTACGGTCCCCTCGCCCATTTGGGGAGAGGGTTTGGGATAGGGGAAAACCCGCACCCTATTTTCAGATCAg                                               <  2:125179/146‑1 (MQ=255)
                                                                                                     tGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAg                                               >  2:51498/1‑150 (MQ=255)
                                                                                                     tGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAAACAGCACCCTATTTTCAGATCAg                                               >  2:1047244/1‑150 (MQ=255)
                                                                                                     tGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCATTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAg                                               >  2:1564341/1‑150 (MQ=255)
                                                                                                      ggTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGt                                              >  1:838687/1‑150 (MQ=255)
                                                                                                      ggTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGt                                              >  2:1209774/1‑150 (MQ=255)
                                                                                                           tttGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGACAACCCGCACCCTATTTTCAGATCAGTTAttt                                         >  1:1113009/1‑150 (MQ=255)
                                                                                                           tttGCATACAAGGCACCATAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCTGATCAGTTAttt                                         >  2:693758/1‑150 (MQ=255)
                                                                                                             tGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGAACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCg                                       >  2:637017/1‑150 (MQ=255)
                                                                                                               cATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCa                                     >  1:432146/1‑150 (MQ=255)
                                                                                                                aTACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCAt                                    <  1:1034472/150‑1 (MQ=255)
                                                                                                                aTACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCAt                                    >  1:523453/1‑150 (MQ=255)
                                                                                                                 tACAAGTCACATTAATTAACGACAAACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCAtt                                   <  2:1028636/150‑1 (MQ=255)
                                                                                                                             taattaaCGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGCAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATCTCcagcag                       >  2:492463/1‑150 (MQ=255)
                                                                                                                             taataacCGACATACCGCCCCCCGTGTGGCTG‑ATGTAAAACCCAAAAACTCCCCCGGACGGCCCCCTCCCCCTTTTGGGGAGAGGGTTAGTGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCcagcag                       <  1:636114/150‑1 (MQ=255)
                                                                                                                               attaaCGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCGTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGAc                     >  1:1417239/1‑150 (MQ=255)
                                                                                                                                  aaCGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCg                  >  2:19811/1‑150 (MQ=255)
                                                                                                                                  aaCGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCg                  >  2:694873/1‑150 (MQ=255)
                                                                                                                                  aaCGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCg                  >  2:158042/1‑150 (MQ=255)
                                                                                                                                   aCGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCACGCAGACACGc                 >  1:370916/1‑150 (MQ=255)
                                                                                                                                      gacatcAGCCTCACGGGAGGCTGTATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTg               >  2:1418320/6‑150 (MQ=39)
                                                                                                                                      aCATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGATAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGc              <  1:275964/150‑1 (MQ=255)
                                                                                                                                        catcAGCCTCACGGGAGGCTGTATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTGTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGCTATTTCGCATTAATATCCAGCAGACCCGCTGCa             >  2:1729906/4‑150 (MQ=255)
                                                                                                                                          aCAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAAcc          >  2:347690/1‑150 (MQ=255)
                                                                                                                                            aGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGc        <  2:48359/150‑1 (MQ=255)
                                                                                                                                             gCCTCCCGTGTGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCt       <  2:1232372/150‑1 (MQ=255)
                                                                                                                                             gCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACACGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCt       >  1:239321/1‑150 (MQ=255)
                                                                                                                                             gCCTCCCGTGAGGCTG‑ATGTCAAACACAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGcg       >  2:689708/1‑149 (MQ=255)
                                                                                                                                              ccTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTg      >  1:1058518/1‑150 (MQ=255)
                                                                                                                                                 cccGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAg   <  2:57375/150‑1 (MQ=255)
                                                                                                                                                  ccGTGTGGCTG‑TTGGCAAACCCAAAAACTCCCCCGGAGGGCCCCCACGCGCCTTGGGGGCGAGGGTTAGGGAGCGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGc  <  1:170667/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
ACCCTGAAAGGTTAATTCAGCCAGACCAGCACCGACGGTCACAATGGTCGCGTGATAATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGC  >  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: ‑