Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,144,958 G→A 7.2% S108N (AGC→AAC)  PP_2760 → ribose ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,144,9580GA 7.2% 276.9 / 6.4 96S108N (AGC→AAC) PP_2760ribose ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (49/40);  new base A (3/4);  total (52/44)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

TTGCCTTCGGCCTGACCACGGTGATCATCGGCGGCGGCTCGAACGTGGTGTCCGGTGGGCTTGATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCTGGTGGCCGGGCTGGAACTGGTGCTGACCGAAAACACCGTGC  >  NC_002947/3144814‑3145107
                                                                                                                                                   |                                                                                                                                                     
ttgccttCGGCCTGACCACGGTGATCATCGGCGGCGGCTCGAACGTGGTGTCCGGTGGGCTTGATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATcg                                                                                                                                                  >  1:325130/1‑150 (MQ=255)
 tgccttCGGCCTGACCACGGTGATCATCGGCGGCGGCTCGAACGTGGTGTCCGGTGGGCTTGATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGCTGACAACATcgc                                                                                                                                                 <  1:1162644/150‑1 (MQ=255)
  gccttCGGCCTGCCCACTGGGATCATCGGCGGCGGCTCGAACGTGGTGTCCGTTGGGCTTGATCTGTCGCTGGCTGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGCTGACAACATcgaa                                                                                                                                                <  2:382453/150‑3 (MQ=255)
    cttCGGCCTGACCACGGTGATCATCGGCGGCGGCTCGAACGTGGTGTCCGGTGGGCTTGATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTAGGCGATGTTGGCAGCATcgctgg                                                                                                                                              >  1:272814/1‑147 (MQ=255)
     ttCGGCCTGACCACGGTGATCATCGGCGGCGGCTCGAACGTGGTGTCCGGTGGGCTTGATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCcgcg                                                                                                                                             >  2:273134/1‑150 (MQ=255)
     ttCGGCCTGACCACGGTGATCATCGGCGGCGGCTCGAACGTGGTGTCCGGTGGGCTTGATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGCTGCAAACATcgccg                                                                                                                                               >  2:839838/1‑148 (MQ=255)
     ttCGGCCTGACCACGGTGATCATCGGCGGCGGCTCGAACGTGGTGTCCGGTGGGCTTGATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGCTGACAACATcgccg                                                                                                                                               <  1:839838/148‑1 (MQ=255)
     ttCGGCCTGACCACGGTGATCATCGGCGGCGGCTCGAACGTGGTGTCCGGGGGGCTTGATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCcgcg                                                                                                                                             >  2:807548/1‑150 (MQ=255)
         gCCTGACCACGGTGATCATCGGCGGCGGCTCGAACGTGGTGTCCGGTGGGCTTGATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGCTGACAACATcgc                                                                                                                                                 >  2:507137/1‑142 (MQ=255)
         gCCTGACCACGGTGATCATCGGCGGCGGCTCGAACGTGGTGTCCGGTGGGCTTGATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGCTGACAACATcgc                                                                                                                                                 <  1:507137/142‑1 (MQ=255)
           cTGACCACGGTGATCATCGGCGGCGGCTCGAACGGGGTGTCCGGTGGGCTTGATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCtt                                                                                                                                       >  1:722490/1‑150 (MQ=255)
             gACCACGGTGATCATCGGCGGCGGCTCGAACGTGGTGTCCGGTGGGCTTGATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTgg                                                                                                                                     <  1:1689099/150‑1 (MQ=255)
                cACGGTGATCATCGGCGGCGGCTCGAACGGGGTGTCCGGTGGGCTTGATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCAc                                                                                                                                  >  2:1553114/1‑150 (MQ=255)
                           tcGGCGGCGGCTCGAACGTGGTGTCCGGTGGGCTTGATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTAGGCGGTGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCa                                                                                                                       >  1:407312/1‑150 (MQ=255)
                              gcggcggcTCGAACGTGGTGTCCGGTGGGCTTGATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCg                                                                                                                    >  2:341730/1‑150 (MQ=255)
                                     cTCGAACGTGGTGTCCGGTGGGCTTGATCTGTCGCTGGCGGCCAACCTCGGGTTGCGGGCGGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCt                                                                                                             <  1:1454268/150‑1 (MQ=255)
                                      tCGAACGTGGTGTCCGGTGGGCTTGATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTAGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTg                                                                                                            >  1:384598/1‑150 (MQ=255)
                                        gAACGTGGTGTCCGGTGGGCTTGATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGaa                                                                                                          <  2:830634/150‑1 (MQ=255)
                                        gAACGTGGTGTCCGGTGGGCTTGATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGaa                                                                                                          >  2:1305016/1‑150 (MQ=255)
                                              ggtgTCCGGTGGGCTTGATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCGCGCTGAGTGGGCt                                                                                                    >  2:1542432/1‑150 (MQ=255)
                                               gtgtCCGGTGGGCTTGATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCtg                                                                                                   <  2:1784169/150‑1 (MQ=255)
                                                 gtCCGGTGGGCTTGATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGCCCGTTGG‑‑AGGCTGAACAGCGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCtggt                                                                                                 <  2:731357/150‑1 (MQ=255)
                                                  tCCGGTGGGCTTGATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCtggtg                                                                                                >  1:658334/1‑150 (MQ=255)
                                                         ggCTTGATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGGCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGGGGTCg                                                                                         >  1:1853501/1‑150 (MQ=255)
                                                          gCTTGATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCgg                                                                                        >  1:251504/1‑150 (MQ=255)
                                                           cTTGATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGGGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTTGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGt                                                                                       >  1:411732/1‑150 (MQ=255)
                                                           cTTGATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATc                                                                                                                     <  2:312498/120‑1 (MQ=255)
                                                           cTTGATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCCCGACGCTTGGCACTGGCCTGGCCATc                                                                                                                     >  1:312498/1‑120 (MQ=255)
                                                            ttGATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCtt                                                                                                                                       <  2:546985/101‑1 (MQ=255)
                                                            ttGATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCtt                                                                                                                                       >  1:546985/1‑101 (MQ=255)
                                                            ttGATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGtt                                                                                      >  1:486422/1‑150 (MQ=255)
                                                              gATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTc                                                                                    >  2:1412254/1‑150 (MQ=255)
                                                              gATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTc                                                                                    >  1:268970/1‑150 (MQ=255)
                                                              gATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTc                                                                                    <  2:909941/150‑1 (MQ=255)
                                                                gcTGTCGCGGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCg                                                                                  >  2:505642/2‑150 (MQ=255)
                                                                 cTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGt                                                                                 <  1:341730/150‑1 (MQ=255)
                                                                  tGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGtt                                                                                >  1:933915/1‑150 (MQ=255)
                                                                  tGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGtt                                                                                >  2:461673/1‑150 (MQ=255)
                                                                  tGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGtt                                                                                >  2:1554875/1‑150 (MQ=255)
                                                                  tGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCACGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGtt                                                                                >  1:1544089/1‑150 (MQ=255)
                                                                   gTCGCTGGCGGCCAACCTCGTGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTTGTGCGCAATTTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGttt                                                                               <  1:1553114/150‑1 (MQ=255)
                                                                          gcggcCAACCTCGGGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGccgccg                                                                        <  2:1231792/150‑1 (MQ=255)
                                                                            ggcCAACCTCGGGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCgct                                                                      >  1:1749447/1‑150 (MQ=255)
                                                                            ggcCAACCTCGGGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCgct                                                                      >  1:1229335/1‑150 (MQ=255)
                                                                               cAACCTCGGGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCgctgct                                                                   >  1:92214/1‑150 (MQ=255)
                                                                               cAACCTCGGGTTGTCCGCTGCG‑‑CCATGCCCGAGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCgctgct                                                                   <  1:941198/150‑1 (MQ=255)
                                                                                aaCCTCGGGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTAAGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGctgctg                                                                  <  1:1367084/150‑1 (MQ=255)
                                                                                  ccTCGGGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGc                                                                >  1:1827325/1‑150 (MQ=255)
                                                                                    tCGGGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCtggtgg                                                                                               <  1:1689832/117‑1 (MQ=255)
                                                                                    tCGGGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCtggtgg                                                                                               >  2:1689832/1‑117 (MQ=255)
                                                                                    tCGGGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGttt                                                                               >  2:794364/1‑133 (MQ=255)
                                                                                    tCGGGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGttt                                                                               >  2:794328/1‑133 (MQ=255)
                                                                                    tCGGGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGttt                                                                               <  1:794364/133‑1 (MQ=255)
                                                                                    tCGGGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGttt                                                                               <  1:794328/133‑1 (MQ=255)
                                                                                    tCGGGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGa                                                              >  1:984489/1‑150 (MQ=255)
                                                                                      ggattGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGCCGGTTTCCGTTTGCCGCCGCTGCTGGCGAcc                                                            >  2:371363/4‑150 (MQ=255)
                                                                                       ggTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGAccc                                                           >  2:1107230/1‑150 (MQ=255)
                                                                                        gTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGttt                                                                                     >  1:1011780/1‑123 (MQ=255)
                                                                                        cttGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGttt                                                                                     <  2:1011780/122‑1 (MQ=255)
                                                                                         ttGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTg                                                         >  1:574323/1‑150 (MQ=255)
                                                                                          tGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTgg                                                        >  2:818794/1‑150 (MQ=255)
                                                                                          tGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTgg                                                        <  1:1412254/150‑1 (MQ=255)
                                                                                          tGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTgg                                                        <  2:1890957/150‑1 (MQ=255)
                                                                                              cGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGAcccgaggttg                                                    <  2:92214/150‑8 (MQ=255)
                                                                                              cGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCAc                                                    <  2:325130/150‑1 (MQ=255)
                                                                                                cTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGGGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGa                                                  >  1:150784/1‑150 (MQ=255)
                                                                                                 tGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGAt                                                 <  1:1258760/150‑1 (MQ=255)
                                                                                                   cGGTCTATGCC‑‑AGCCTGAACAACGCCGTTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGa                                               <  2:891227/150‑1 (MQ=255)
                                                                                                   cGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGa                                               <  2:1747754/150‑1 (MQ=255)
                                                                                                    ggTCTATGCC‑‑AGCCTGAACAACGCCGGTATCGCTGATGTTTTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGaa                                              <  2:150784/150‑1 (MQ=255)
                                                                                                      tCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCCTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAAcc                                            >  2:889652/1‑150 (MQ=255)
                                                                                                       cTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCTCCTGGTGGCGACCCTGGCCACGATGAACCt                                           <  1:51866/150‑1 (MQ=255)
                                                                                                        tATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGa                                                  >  1:1171996/1‑142 (MQ=255)
                                                                                                        tATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGa                                                  <  2:1171996/142‑1 (MQ=255)
                                                                                                        tATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCtg                                          <  1:1826251/150‑1 (MQ=255)
                                                                                                         aTGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCtgg                                         <  1:807548/150‑1 (MQ=255)
                                                                                                           gCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCtggtg                                       <  2:486422/150‑1 (MQ=255)
                                                                                                           gCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTAAGCGTAGGTGCGACGTTTGTCACTGTCCGGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCtggtg                                       <  2:1366158/150‑1 (MQ=255)
                                                                                                            cc‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCtggtgg                                      >  1:596455/1‑150 (MQ=255)
                                                                                                                   cTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCTGGTGGCCggg                                 <  1:824361/150‑1 (MQ=255)
                                                                                                                    tGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCTGGTGGCCGGGc                                >  2:1738214/1‑150 (MQ=255)
                                                                                                                         aacGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGTCTACCGGCGCATTCCTTTTTCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCTGGTGGCCGGGCTGGaa                           >  1:974866/1‑150 (MQ=255)
                                                                                                                          acGCCGGTTTCGGCTATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCTGGTGGCCGGGCTGGAAc                          <  1:1271277/150‑1 (MQ=255)
                                                                                                                            gCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTCCCTATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCTGGTGGCCGGGCTGGAACTg                        <  1:1366481/150‑1 (MQ=255)
                                                                                                                                 tttCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCTGGTGGCCGGGCTGGAACTGGTGCt                   >  2:543150/1‑150 (MQ=255)
                                                                                                                                 tttCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGACCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCTGGTGGCCGGGCTGGAACTg                        <  2:1067673/145‑1 (MQ=255)
                                                                                                                                 tttCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGACCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCTGGTGGCCGGGCTGGAACTg                        >  1:1067673/1‑145 (MQ=255)
                                                                                                                                    cGGCGATGCTGACAACATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCTGGTGGCCGGGCTGGAACTGGTGCTGAc                >  2:827517/1‑150 (MQ=255)
                                                                                                                                       gCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCTGGTGGCCGGGCTGGAACTGGTGCTGACCg              >  2:137554/1‑150 (MQ=255)
                                                                                                                                         gATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCTGGTGGCCGGGCTGGAACTGGTGCTGACCGaa            <  2:251504/150‑1 (MQ=255)
                                                                                                                                         gATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCTGGTGGCCGGGCTGGAACTGGTGCTGACCGaa            <  1:461673/150‑1 (MQ=255)
                                                                                                                                             ttgtCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCTGGTGGCCGGGCTGGAACTGGTGCTGACCGAAAaca        >  2:1546255/1‑150 (MQ=255)
                                                                                                                                             ttgtCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCTGGTGGCCGGGCTGGAACTGGTGCTGACCGAAAaca        <  1:1738214/150‑1 (MQ=255)
                                                                                                                                                  aGCATCGCCGCGACGCTTGGCAGTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGGTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCTGGTGGCCGGGCTGGATCTGGTGCTGACCGAAAACACCGtg   <  2:272814/150‑1 (MQ=255)
                                                                                                                                                  aGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCTGGTGGCCGGGCTGGAACTGGTGCTGACCGAAAACACCGtg   >  1:937013/1‑150 (MQ=255)
                                                                                                                                                  aGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCTGGTGGCCGGGCTGGAACTGGTGCTGACCGAAAACACCGtg   <  1:1546255/150‑1 (MQ=255)
                                                                                                                                                   gCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCTGGTGGCCGGGCTGGAACTGGTGCTGACCGAAAACACCGtgc  >  2:840598/1‑150 (MQ=255)
                                                                                                                                                   gCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCCATGCTGGCGACCCTGGCCACGATGAACCGGGTGGCCGGGCTGGAACTGGTGCTGACCGAAAACACCGtgc  >  2:686859/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TTGCCTTCGGCCTGACCACGGTGATCATCGGCGGCGGCTCGAACGTGGTGTCCGGTGGGCTTGATCTGTCGCTGGCGGCCAACCTCGGGTTGTCCGCTGCGGTCTATGCC‑‑AGCCTGAACAACGCCGGTTTCGGCGATGTTGTCAGCATCGCCGCGACGCTTGGCACTGGCCTGGCCATCGGCACGCTGAATGCGCTGGTGGTGGTCGGTTTCCGTTTGCCGCCGCTGCTGGCGACCCTGGCCACGATGAACCTGGTGGCCGGGCTGGAACTGGTGCTGACCGAAAACACCGTGC  >  NC_002947/3144814‑3145107

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