Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 503,597 G→A 8.4% D209N (GAC→AAC)  PP_0414 → polyamine ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947503,5970GA 8.4% 259.7 / 11.5 95D209N (GAC→AAC) PP_0414polyamine ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (51/35);  new base A (4/4);  total (56/39)
Fisher's exact test for biased strand distribution p-value = 7.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.87e-01

GACGCTGCAGGGGTTCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGA  >  NC_002947/503452‑503744
                                                                                                                                                 |                                                                                                                                                   
gACGCTGCAGGGGTTCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGa                                                                                                                                                 >  1:1899631/1‑150 (MQ=255)
    cTGCACGTGTTCAACTACAATCTTGTTCGTGCTGCGCCCAGGCTGGGCGCTTCATCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAggtg                                                                                                                                             <  2:728396/150‑1 (MQ=255)
       cAGGGGTTCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAggtggtg                                                                                                                                          >  2:237460/1‑150 (MQ=255)
         ggggTTCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAggtggtggt                                                                                                                                        >  2:379496/1‑150 (MQ=255)
         ggggTTCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAggtggtaaa                                                                                                                                        <  1:1498791/150‑4 (MQ=255)
         ggggTTCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCACAGTGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTCAAACGAggtggcaaa                                                                                                                                        <  1:747165/150‑5 (MQ=255)
          gggTTCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTCAAACGAggtggcgaag                                                                                                                                       <  1:2284890/150‑6 (MQ=255)
           ggTTCAACTACAATCTGGTGGGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGa                                                                                                                                      <  1:973051/150‑1 (MQ=255)
            gTTCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGAc                                                                                                                                     >  2:1551276/1‑150 (MQ=255)
             ttCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGGGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACg                                                                                                                                    >  1:472330/1‑150 (MQ=255)
                aaCTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTg                                                                                                                                 >  2:1042503/1‑150 (MQ=255)
                  cTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTCAAACGAggtggc                                                                                                                                           >  1:1027905/1‑137 (MQ=255)
                  cTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTCAAACGAggtggc                                                                                                                                           <  2:1027905/138‑2 (MQ=255)
                   tACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTc                                                                                                                              <  2:318046/150‑1 (MQ=255)
                             tgcgtgcgGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCgccgggc                                                                                                                    >  1:2307916/1‑150 (MQ=255)
                              gcgtgcgGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCgccgggcc                                                                                                                   <  2:651986/150‑1 (MQ=255)
                                    ggcggcCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAggt                                                                                                                                              >  2:1962037/1‑117 (MQ=255)
                                    ggcggcCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAggt                                                                                                                                              <  1:1962037/117‑1 (MQ=255)
                                        gcCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTACCACCTCGTCAAACGAggtggcaaag                                                                                                                                       >  2:611426/1‑115 (MQ=255)
                                        gcCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTACCACCTCGTCAAACGAggtggcaaag                                                                                                                                       <  1:611426/120‑6 (MQ=255)
                                             ccTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACt                                                                                                    >  2:2281028/1‑150 (MQ=255)
                                               tGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGc                                                                                                  >  2:2089593/1‑150 (MQ=255)
                                                 ggCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCcg                                                                                                >  1:1691729/1‑150 (MQ=255)
                                                  gcgcTTCACCTCTGCTGGCTTTCCGCCGTGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCcgc                                                                                               >  2:1256202/1‑150 (MQ=255)
                                                   cgcTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCcgcg                                                                                              <  1:772837/150‑1 (MQ=255)
                                                      ttCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCa                                                                                           <  2:1929137/150‑1 (MQ=255)
                                                      ttCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCa                                                                                           >  1:160731/1‑150 (MQ=255)
                                                          ccTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATg                                                                                       >  2:640511/1‑150 (MQ=255)
                                                           ctctGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCAACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGt                                                                                      >  1:2375489/1‑150 (MQ=255)
                                                            tctGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGtt                                                                                     >  2:1709488/1‑150 (MQ=255)
                                                              tgctgGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGGGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGGGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCa                                                                                   >  2:251272/1‑150 (MQ=255)
                                                               gctgGCTTTCCGCCGGGTGACCCTGCCGCTGATTGTGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAg                                                                                  <  1:1709532/150‑1 (MQ=255)
                                                                ctgGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAg                                                                                                                                                >  2:1646917/1‑87 (MQ=255)
                                                                ctgGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAg                                                                                                                                                <  1:1646917/87‑1 (MQ=255)
                                                                 tgGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCg                                                                                >  1:1783667/1‑150 (MQ=255)
                                                                 tgGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCAGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCg                                                                                >  2:806379/1‑150 (MQ=255)
                                                                         cgccgGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTcgc                                                                        <  1:2089593/150‑1 (MQ=255)
                                                                          gccgGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTTACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTcgcg                                                                       >  1:166373/1‑150 (MQ=255)
                                                                           ccgGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCga                                                                      <  2:458190/150‑1 (MQ=255)
                                                                            cgGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCgag                                                                     >  1:1374474/1‑150 (MQ=255)
                                                                                 gACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCt                                                                <  2:1899631/150‑1 (MQ=255)
                                                                                  aCCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCtc                                                               >  2:1494700/1‑149 (MQ=255)
                                                                                   cccTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGa                                                              <  2:472330/150‑1 (MQ=255)
                                                                                    ccTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAg                                                             >  2:689872/1‑150 (MQ=255)
                                                                                     cTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCt                                                                >  1:1652218/1‑146 (MQ=255)
                                                                                     cTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCt                                                                <  2:1652218/146‑1 (MQ=255)
                                                                                     cTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGc                                                            <  1:237460/150‑1 (MQ=255)
                                                                                     cTGCCGCTGATTGCGCCCGGTGTCATCTAGGGAGCCCTGTTTTCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGc                                                            >  2:521266/1‑150 (MQ=255)
                                                                                          gCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGAc                                                       >  1:865055/1‑150 (MQ=255)
                                                                                               ttGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTg                                                  <  1:1709488/150‑1 (MQ=255)
                                                                                                tGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGc                                                 <  2:166373/150‑1 (MQ=255)
                                                                                                 gcgcCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGcc                                                >  1:1710464/1‑150 (MQ=255)
                                                                                                     ccGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCgcgg                                            >  2:1748044/1‑150 (MQ=255)
                                                                                                     ccGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCgcgg                                            >  1:793161/1‑150 (MQ=255)
                                                                                                         tgtCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCgcggcggc                                        <  2:1691729/150‑1 (MQ=255)
                                                                                                             aTCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGAcc                                    >  2:641071/1‑150 (MQ=255)
                                                                                                              tctcGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCgccgg                                                                                                                      >  1:1030046/1‑67 (MQ=255)
                                                                                                              tctcGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCgccgg                                                                                                                      <  2:1030046/67‑1 (MQ=255)
                                                                                                                  gggAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTNGCt                               <  2:1512816/150‑1 (MQ=255)
                                                                                                                      gCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGAtt                           <  2:147955/150‑1 (MQ=255)
                                                                                                                       cccTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCgg                                                                               >  1:2194536/1‑97 (MQ=255)
                                                                                                                       cccTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCgg                                                                               <  2:2194536/97‑1 (MQ=255)
                                                                                                                       cccTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTg                          >  1:2385710/1‑150 (MQ=255)
                                                                                                                        ccTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACTTTGCTGATTGc                         <  2:1570913/150‑1 (MQ=255)
                                                                                                                         cTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCg                        >  1:381032/1‑150 (MQ=255)
                                                                                                                         cTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCAGCGGCGACCTTGCTGATTGCg                        >  2:1984995/1‑150 (MQ=255)
                                                                                                                          tGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGt                       >  1:1056911/1‑150 (MQ=255)
                                                                                                                          tGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGt                       >  1:2019855/1‑150 (MQ=255)
                                                                                                                            tttgcctttgccACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTtc                     >  1:939132/1‑150 (MQ=255)
                                                                                                                               gcctttgccACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCg                  <  2:865055/150‑1 (MQ=255)
                                                                                                                                cctttgccACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCgg                 >  1:2206581/1‑150 (MQ=255)
                                                                                                                                cctttgccACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCgg                 <  1:521266/150‑1 (MQ=255)
                                                                                                                                 ctttgccACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCAtt                                                   <  1:1763491/115‑1 (MQ=255)
                                                                                                                                 ctttgccACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCAtt                                                   >  2:1763491/1‑115 (MQ=255)
                                                                                                                                  tttgccACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAgg                                                                                                           <  1:1505922/58‑1 (MQ=255)
                                                                                                                                  tttgccACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAgg                                                                                                           >  2:1505922/1‑58 (MQ=255)
                                                                                                                                  tttgccACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGgtg               >  1:2309215/1‑150 (MQ=255)
                                                                                                                                  tttgccACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGgtg               <  2:1710464/150‑1 (MQ=255)
                                                                                                                                   ttgccACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACCCTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGgtgt              <  2:2307916/150‑1 (MQ=255)
                                                                                                                                   ccaccACCTCGTCAAACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGgtgt              >  2:1420874/4‑150 (MQ=255)
                                                                                                                                    caccACCTCGTCAAACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAg                                                             >  2:1270459/3‑102 (MQ=255)
                                                                                                                                     gccACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGtgttg            <  2:2019855/150‑1 (MQ=255)
                                                                                                                                     gccACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGtgttg            >  2:1438878/1‑150 (MQ=255)
                                                                                                                                      ccACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGtgttgt           >  1:1794319/1‑150 (MQ=255)
                                                                                                                                       cACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTgttgtt          >  2:666297/1‑150 (MQ=255)
                                                                                                                                        aCCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCCGGACATTCTTGCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGttgttg         >  2:1797174/1‑150 (MQ=255)
                                                                                                                                        aCCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCAtt                                                   <  1:114465/108‑1 (MQ=255)
                                                                                                                                        aCCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCAtt                                                   >  2:114465/1‑108 (MQ=255)
                                                                                                                                        aCCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGttgttg         >  2:1868753/1‑150 (MQ=255)
                                                                                                                                         ccTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTtgc        >  2:1886205/1‑150 (MQ=255)
                                                                                                                                           tCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATg                                                                                       <  1:710156/69‑1 (MQ=255)
                                                                                                                                           tCGTTTGACGAGGTGGTAGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATg                                                                                       >  2:710156/1‑69 (MQ=255)
                                                                                                                                             gTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTtgc        >  1:1660058/1‑146 (MQ=255)
                                                                                                                                             gTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTCCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTtgc        <  2:1660058/146‑1 (MQ=255)
                                                                                                                                               ttGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTCgctgctta  <  2:1056911/150‑3 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GACGCTGCAGGGGTTCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGCTGA  >  NC_002947/503452‑503744

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