Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 503,594 G→A 7.3% 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,5940GA 7.3% 196.7 / 4.2 68D209N (GAC→AAC) PP_0414polyamine ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (33/30);  new base A (2/3);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.61e-01

CGACGCTGCAGGGGTTCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGC  >  NC_002947/503448‑503738
                                                                                                                                                  |                                                                                                                                                
cGACGCTGCAGGGGTTCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGAc                                                                                                                                                <  1:169522/149‑1 (MQ=255)
        cAGGGGTTCAACTACAATCTGGTGCGTGCGTCGGCCAGCATGTGCGCTTCACCTCTGCTGTCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGATCCCTGTTTGCCTTTGACACCTCGTTTGACGAggtggt                                                                                                                                        >  8:156344/1‑149 (MQ=255)
        cAGGGGTTCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTACCACCTCGTCAAACGAggtggc                                                                                                                                        <  6:199978/149‑2 (MQ=255)
          ggggTTCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAggtggtgg                                                                                                                                      <  2:146696/149‑1 (MQ=255)
           gggTTCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAggtggtggt                                                                                                                                     <  3:95985/149‑1 (MQ=255)
            ggTTCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGTGAGCCCTGTTTGCTTTTACCACCTCGTCAAACGAggtggcaaag                                                                                                                                    <  1:24503/149‑6 (MQ=255)
            ggTTCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTACCACCTCGTCAAACGAggtggcaaag                                                                                                                                    <  2:212792/149‑6 (MQ=255)
              ttCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGAc                                                                                                                                  >  6:164634/1‑149 (MQ=255)
                  aCTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTg                                                                                                                              <  6:126212/149‑1 (MQ=255)
                   gttcAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCGTTCCGCCGGTAGACCCTGCCGCTGATTGCGCCGGGTGTCATCTCGGGAGCCGTGTTTGCCTTTGCCACGTCGTTTGACGAGGGGGTGGGGCCGCGGt                                                                                                                             >  5:187904/4‑149 (MQ=255)
                   cTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGa                                                                                                                                   >  5:203380/1‑143 (MQ=255)
                   cTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGa                                                                                                                                   <  6:203380/143‑1 (MQ=255)
                   cTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCATTACGCCGGGTGACCCTGCCGCTGATGGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGt                                                                                                                             >  5:187903/1‑149 (MQ=255)
                    tACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCAGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGtt                                                                                                                            >  5:216757/1‑149 (MQ=255)
                    tACAATCAGGTGCGTGCGGCGGCCAGCCTCGGCGCTTCCCCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATAGCGCCCGGAGTCATCTCGGGAGCCCTGTTTGCCTTAGCCACCTCGTTAGACGAGGTGGTGGTGACGCTGtt                                                                                                                            >  5:198749/1‑149 (MQ=255)
                         tCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCg                                                                                                                       >  3:164138/1‑149 (MQ=255)
                                 gtgcgGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCgccgggccg                                                                                                               >  6:24936/1‑149 (MQ=255)
                                   gcggcgGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGa                                                                                                             <  4:158958/149‑1 (MQ=255)
                                      gcggcCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCa                                                                                                          <  6:216757/149‑1 (MQ=255)
                                        ggcCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAgg                                                                                                        <  3:80337/149‑1 (MQ=255)
                                           cAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCaa                                                                                                     >  3:38070/1‑149 (MQ=255)
                                           cAGCCTGGGCGCTTCACCTATGCTGGCTTTACGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTACTCGCCGGGCCGGAGCAGGCaa                                                                                                     >  5:56871/1‑149 (MQ=255)
                                                       ttCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCaa                                                                                                     <  5:63397/137‑1 (MQ=255)
                                                       ttCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCaa                                                                                                     >  6:63397/1‑137 (MQ=255)
                                                         cACCTCTTCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGC                                                                                         <  1:245715‑M1/147‑1 (MQ=255)
                                                                    gCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGG                                                                                          <  3:41356‑M1/135‑1 (MQ=255)
                                                                    gCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGG                                                                                          >  4:41356‑M1/1‑135 (MQ=255)
                                                                    gCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGACTTTGCCACATCGTTTGACGAGGTGGTGGTGACGCTCTTCCTCGCCGGGACGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCgc                                                                            >  4:202096/1‑148 (MQ=255)
                                                                         ccgccgGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTc                                                                       >  3:265428/1‑149 (MQ=255)
                                                                          cgccgGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCG                                                                                             <  7:182746‑M1/126‑1 (MQ=255)
                                                                          cgccgGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCG                                                                                             >  8:182746‑M1/1‑126 (MQ=255)
                                                                           gccgGGTGACACTGCCGCTGATTGCGCCCGGTGTCATCTCTGGAGCCCTGTTTGCCTTTGCAACCTCGTTTGACGAGGGTGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCTGCAGATGTTCAGCGGTATTcgc                                                                     <  7:156344/149‑1 (MQ=255)
                                                                             cgGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCga                                                                   <  4:164138/149‑1 (MQ=255)
                                                                              gggTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTcgc                                                                     >  5:42100/1‑146 (MQ=255)
                                                                                 tGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAAc                                                               >  1:114746/1‑149 (MQ=255)
                                                                                    cccTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCt                                                             >  5:16995/1‑148 (MQ=255)
                                                                                     ccTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCt                                                             >  4:146701/1‑147 (MQ=255)
                                                                                      cgtccGCTGATTGAGCCCGGGGTGAGCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAg                                                          <  6:187904/146‑1 (MQ=255)
                                                                                      cTGCCGCTGCTTGCGCACGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAg                                                          <  6:187903/149‑1 (MQ=255)
                                                                                       tGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGc                                                         >  2:125407/1‑149 (MQ=255)
                                                                                       tGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGc                                                         >  8:243391/1‑149 (MQ=255)
                                                                                         ccGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGccc                                                       >  5:6785/1‑149 (MQ=255)
                                                                                         ccGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGc                                                         <  4:230121/147‑1 (MQ=255)
                                                                                            cTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGAc                                                    <  2:114746/149‑1 (MQ=255)
                                                                                                ttGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATt                                                <  6:2182‑M2/149‑2 (MQ=255)
                                                                                                      ccGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGCCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCCACCATtgccgag                                          <  8:70523‑M2/149‑8 (MQ=255)
                                                                                                              aTCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGAAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGAc                                  >  1:263848/1‑149 (MQ=255)
                                                                                                               tctcGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGAAGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGAcc                                 >  4:239969/1‑149 (MQ=255)
                                                                                                                   gggAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAatggtcgggctcaggttctcgcgaataccgctgaacatctgccgcggc                                        <  8:38863‑M1/138‑49 (MQ=255)
                                                                                                                   gggAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAatggtcgggctcaggttctcgcgaataccgctgaacatctgccgcggc                                        >  7:38863‑M1/1‑90 (MQ=255)
                                                                                                                      aGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGTCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTCGCTGa                          <  7:243391/149‑1 (MQ=255)
                                                                                                                         ccTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAatggtcgggctcaggttctcgcgaataccgctgaacatctgccgcggcagtgtt                                  <  1:235429‑M1/138‑55 (MQ=255)
                                                                                                                         ccTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAatggtcgggctcaggttctcgcgaataccgctgaacatctgccgcggcagtgtt                                  >  2:235429‑M1/1‑84 (MQ=255)
                                                                                                                          cTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGc                      <  6:42100/149‑1 (MQ=255)
                                                                                                                           tGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCg                     >  3:47628/1‑149 (MQ=255)
                                                                                                                               tgcctttgccACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTtct                 >  5:217887/1‑149 (MQ=255)
                                                                                                                               tgcctttgccACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTtc                  <  3:146701/148‑1 (MQ=255)
                                                                                                                                 cctttgccACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTtctc                <  4:265428/148‑1 (MQ=255)
                                                                                                                                  ctttgccACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTAGCTGATTGCGTTCTCgg              >  1:240608/1‑149 (MQ=255)
                                                                                                                                     tgccACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGgt             <  6:217887/147‑1 (MQ=255)
                                                                                                                                     caccACCTCGTCAAACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGgtgt           >  6:226289/3‑149 (MQ=255)
                                                                                                                                     caccACCTCGTCAAACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGgtgt           >  8:217071/3‑149 (MQ=255)
                                                                                                                                        cACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGa                                                           <  1:23478/98‑1 (MQ=255)
                                                                                                                                        cACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGa                                                           >  2:23478/1‑98 (MQ=255)
                                                                                                                                        cACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGtgttgt        >  2:29533/1‑149 (MQ=255)
                                                                                                                                              gTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGc                             >  7:224290/1‑122 (MQ=255)
                                                                                                                                              gTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGc                             <  8:224290/122‑1 (MQ=255)
                                                                                                                                              gTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTtgctgc  <  6:16995/149‑1 (MQ=255)
                                                                                                                                                  gACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATtgcc                                             >  2:249711‑M2/1‑98 (MQ=255)
                                                                                                                                                  gACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATtgcc                                             <  1:249711‑M2/102‑5 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
CGACGCTGCAGGGGTTCAACTACAATCTGGTGCGTGCGGCGGCCAGCCTGGGCGCTTCACCTCTGCTGGCTTTCCGCCGGGTGACCCTGCCGCTGATTGCGCCCGGTGTCATCTCGGGAGCCCTGTTTGCCTTTGCCACCTCGTTTGACGAGGTGGTGGTGACGCTGTTCCTCGCCGGGCCGGAGCAGGCAACACTGCCGCGGCAGATGTTCAGCGGTATTCGCGAGAACCTGAGCCCGACCATTGCCGCGGCGGCGACCTTGCTGATTGCGTTCTCGGTGTTGTTGCTGC  >  NC_002947/503448‑503738

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