Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,019,865 T→C 6.3% R120R (CGT→CGC PP_2638 → cellulose synthase operon protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,019,8650TC 6.3% 232.8 / 4.2 79R120R (CGT→CGCPP_2638cellulose synthase operon protein C
Reads supporting (aligned to +/- strand):  ref base T (35/39);  new base C (3/2);  total (38/41)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.87e-01

GAGAAGAACAGCAGTGAAGCCGAGCAGTTGAGTGCGCGGATCAGTGCCGTGGCGCCTGGCAGCCTGCAAGCACGTCAGGCCGCACGCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGAACAGTTGTTTGCCGGCGAACCACCCAGCCTTGAATTGGCCCTGGACTACTGGCGCGTGCGCGGCAATCTGCCAGGCCAACGCCCTGAGGCGATCCGCCA  >  NC_002947/3019728‑3020011
                                                                                                                                         |                                                                                                                                                  
gagaAGAACAGCAGTGAAGCCGAGCAGTTGAGTGCGCGGATCAGTGCCGTGGCGCCTGGCAGCCTGCAAGCACGTCAGGCCGCACGCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGtt                                                                                                                                        <  1:1482315/150‑1 (MQ=255)
gagaAGAACAGCAGTGAAGCCGAGCAGTGGAGTGCGCGGATCAGTGCCGTGGCGCCTGGCAGCCTGCAAGCACGTCAGGCCGCACGCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCTCGTCTGCTGGCAGtt                                                                                                                                        <  1:1716763/150‑1 (MQ=255)
   aagaaCAGCAGTGAAGCCGAGCAGTTGAGTGCGCGGATCAGTGCCGTGGCGCCTGGCAGCCTGCAAGCACGTCATGCCGCACGCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTAcc                                                                                                                                     >  2:1576861/1‑150 (MQ=255)
   aagaaCAGCAGTGAAGCCGAGCAGTTGAGTGCGCGGATCAGTGCCGTGGCGCCTGGCAGCCTGCAAGCACGTCAGGCCGCACGCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTAcc                                                                                                                                     >  2:1343530/1‑150 (MQ=255)
     gaaCAGCAGTGAAGCCGAGCAGTTGAGTGCGCGGATCAGTGCCGTGGCGCCTGGCAGCCTGCAAGCACGTCAGGCCGCACGCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCgg                                                                                                                                   >  1:379898/1‑150 (MQ=255)
         agcagTGAAGCCGAGCAGTTGAGTGCGCGGATCAGTGCCGTGGCGCCTGGCAGCCTGCAAGCACGTCAGGCCGCACGCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGcc                                                                                                                                 >  2:877533/1‑148 (MQ=255)
         agcagTGAAGCCGAGCAGTTGAGTGCGCGGATCAGTGCCGTGGCGCCTGGCAGCCTGCAAGCACGTCAGGCCGCACGCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGcc                                                                                                                                 <  1:877533/148‑1 (MQ=255)
         agcagTGAAGCCGAGCAGTTGAGTGCGCGGATCAGTGCCGTGGCGCCTGGCAGCCTGCAAGCACGTCAGGCCGCACGCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGt                                                                                                                               >  2:725605/1‑150 (MQ=255)
           cagTGAAGCCGAGCAGTTGAGTGCGCGGATCAGTGCCGTGGCGCCTGGCAGCCTGCAAGCACGTCAGGCCGCACGCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTaa                                                                                                                             <  1:1223942/150‑1 (MQ=255)
               gAAGCCGAGCAGTTGAGTGCGCGGATCAGTGCCGTGGCGCCTGGCAGCCTGCAAGCACGTCAGGCCGCACGCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGaa                                                                                                                         <  2:379898/150‑1 (MQ=255)
                aaGCCGAGCAGTTGAGTGCGCGGATCAGTGCCGTGGCGCCTGGCAGCCTGCAAGCACGTCAGGCCGCACGCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAag                                                                                                                        >  2:1659642/1‑150 (MQ=255)
                    cGAGCAGTTGAGTGCGCGGATCAGTGCCGTGGCGCCTGGCAGCCTGCAAGCACGTCAGGTCGCACGCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGgc                                                                                                                    >  2:1698705/1‑150 (MQ=255)
                           ttGAGTGCGCGGATCAGTGCCGTGGCGCCTGGCAGCCTGCAAGCACGTCAGGCCGCACGCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCg                                                                                                             <  1:924397/150‑1 (MQ=255)
                             gAGTGCGCGGATCAGTGCCGTGGCGCCTGGCAGCCTGCAAGCACGTCAGGCCGCACGCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGt                                                                                                           <  1:694507/150‑1 (MQ=255)
                             gAGTGCGCGGATCAGTGCCGTGGCGCCTGGCAGCCTGCAAGCACGTCAGGCCGCACGCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGt                                                                                                           <  1:1936878/150‑1 (MQ=255)
                              aGTGCGCGGATCAGTGCCGTGGCGCCTGGCAGCCTGCAAGCACGTCAGGCCGCACGCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTGATGAAGAAGCGCTGGCGGta                                                                                                          >  2:2086267/1‑150 (MQ=255)
                                 gcgcGGATCAGTGCCGTGGCGCCTGGCAGCCTGCAAGCACGTCAGGCCGCACGCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGACGCGCctgctgcaa                                                                                                                                           >  2:1100169/1‑111 (MQ=255)
                                 gcgcGGATCAGTGCCGTGGCGCCTGGCAGCCTGCAAGCACGTCAGGCCGCACGCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGACGCGCctgctgcaa                                                                                                                                           <  1:1100169/114‑4 (MQ=255)
                                   gcgGATCAGTGCCGTGGCGCCTGGCAGCCTGCAAGCACGTCAGGCCGCACGCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGa                                                                                                     <  2:2036190/150‑1 (MQ=255)
                                     ggATCAGTGCCGTGGCGCCTGGCAGCCTGCAAGCACGTCAGGCCGCACGCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGaa                                                                                                    >  2:2150628/1‑149 (MQ=255)
                                     ggATCAGTGCCGTGGCGCCTGGCAGCCTGCAAGCACGTCAGGCCGCACGCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGaa                                                                                                    <  1:2150628/149‑1 (MQ=255)
                                     ggATCAGTGCCGTGGCGCCTGGCAGCCTGCAAGCACGTCAGGCCGCACGCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGAAc                                                                                                   >  1:1031752/1‑150 (MQ=255)
                                      gATCAGTGCCGTGGCGCCTGGCAGCCTGCAAGCACGTCAGGCCGCACGCCTGCTTGAGCTTCAAGAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGAACa                                                                                                  <  2:938033/150‑1 (MQ=255)
                                       aTCAGTGCCGTGGCGCCTGGCAGCCTGCAAGCACGTCAGGCCGCACGCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGAACAg                                                                                                 >  2:2268998/1‑150 (MQ=255)
                                         cAGTGCCGTGGCGCCTGGCAGCCTGCAAGCACGTCAGGCCGCACGCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGAACAgtt                                                                                               <  2:1011796/150‑1 (MQ=255)
                                           gTGCCGTGGCGCCTGGCAGCCTGCAAGCACGTCAGGCCGCACGCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGAACAgttgt                                                                                             >  1:917945/1‑150 (MQ=255)
                                            tGCCGTGGCGCCTGGCAGCCTGCAAGCACGTCAGGCCGCACGCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTaa                                                                                                                             >  1:210434/1‑117 (MQ=255)
                                            tGCCGTGGCGCCTGGCAGCCTGCAAGCACGTCAGGCCGCACGCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTaa                                                                                                                             <  2:210434/117‑1 (MQ=255)
                                              ccGTGGCGCCTGGCAGCCTGCAAGCACGTCAGGCCGCACGCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGAACAGTTGTTTg                                                                                          <  1:725605/150‑1 (MQ=255)
                                              ccGTGGCGCCTGGCAGCCTGCAAGCACGTCAGGCCGCACGCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGAACAGTTGTTTg                                                                                          >  2:1577123/1‑150 (MQ=255)
                                                 tGGCGCCTGGCAGCCTGCAAGCACGTCAGGCCGCACGCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGAACAGTTGTTTGCCg                                                                                       >  1:1920185/1‑150 (MQ=255)
                                                         ggCAGCCTGCAAGCACGTCAGGCCGCACGCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGAACAGTTGTTTGCCGGCGAacca                                                                               <  1:2268998/150‑1 (MQ=255)
                                                              ccTGCAAGCACGTCAGGCCGCACGCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGAACAGTTGTTTGCCGGCGAACCACCCAg                                                                          <  1:380118/150‑1 (MQ=255)
                                                                tGCAAGCACGTCAGGCCGCACGCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCg                                                                                                                                    <  1:1467089/90‑1 (MQ=255)
                                                                tGCAAGCACGTCAGGCCGCACGCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCg                                                                                                                                    >  2:1467089/1‑90 (MQ=255)
                                                                tGCAAGCACGTCAGGCCGCACGCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGAACAGTTGTTTGCCGGCGAACCACCCAGcc                                                                        >  2:1279353/1‑150 (MQ=255)
                                                                  cAAGCACGTCAGGTCGCACGCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGAACAGTTGTTTGCCGGCGAACCACCCAGCCtt                                                                      <  1:1698705/150‑1 (MQ=255)
                                                                     gCACGTCAGGCCGCACGCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGa                                                                                                                          <  2:1607491/95‑1 (MQ=255)
                                                                     gCACGTCAGGCCGCACGCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGa                                                                                                                          >  1:443561/1‑95 (MQ=255)
                                                                     gCACGTCAGGCCGCACGCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGa                                                                                                                          <  2:443561/95‑1 (MQ=255)
                                                                     gCACGTCAGGCCGCACGCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGa                                                                                                                          >  1:1607491/1‑95 (MQ=255)
                                                                      cACGTCAGGCCGCACGCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGAACAGTTGTTTGCCGGCGAACCACCCAGCCTTGAAt                                                                  <  1:1866890/150‑1 (MQ=255)
                                                                       aCGTCAGGCCGCACGCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGAACAGTTGTTTGCCGGCGAACCACCCAGCCTTGAAtt                                                                 >  2:1471918/1‑150 (MQ=255)
                                                                       aCGTCAGGCCGCACGCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGAACAGTTGTTTGCCGGCGAACCACCCAGCCTTGAAtt                                                                 >  1:366388/1‑150 (MQ=255)
                                                                          tCAGGCCGCACGCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGAACAGTTGTTTGCCGGCGAACCACCCAGCCTTGAATTGGc                                                              <  2:596921/150‑1 (MQ=255)
                                                                            aGGCCGCACGCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCTCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGAACAGTTGTTTGCCGGCGAACCACCCAGCCTTGAATTGGccc                                                            <  2:1833418/150‑1 (MQ=255)
                                                                              gCCGCACGCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGAACAGTTGTTTGCCGGCGAACCACCCAGCCTTGAATTGGCCCTg                                                          >  2:1266121/1‑150 (MQ=255)
                                                                                cGCACGCCTTCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGAACAGTTGTTTGCCGGCGAACCACCCAGCCTTGAATTGGCCCTGGa                                                        <  1:234868/150‑1 (MQ=255)
                                                                                cGCACGCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGAACAGTTGTTTGCCGGCGAac                                                                                 <  2:1499922/125‑1 (MQ=255)
                                                                                cGCACGCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGAACAGTTGTTTGCCGGCGAac                                                                                 >  1:1499922/1‑125 (MQ=255)
                                                                                     gCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGAACAGTTGTTTGCCg                                                                                       >  2:1199091/1‑114 (MQ=255)
                                                                                     gCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGAACAGTTGTTTGCCg                                                                                       <  1:1199091/114‑1 (MQ=255)
                                                                                     gCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGAACAGTTGTTTGCCGGCGAACCACCCAGCCTTGAATTGGcc                                                             <  2:1417422/140‑1 (MQ=255)
                                                                                     gCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGAACAGTTGTTTGCCGGCGAACCACCCAGCCTTGAATTGGcc                                                             >  1:1417422/1‑140 (MQ=255)
                                                                                     gCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGAACAGTTGTTTGCCGGCGAACCACCCAGCCTTGAATTGGCCCTGGactact                                                   <  1:1471918/150‑1 (MQ=255)
                                                                                     gCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGAACAGTTGTTTGCCGGCGAACCACCCAGCCTTGAATTGGCCCTGGa                                                        <  2:1712354/145‑1 (MQ=255)
                                                                                     gCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGAACAGTTGTTTGCCGGCGAACCACCCAGCCTTGAATTGGCCCTGGa                                                        >  1:1712354/1‑145 (MQ=255)
                                                                                     cccTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGAACAGTTGTTTGCCGGCGAACCACCCAGCCTTGAATTGGCCCTGGactact                                                   >  2:1603487/2‑150 (MQ=255)
                                                                                            tGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGAACAGTTGTTTGCCGGCGAACCACCCAGCCTTGAATTGGCCCTGGACTACTGGCGCGt                                            <  1:381721/150‑1 (MQ=255)
                                                                                                cTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGGAATGAAGAGGCGCTGGCGGTATACGAACAGTTGTTTGCCGGCGAACCACCCAGCCCGGAATTGGCCCTGGACTACTAGGGAGTgcgc                                        >  1:1065229/1‑150 (MQ=255)
                                                                                                    aacaacCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGAACAGTTGTTTGCCGGCGAACCACCCAGCCTTGAATTGGc                                                              <  1:793387/124‑1 (MQ=255)
                                                                                                    aacaacCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGAACAGTTGTTTGCCGGCGAACCACCCAGCCTTGAATTGGc                                                              >  2:793387/1‑124 (MQ=255)
                                                                                                            gcgcCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGAACAGTTGTTTGCCGGCGAACCACCCAGCCTTGAATTGGCCCTGGACTACTGGCGCGTGCGCGGCAATCTGCCa                            <  1:1603487/150‑1 (MQ=255)
                                                                                                              gcCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGAACAGTTGTTTGCCGGCGAACCACCCAGCCTTGAATTGGCCCTGGACTACTGGCGCGTGCGCGGCAATCTGCCAgg                          <  1:2078497/150‑1 (MQ=255)
                                                                                                                cTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGAACAGTTGTTTGCCGGCGAACCACCCAGCCTTGAATTGGCCCTGGACTACTGGCGCGTGCGCGGCAATCTGCCAGGcc                        >  2:1825798/1‑150 (MQ=255)
                                                                                                                cTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGAACAGTTGTTTGCCGGCGAACCACCCAGCCTTGAATTGGCCCTGGACTACTGGCGCGTGCGCGGCAATCTGCCAGGcc                        >  2:391966/1‑150 (MQ=255)
                                                                                                                 tGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGAACAGTTGTTTGCCGGCGAACCACCCAGCCTTGAATTGGCCCTGGACTACTGGCGCGTgcg                                         <  2:100007/132‑1 (MQ=255)
                                                                                                                 tGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGAACAGTTGTTTGCCGGCGAACCACCCAGCCTTGAATTGGCCCTGGACTACTGGCGCGTgcg                                         >  1:100007/1‑132 (MQ=255)
                                                                                                                    aCGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGAACAGTTGTTTGCCGGCGAACCACCCAGCCTTGAATTGGCCCTGGACTACTGGCGCGTGCGCGGCAATCTGCCAGGCCAACg                    >  1:288030/1‑150 (MQ=255)
                                                                                                                      ggCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGAACAGTTGTTTGCCGGCGAACCACCCAGCCTTGAATTGGCCCTGGACTACTGGCGCGTGCGCGGCAATCTGCCAGGCCAACGcc                  <  2:1031752/150‑1 (MQ=255)
                                                                                                                         gcTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGAACAGTTGTTTGCCGGCGAACCACCCAGCCTTGAATTGGCCCTGGACTACTGGCGCGTGCGCGGCAATCTGCCAGGCCAACGCCCTg               <  1:391966/150‑1 (MQ=255)
                                                                                                                          cTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGAACAGTTGTTTGCCGGCGAACCACCCAGCCTTGAATTGGCCCTGGACTACTGGCGCGTGCGCGGCAATCTGCCAGGCCAACGCCCTGa              <  1:1343530/150‑1 (MQ=255)
                                                                                                                           tgccagcagACGCGCCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGAACAGTTGTTTGCCGGCGAACCACCCAGCCTTGAATTGGCCCTGGACTACTGGCGCGTGCGCGGCAATCTGCCAGGCCAACGCCCTGAg             >  2:1860791/4‑150 (MQ=255)
                                                                                                                            gccagcagACGCGCCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGAACAGTTGTTTGCCGGCGAACCACCCAGCCTTGAATTGGCCCTg                                                          <  1:1549343/102‑1 (MQ=255)
                                                                                                                            gccagcagACGCGCCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGAACAGTTGTTTGCCGGCGAACCACCCAGCCTTGAATTGGCCCTg                                                          >  2:1549343/3‑104 (MQ=255)
                                                                                                                               agcagGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGAACAGTTGTTTGCCGGCGAACCACCCAGCCTTGAATTGGCCCTGGACTACTGGCGCGTGCGCGGCAATCTGCCAGGCCAACGCCCTGAGGCGa         >  2:1613503/1‑150 (MQ=255)
                                                                                                                                    gcgcGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGAACAGTTGTTTGCCGGCGAACCACCCAGCCTTGAATTGGCCCTGGACTACTGGCGCGTGCGCGGCAATCTGCCAGGCCAACGCCCTGAGGCGATCCGc    >  2:1610569/1‑150 (MQ=255)
                                                                                                                                    gcgcGTCTGCTGGCAGTTAACGGCCGTAATGAAGCGGCGCTGGCGGTATACGAACAGTTGTTTGCCGGCGAACCACCCAGCCTTGAATTGGCCCTGGACTACTGGCGCGTGCGCGGCAATCTGCCAGGCCAACGCCCTGAGGCGATCCGc    <  2:1092204/150‑1 (MQ=255)
                                                                                                                                      gcgTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGAACAGTTGTTTGCCGGCGAACCACCCAGCCTTGAATTGGCCCTGGACTACTGGCGCGTGCGCGGCAATCTGCCAGGCCAACGCCCTGAGGCGATCCGCCa  <  1:1659642/150‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                  
GAGAAGAACAGCAGTGAAGCCGAGCAGTTGAGTGCGCGGATCAGTGCCGTGGCGCCTGGCAGCCTGCAAGCACGTCAGGCCGCACGCCTGCTTGAGCTTCAACAACCGGCGCCTGCACGGCGCTTGCAGCAGGCGCGTCTGCTGGCAGTTACCGGCCGTAATGAAGAGGCGCTGGCGGTATACGAACAGTTGTTTGCCGGCGAACCACCCAGCCTTGAATTGGCCCTGGACTACTGGCGCGTGCGCGGCAATCTGCCAGGCCAACGCCCTGAGGCGATCCGCCA  >  NC_002947/3019728‑3020011

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