Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,873,799:1 +T 7.1% intergenic (+62/+20) cigR → / ← PP_2527 membrane protein/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,873,7991.T 7.1% 380.2 / 24.0 98intergenic (+62/+20)cigR/PP_2527membrane protein/membrane protein
Reads supporting (aligned to +/- strand):  ref base . (48/43);  new base T (2/5);  total (50/48)
Fisher's exact test for biased strand distribution p-value = 2.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.05e-01

Only 100 of 107 total aligned reads displayed.

AGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGCACCACAGCGCACAGCATGGCCAGCCCGATGGTCACGCCAACGCCCCACGTTTCCAGCAGCCCTGCCAATACCAGATAGAACGCAATCAC  >  NC_002947/2873652‑2873947
                                                                                                                                                    |                                                                                                                                                    
aGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTGTGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTCCGATGGGCTgcgc                                                                                                                                                    >  2:1321042/1‑150 (MQ=255)
    gATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCgcagca                                                                                                                                                <  2:2796513/150‑1 (MQ=255)
    gATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCgcagca                                                                                                                                                <  1:1978798/150‑1 (MQ=255)
    gATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCgcagca                                                                                                                                                <  2:1107403/150‑1 (MQ=255)
    gATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGCCGCTGGCGTGCCTGCGATGGGCTGCgcagca                                                                                                                                                >  2:574250/1‑150 (MQ=255)
    gATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGAGCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCgcagca                                                                                                                                                <  1:1098245/150‑1 (MQ=255)
    gATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGCGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCATGCGCTGGCGTGCCTGCGATGGGCTGCgcagca                                                                                                                                                <  1:1491816/150‑1 (MQ=255)
          aGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCaa                                                                                                                                          >  1:1319842/1‑150 (MQ=255)
           gCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGCCTAAAGAAGCAGGTGCGCGCAGTCCTCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCTGCGCAGCCCatc                                                                                                                                         <  1:375302/150‑3 (MQ=255)
           gCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGTTGGGCTGCTGCGCAGCCCatc                                                                                                                                         <  2:1367771/150‑3 (MQ=255)
           gCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAAc                                                                                                                                         >  1:1366641/1‑150 (MQ=255)
           gCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAAc                                                                                                                                         <  2:6375/150‑1 (MQ=255)
             gggCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACgg                                                                                                                                       >  1:2109811/1‑150 (MQ=255)
               gcgcAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGTTGGGCTGCTGCGCAGCCCatcgcag                                                                                                                                     <  2:1562801/150‑7 (MQ=255)
                 gcAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTa                                                                                                                                   <  1:773721/150‑1 (MQ=255)
                      ccTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACTGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAAc                                                                                                                              <  1:1203313/150‑1 (MQ=255)
                       atGATTATGGCCGCCCCGGCTGCCGGGAGCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACg                                                                                                                             >  2:360025/2‑150 (MQ=255)
                           ttCTGGTCGCCGTGTCGACCGGTATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGc                                                                                                                         >  2:2349017/1‑150 (MQ=255)
                            tCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGcc                                                                                                                        <  1:2751873/150‑1 (MQ=255)
                              tGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGt                                                                                                                      >  2:202394/1‑150 (MQ=255)
                                   gCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCTGCGCAGCCCatc                                                                                                                                         >  2:2591839/1‑124 (MQ=255)
                                   gCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCTGCGCAGCCCatc                                                                                                                                         <  1:2591839/126‑3 (MQ=255)
                                      gTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTACGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAATGGCCTATCAACGCTGCCGTAACACcgc                                                                                                              <  1:1440102/150‑1 (MQ=255)
                                       tGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACcgcg                                                                                                             <  2:601027/150‑1 (MQ=255)
                                           gACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGTTGGGCTGCTGCGCAGCCCatcgcag                                                                                                                                     >  2:2437183/1‑116 (MQ=255)
                                              cGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGcc                                                                                                      >  1:180450/1‑150 (MQ=255)
                                                 gATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCTAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCaa                                                                                                   <  1:2690166/150‑1 (MQ=255)
                                                   tcatcTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAgg                                                                                                 >  1:1681967/1‑150 (MQ=255)
                                                        tATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAg                                                                                            >  2:1547103/1‑150 (MQ=255)
                                                        tATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAg                                                                                            <  2:2630371/150‑1 (MQ=255)
                                                          tGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAg                                                                                                  >  2:1806353/1‑142 (MQ=255)
                                                          tGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAg                                                                                                  <  1:1806353/142‑1 (MQ=255)
                                                          tGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGca                                                                                          >  1:309681/1‑150 (MQ=255)
                                                            aGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGcacc                                                                                        >  1:1526790/1‑150 (MQ=255)
                                                             ggTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGcc                                                                                                                        <  2:2482589/117‑1 (MQ=255)
                                                             ggTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGcc                                                                                                                        >  1:2482589/1‑117 (MQ=255)
                                                                gtTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGCACCACAg                                                                                    >  2:1691761/1‑150 (MQ=255)
                                                                 ttGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGCACCACAgc                                                                                   >  2:1292530/1‑150 (MQ=255)
                                                                 ttGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGCACCACAgc                                                                                   <  1:134022/150‑1 (MQ=255)
                                                                 ttGAACAGGGCATGCAACTAAGGCCGCTGGGGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCCACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGCACCACAgc                                                                                   >  2:2491688/1‑150 (MQ=255)
                                                                     aCGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGCACCACAGCGcac                                                                               <  2:1681967/150‑1 (MQ=255)
                                                                     aCGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCACGGCCAGCCCAAGGGTCAGCACCACAGCGcac                                                                               >  1:1965455/1‑150 (MQ=255)
                                                                      cGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGCACCACAGCGcaca                                                                              <  1:360025/150‑1 (MQ=255)
                                                                      cGGGGCATTCGACTAAGGCCGCAGGTGCAAGCTGCCCGCAGGGTCCGTTTTTGCGCTGGCGTGCCTGCGATGGTTTTCTCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGCACCACAGCGcaca                                                                              <  1:1452952/150‑1 (MQ=255)
                                                                       ggggCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGCACCACAGCGCACAg                                                                             >  1:2566918/1‑150 (MQ=255)
                                                                       ggggCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGCACCACAGCGCACAg                                                                             <  1:1292530/150‑1 (MQ=255)
                                                                        gggCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGCACCACAGCGCACAGc                                                                            <  2:309681/150‑1 (MQ=255)
                                                                            aTTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGCACCACAGCGCACAGCATgg                                                                        >  1:669530/1‑150 (MQ=255)
                                                                                       gCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGCACCACAGCGCACAGCATGGCCAGCCCGATg                                                             <  2:180450/150‑1 (MQ=255)
                                                                                        ccGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGCACCACAGCGCACAGCATGGCCAGCCCGATgg                                                            >  1:111909/1‑150 (MQ=255)
                                                                                            aGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGCACCACAGCGCACAGCATGGCCAGCCCGATGGTCAc                                                        <  2:1526790/150‑1 (MQ=255)
                                                                                             ggTGCGAGCAGGCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGCACCACAGCGCACAGCATGGCCAGCCCGATGGTCACg                                                       <  2:1136954/150‑1 (MQ=255)
                                                                                              gTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGCACCACAGCGCACAGCATGGCCAGCCCGATGGTCACGc                                                      >  2:841298/1‑150 (MQ=255)
                                                                                               tGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGCACCACAGCGCACAGCATGGCCAGCCCGATGGTCACGcc                                                     >  1:1904352/1‑150 (MQ=255)
                                                                                                gCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAacac                                                                                                                 <  1:10303/89‑1 (MQ=255)
                                                                                                gCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAacac                                                                                                                 >  2:10303/1‑89 (MQ=255)
                                                                                                gCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGCACCACAGCGCACAGCATGGCCAGCCCGATGGTCACGCCa                                                    <  2:1904352/150‑1 (MQ=255)
                                                                                                 cGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGCACCACAGCGCACAGCATGGCCAGCCCGATGGTCACGCCaa                                                   <  2:111909/150‑1 (MQ=255)
                                                                                                    gCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCACGGCCAGCCCAAGGGTCAGCACCACAGCGCACAGCATGGCCAGCCCGATGGTCACGCCAACGc                                                <  2:1965455/150‑1 (MQ=255)
                                                                                                    gCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGATGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGCACCACAGCGCACAGCATGGCCAGCCCGATGGTCACGCCAACGc                                                >  2:1849399/1‑150 (MQ=255)
                                                                                                      aGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGCACCACAGCGCACAGCATGGCCAGCCCGATGGTCACGCCAACGccc                                              <  2:909015/150‑1 (MQ=255)
                                                                                                        tCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGCACCACAgc                                                                                   >  1:1591876/1‑111 (MQ=255)
                                                                                                        tCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGCACCACAgc                                                                                   <  2:1591876/111‑1 (MQ=255)
                                                                                                        tCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGCACCACAGCGCACAGCATGGCCAGCCCGATGGTCACGCCAACGCCCCa                                            >  2:3028110/1‑150 (MQ=255)
                                                                                                        tCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGCACCACAGCGCACAGCATGGCCAGCCCGATGGTCACGCCAACGCCCCa                                            >  1:846290/1‑150 (MQ=255)
                                                                                                           gCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCa                                                                                             <  1:1693670/98‑1 (MQ=255)
                                                                                                           gCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCa                                                                                             >  2:1693670/1‑98 (MQ=255)
                                                                                                           gCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCAGTAACAACGCGGCCAGCCCAAGGGTCAGCACCACAGCGCACAGCATGGCCAGCCCGATGGTCACGCCAACGCCCCACGt                                         >  1:370911/1‑150 (MQ=255)
                                                                                                              cGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGCACCACAGCGCACAGCACGGCCAGCCCGATGGTCACGCCAACGCCCCACGTTTc                                      <  1:798320/150‑1 (MQ=255)
                                                                                                               ggTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGCACCACAGCGCACAGCATGGCCAGCCCGATGGTCACGCCAACGCCCCACGTTTcc                                     >  1:93138/1‑150 (MQ=255)
                                                                                                                gTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCa                                                                                             >  2:426239/1‑93 (MQ=255)
                                                                                                                gTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCa                                                                                             <  1:426239/93‑1 (MQ=255)
                                                                                                                 tCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGcacca                                                                                       <  1:2301401/98‑1 (MQ=255)
                                                                                                                 tCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGcacca                                                                                       >  2:2301401/1‑98 (MQ=255)
                                                                                                                   ccTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGCACCACAGCGCACAGCATGGCCAGCCCGATGGTCACGCCCACGCCCCACGTTTCcagca                                 >  1:1080082/1‑150 (MQ=255)
                                                                                                                   ccTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGCACCACAGCGCACAGCATGGCCAGCCCGATGGTCACGCCAACGCCCCACGTTTCcagca                                 <  1:1691761/150‑1 (MQ=255)
                                                                                                                    cTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGcc                                                                                                                                     >  2:472123/1‑49 (MQ=255)
                                                                                                                    cTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGcc                                                                                                                                     <  1:472123/49‑1 (MQ=255)
                                                                                                                    cTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGCACCACAGCGCACAGCATGGCCAGCCCGATGGTCACGCCAACGCCCCACGTTTCcagcag                                >  1:1521351/1‑150 (MQ=255)
                                                                                                                     tGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGCACCACAGCGCACAGCATGGCCAGCCCGATGGTCACGCCAACGCCCCACGTTTCCagcagc                               >  1:775018/1‑150 (MQ=255)
                                                                                                                     tGCAGGCGCTGGCGGGCCTGCGATGGGCTGCGCAGCAGCCCAACGGACTGGCTGCGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGCACCACAGCGCACAGCATGGCCAGCCCGATGGTCACGCCAACGCCCCACGTTTCCagcagc                               <  1:2491688/150‑1 (MQ=255)
                                                                                                                          gcgcTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGCACCACAGCGCACAGCATGGCCAGCCCGATGGTCACGCCAACGCCCCACGTTTCCAGCAGCCCTGc                          <  2:1319842/150‑1 (MQ=255)
                                                                                                                           cgcTGGCGTGCCTGCGATGGGCTGCGCAGCAGCACAACGGGCTTGCAACGCTGCCGTAACACCGCGGCCCGCCCAATTGTCAGCACCACAGCGCACAGCATGGCCAGCCCGATGGTCACGCCAACGCCCCACGTTTCCAGCAGCCCTGcc                         <  1:202394/150‑1 (MQ=255)
                                                                                                                            gcTGGCGCGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTACCGTAACACCGCGGCCAGCCGAAGGGTCAGCACCACAGCGCACAGCATGGCCAGCCCGATGGTCACGCCAACGCCCCACGTTTCCAGCAGCCCTGCCa                        >  2:2208000/1‑150 (MQ=255)
                                                                                                                              tGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGCACCACAGCGCACAGCATGGCCAGCCCGATGGTCACGCCAACGCCCCCCGTTTCCAGCAGCCCTGCCAAt                      >  2:687200/1‑150 (MQ=255)
                                                                                                                                  gTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGCACCACAGCGCACAGCATGGCCAGCCCGATGGTCACGCCAACGCCCCACGTTTCCAGCAGCCCTGCCAATACCa                  <  2:1648341/150‑1 (MQ=255)
                                                                                                                                    gCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGTGTCAGCACCACAGCGCACAGCATGGCCAGCCCGATGGTCACGCCAACGCCCCACGTTTCCAGCAGCCCTGCCAATACCAGa                <  1:1591721/150‑1 (MQ=255)
                                                                                                                                    gCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGCACCACAGCGCACAGCATGGCCAGCCCGATGGTCACGCCAACGCCCCACGTTTCCAGCAGCCCTGCCAATACCAGa                <  2:1689263/150‑1 (MQ=255)
                                                                                                                                      cTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGCACCACAGCGCACAGCATGGCCAGCCCGATg                                                             <  1:711603/103‑1 (MQ=255)
                                                                                                                                      cTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGCACCACAGCGCACAGCATGGCCAGCCCGATg                                                             >  2:711603/1‑103 (MQ=255)
                                                                                                                                      cTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGCACCACAGCGCACAGCATGGCCAGCCCGATGGTCACGCCAACGCCCCACGTTTCCAGCAGCCCTGCCAATACCAGATa              >  1:1764811/1‑150 (MQ=255)
                                                                                                                                       tGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGCACCACAGCGCACAGCATGGCCAGCCCGATGGTCACGCCAACGCCCCACGTTTCCAGCAGCCCTGCCAATACCAGATAg             >  1:149766/1‑150 (MQ=255)
                                                                                                                                       tGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGCACCACAGCGCACAGCATGGCCAGCCCGATGGTCACGCCAACGCCCCACGTTTCCAGCAGCCCTGCCAATACCAGATAg             >  2:3008891/1‑150 (MQ=255)
                                                                                                                                          gATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGCACCACAGCGCACAGCATGGCCAGCCCGATGGTCACGCCAACGCCCCACGTTTCCAGCAGCCCTGCCAATACCAGATAGAAc          >  1:2988181/1‑150 (MQ=255)
                                                                                                                                          aaTGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGCACCACAGCGCACAGCATGGCCAGCCCGATg                                                             >  2:2755623/2‑99 (MQ=255)
                                                                                                                                          aaTGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGCACCACAGCGCACAGCATGGCCAGCCCGATg                                                             <  1:2755623/98‑1 (MQ=255)
                                                                                                                                              ggCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGCACCACAGCGCACAGCATGGCCAGCCCGATGGTCACGCCAACGCCCCACGTTTCCAGCAGCCCTGCCAATACCAGATAGAACGCaa      <  2:1764811/150‑1 (MQ=255)
                                                                                                                                                cTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGCACCACAGCGCACAGCATGGCCAGCCCGATGGTCACGCCAACGCCCCACGTTTCCAGCAGCCCTGCCAATACCAGATAGAACGCAATc    >  2:771716/1‑150 (MQ=255)
                                                                                                                                                  gcgcAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGCACCACAGCGCACAGCATGGCCAGCCCGATGGTCACGCCAACGCCCCACGTTTCCAGCAGCCCTGCCAATACCAGATAGAACGCAATCAc  <  2:149766/150‑1 (MQ=255)
                                                                                                                                                  gcgcAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGCACCACAGCGCACAGCATGGCCAGCCCGATGGTCACGCCAACGCCCCACGTTTCCAGCAGCCCTGCCAATACCAGATAGAACGCAATCAc  <  1:841298/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGTGGATGCAAGCGGGCGCAGACCTGATTCTGGTCGCCGTGGCGACCGGGATCATCTATGAGGTGTTGAACGGGGCATTCGACTAAGGCCGCAGGTGCGAGCAGTCCGCTCGGTCCCTGCAGGCGCTGGCGTGCCTGCGATGGGCTGCGCAGCAGCCCAACGGCCTATCAACGCTGCCGTAACACCGCGGCCAGCCCAAGGGTCAGCACCACAGCGCACAGCATGGCCAGCCCGATGGTCACGCCAACGCCCCACGTTTCCAGCAGCCCTGCCAATACCAGATAGAACGCAATCAC  >  NC_002947/2873652‑2873947

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: