Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,542,321 C→T 6.6% P134P (CCG→CCA galE ← UDP‑glucose 4‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,542,3210CT 6.6% 269.1 / 6.3 91P134P (CCG→CCAgalEUDP‑glucose 4‑epimerase
Reads supporting (aligned to +/- strand):  ref base C (44/41);  new base T (2/4);  total (46/45)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.23e-01

GATGGGTTTCCGGCTCGTGACGCTCGCCTAGCTGGCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGCAGCAGGGAAAGGGTGGCTGCGACGTTGTTCT  >  NC_002947/3542173‑3542469
                                                                                                                                                    |                                                                                                                                                      
gATGGGTTTCCGGCTCGTGACGCTCGCCTAGCTGGCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATAc                                                                                                                                                        >  3:95401/1‑149 (MQ=255)
 aTGGGTTTCCGGCTCGTGACGCTCGCCTAGCTGGCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACAGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACg                                                                                                                                                       >  1:218252/1‑149 (MQ=255)
       ttCCGGCTCGTGACGCTCGCCTAGCTGGCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCGTATGGGtt                                                                                                                                                   >  4:19246/1‑147 (MQ=255)
         ccGGCTCGTGACGCTCGCCTAGCTGGCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCGTATGGGTTGag                                                                                                                                                <  3:19246/148‑2 (MQ=255)
         ccGGCTCGTGACGCTCGCCTAGCTGGCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACgg                                                                                                                                                      <  5:372560/142‑1 (MQ=255)
         ccGGCTCGTGACGCTCGCCTAGCTGGCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACgg                                                                                                                                                      >  6:372560/1‑142 (MQ=255)
          cGGCTCGTGACGCTCGCCTAGCTGGCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCGTATGGGTTGagca                                                                                                                                              <  6:300142/149‑4 (MQ=255)
          cGGCTCGTGACGCTCGCCTAGCTGGCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCGTATGGGTTGagc                                                                                                                                               <  3:193856/148‑3 (MQ=255)
          cGGCTCGTGACGCTCGCCTAGCTGGCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGt                                                                                                                                                    <  3:6667/143‑1 (MQ=255)
          cGGCTCGTGACGCTCGCCTAGCTGGCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGt                                                                                                                                                    >  4:6667/1‑143 (MQ=255)
            gCTCGTGACGCTCGCCTAGCTGGCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCGTATGGGTTGagcaaa                                                                                                                                            >  1:122766/1‑144 (MQ=255)
             cTCGTGACGCTCGCCTAGCTGTCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGcc                                                                                                                                           <  4:214354/149‑1 (MQ=255)
             cTCGTGACGCTCGCCTAGCTGGCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCGTATGGGTTGagcaaat                                                                                                                                           <  2:122766/149‑7 (MQ=255)
                   aCGCTCGCCTAGCTGGCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTg                                                                                                                                     <  5:319823/148‑1 (MQ=255)
                    cGCTCGCCTAGCTGGCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGgt                                                                                                                                   <  4:95401/149‑1 (MQ=255)
                    cGCTCGCCTAGCTGGCCTTCCGGGGCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGgt                                                                                                                                   <  8:347563/149‑1 (MQ=255)
                      cTCGCCTAGCTGGCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGgtgt                                                                                                                                 <  7:8292/149‑1 (MQ=255)
                          ccTAGCTGGCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGATTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGCTGTGTGc                                                                                                                            >  4:223944/1‑149 (MQ=255)
                            tAGCTGGCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCtt                                                                                                                          >  8:111645/1‑149 (MQ=255)
                             aGCTGGCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAACCCATAGGAACGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTCCtt                                                                                                                          >  3:319630/1‑148 (MQ=255)
                             aGCTGGCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCtt                                                                                                                          >  6:107591/1‑148 (MQ=255)
                                  gCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCatcga                                                                                                                    >  2:247044/1‑149 (MQ=255)
                                    cTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATAc                                                                                                                                                        <  4:149010/113‑1 (MQ=255)
                                    cTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATAc                                                                                                                                                        >  3:149010/1‑113 (MQ=255)
                                    cTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCatcgatc                                                                                                                  >  3:278952/1‑149 (MQ=255)
                                     ttCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCatcgatcg                                                                                                                 >  4:111589/1‑149 (MQ=255)
                                               gcgcCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTg                                                                                                       >  2:116472/1‑149 (MQ=255)
                                               gcgcCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTg                                                                                                       >  2:46291/1‑149 (MQ=255)
                                                  cccGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCgg                                                                                                    >  8:170677/1‑149 (MQ=255)
                                                   ccGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGAACGGCCCCTTGGTGTGTGCTTCATAGATAGGAACGTACTACggg                                                                                                   >  7:232267/1‑149 (MQ=255)
                                                   ccGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCACGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCggg                                                                                                   <  1:116472/149‑1 (MQ=255)
                                                         cATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTg                                                                                                >  3:18534/1‑146 (MQ=255)
                                                         cATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTg                                                                                                <  4:18534/146‑1 (MQ=255)
                                                         cATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGcc                                                                                              <  1:46291/148‑1 (MQ=255)
                                                          aTTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCatc                                                                                                                      >  7:70665/1‑123 (MQ=255)
                                                          aTTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCatc                                                                                                                      <  8:70665/123‑1 (MQ=255)
                                                          aTTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCg                                                                                             >  2:212736/1‑148 (MQ=255)
                                                                   gcgcAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATAcggcgg                                                                                   <  1:212736/149‑1 (MQ=255)
                                                                   gcgcAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATAcggcgg                                                                                   <  2:362941/149‑1 (MQ=255)
                                                                      cAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCAc                                                                                <  3:111589/149‑1 (MQ=255)
                                                                         gCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTc                                                                              >  7:261716/1‑148 (MQ=255)
                                                                          cTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCATGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCga                                                                            <  7:111645/149‑1 (MQ=255)
                                                                           ttACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCg                                                                             >  7:245437/1‑147 (MQ=255)
                                                                           ttACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCg                                                                             >  3:140120/1‑147 (MQ=255)
                                                                           ttACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCg                                                                             >  6:43942/1‑147 (MQ=255)
                                                                                gACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGaaaa                                                                      <  4:140120/149‑1 (MQ=255)
                                                                                gACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGaaaa                                                                      <  3:223944/149‑1 (MQ=255)
                                                                                 aatttaaaCCATACGCCCGGTCGAATTCCTCAAGAATCTGCTCGTCCATCCATTTGCTCAACCCATACTGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAac                                                                     <  8:232267/144‑1 (MQ=255)
                                                                                 aCTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTGCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAac                                                                     <  1:240221/149‑1 (MQ=255)
                                                                                 aCTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTACATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGaaaa                                                                      >  1:303673/1‑148 (MQ=255)
                                                                                   ttCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAacac                                                                   >  7:207478/1‑149 (MQ=255)
                                                                                   ttCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAacac                                                                   >  6:235822/1‑149 (MQ=255)
                                                                                   ttCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGATGATCGGCCCCTAGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAaca                                                                    >  6:360884/1‑148 (MQ=255)
                                                                                     cAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATAcggcgg                                                                                   >  8:126326/1‑131 (MQ=255)
                                                                                     cAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATAcggcgg                                                                                   <  7:126326/131‑1 (MQ=255)
                                                                                     cAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAacac                                                                   >  6:92254/1‑147 (MQ=255)
                                                                                     cAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAacac                                                                   >  6:174559/1‑147 (MQ=255)
                                                                                     cAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGGTGCCGTATACGGCGGCACTCGAAGAAAacac                                                                   <  1:247044/147‑1 (MQ=255)
                                                                                      aaaCCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCa                                                                 <  8:245437/148‑1 (MQ=255)
                                                                                       aaCCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCAATACGGGTTGATAGGCCCCTTGGTGTGTGCTTCATAGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAgg                                                               >  1:352390/1‑149 (MQ=255)
                                                                                          cATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGc                                                            >  5:89219/1‑149 (MQ=255)
                                                                                           aTAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGcc                                                           <  5:167247/149‑1 (MQ=255)
                                                                                               gCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACGCGAAGAAAACACCAGGTGCCTGAt                                                       >  4:365298/1‑149 (MQ=255)
                                                                                                 ccGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGc                                                     >  6:334172/1‑149 (MQ=255)
                                                                                                   ggTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGAACTCGAAGAAAACACACGGTGCCTGATGcct                                                   >  8:327611/1‑148 (MQ=255)
                                                                                                       gAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCAt                                               >  4:56562/1‑149 (MQ=255)
                                                                                                        aaGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCGCCTTGGTGTGTGCTTAATAGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCAtt                                              <  4:319630/149‑1 (MQ=255)
                                                                                                        aaGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCAtt                                              <  2:218252/149‑1 (MQ=255)
                                                                                                          gTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGa                                            <  5:235822/149‑1 (MQ=255)
                                                                                                            ccTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGAc                                           <  7:170677/148‑1 (MQ=255)
                                                                                                             cTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCg                                                                                                     <  7:356061/89‑1 (MQ=255)
                                                                                                             cTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCg                                                                                                     >  8:356061/1‑89 (MQ=255)
                                                                                                             cTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCa                                         >  4:351648/1‑149 (MQ=255)
                                                                                                             cTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCa                                         >  5:117927/1‑149 (MQ=255)
                                                                                                                aaGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATg                                       <  6:89219/148‑1 (MQ=255)
                                                                                                                       tGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAg                                                                >  6:261490/1‑116 (MQ=255)
                                                                                                                         cTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTgcagc                              <  8:261716/148‑1 (MQ=255)
                                                                                                                           cGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGCAGCAgg                           <  3:351648/149‑1 (MQ=255)
                                                                                                                            gACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCa                                                                 >  1:125652/1‑110 (MQ=255)
                                                                                                                            gACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCa                                                                 <  2:125652/110‑1 (MQ=255)
                                                                                                                            gACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGCAGCAggg                          <  5:92254/149‑1 (MQ=255)
                                                                                                                            gACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGCAGCAgg                           <  5:174559/148‑1 (MQ=255)
                                                                                                                             aCCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGCAGCAggg                          >  4:110820/1‑148 (MQ=255)
                                                                                                                               catccatTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGCAGCAGGGa                         >  1:76618/1‑148 (MQ=255)
                                                                                                                                   catTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTTCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGCAGCAGGGAAAGGGt                   <  7:327611/149‑1 (MQ=255)
                                                                                                                                   catTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGCAGCAGGGAAAGGGt                   <  3:222466/149‑1 (MQ=255)
                                                                                                                                    atTTGCTCAACCCATACGTGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGCAGCAGGGAAAGGGTg                  <  6:117927/149‑1 (MQ=255)
                                                                                                                                    atTTGCTCAACCCATACGGGTTGATCGGCCCCTTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTACAGCAGGGAAAgg                     >  1:44543/1‑147 (MQ=255)
                                                                                                                                    atTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGCAGCAGGGAAAGGGTg                  <  3:56562/149‑1 (MQ=255)
                                                                                                                                    atTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTACAGCAGGGAAAgg                     <  2:44543/146‑1 (MQ=255)
                                                                                                                                            aaCCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGCAGCAGGGAAAGGGTGGCTGCGa           <  5:26209/148‑1 (MQ=255)
                                                                                                                                               ccATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGCAGCAGGGAAAGGGTGGCTGCGACgtt       <  8:207478/149‑1 (MQ=255)
                                                                                                                                                cATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGAAGCAGGGAAAGGGTCGCTGCGACgttg      >  4:119640/1‑149 (MQ=255)
                                                                                                                                                   aCGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGCCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGCAGCAGGGAAAGGGTGGCTGCGACGTTGTTc   >  2:137009/1‑149 (MQ=255)
                                                                                                                                                    cGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGCAGCAGGGAAAGGGTGGCTGCGACGTTGTTCt  <  5:334172/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
GATGGGTTTCCGGCTCGTGACGCTCGCCTAGCTGGCCTTCCGGGTCAGCGCCCGCTGCATTGAAATAGCGCAGGCTTACCGACTTCAAACCATAGGCCCGGTCGAAGTCCTCAAGAATCTGCTCGACCATCCATTTGCTCAACCCATACGGGTTGATCGGCCCCTTGGTGTGTGCTTCATCGATCGGCACGTACTGCGGGTTGCCGTATACGGCGGCACTCGAAGAAAACACCAGGTGCCTGATGCCGGCATTGACCATGGCCTGCAGCAGGGAAAGGGTGGCTGCGACGTTGTTCT  >  NC_002947/3542173‑3542469

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: