Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,934,031 A→T 9.3% C476S (TGC→AGC)  PP_3467 ← sigma‑54 dependent transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,934,0310AT 9.3% 284.0 / 16.3 96C476S (TGC→AGC) PP_3467sigma‑54 dependent transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (41/46);  new base T (4/5);  total (45/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.42e-01

CTTGTCCAGCAGGCCCTGGATGATCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAGGAACAGCGTGCCACCATGGGCCTGCATC  >  NC_002947/3933882‑3934179
                                                                                                                                                     |                                                                                                                                                     
cTTGTCCAGCAGGCCCTGGATGATCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCa                                                                                                                                                       <  1:784745/150‑1 (MQ=255)
cTTGGCCAGCAGGCCCTGGATGATCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCa                                                                                                                                                       >  2:1603821/1‑150 (MQ=255)
 ttGTCCAGCAGGCCCTGGATGATCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAg                                                                                                                                                      >  1:667912/1‑150 (MQ=255)
   gTCCAGCAGGCCCTGGATGATCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGAt                                                                                                                                                    <  1:1232450/150‑1 (MQ=255)
       agcagGCCCTGGATGATCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGATTAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAgc                                                                                                                                                <  1:319156/150‑1 (MQ=255)
       agcagGCCCTGGATGATCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAgc                                                                                                                                                <  1:318226/150‑1 (MQ=255)
       agcagGCCCTGGATGATCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAgc                                                                                                                                                <  2:1558195/150‑1 (MQ=255)
           ggCCCTGGATGATCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCTGATCTGCGCCa                                                                                                                                            <  1:1590546/150‑1 (MQ=255)
              ccTGGATGATCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGt                                                                                                                                         >  1:756166/1‑150 (MQ=255)
               cTGGATGATCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTc                                                                                                                                        <  1:1103493/150‑1 (MQ=255)
               cTGGATGATCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTc                                                                                                                                        >  2:1648728/1‑150 (MQ=255)
                 ggATGATCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGa                                                                                                                                      <  2:1581947/150‑1 (MQ=255)
                  gatgatCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGAt                                                                                                                                     <  1:844119/150‑1 (MQ=255)
                  gatgatCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGAt                                                                                                                                     <  1:1603821/150‑1 (MQ=255)
                   atgatCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGCAACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATg                                                                                                                                    >  1:2273123/1‑150 (MQ=255)
                   atgatCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATg                                                                                                                                    <  2:756166/150‑1 (MQ=255)
                      atCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTcc                                                                                                                                 <  1:1643169/150‑1 (MQ=255)
                       tCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCt                                                                                                                                >  1:1539726/1‑150 (MQ=255)
                         gccgccAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGc                                                                                                                              >  1:807229/1‑150 (MQ=255)
                          ccgccAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCt                                                                                                                             >  1:2026187/1‑150 (MQ=255)
                            gccAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCGCGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCg                                                                                                                           >  2:572590/1‑150 (MQ=255)
                                    gtcgCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGc                                                                                                                   >  1:2178029/1‑150 (MQ=255)
                                     tcgCGTTCACGCAGGGCCGGCAAGCGCAGCCACACACCGTTGACCCGGTAGACCAGGTCCTAACGGCAGTGGTGCTCCTGGAGTACGGGCTTGAGGTCGCGGTGTGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGcc                                                                                                                  <  1:1306643/150‑1 (MQ=255)
                                      cgcgTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGTGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGccc                                                                                                                 <  1:171992/150‑1 (MQ=255)
                                      cgcgTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGccc                                                                                                                 <  1:681146/150‑1 (MQ=255)
                                           tCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAg                                                                                                            <  2:2178029/150‑1 (MQ=255)
                                            cACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAgg                                                                                                           <  1:49059/150‑1 (MQ=255)
                                              cGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGt                                                                                                         >  1:552480/1‑150 (MQ=255)
                                               gCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTg                                                                                                        <  2:74263/150‑1 (MQ=255)
                                                      cGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAggg                                                                                                          <  1:270468/141‑1 (MQ=255)
                                                      cGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAggg                                                                                                          >  2:270468/1‑141 (MQ=255)
                                                      cGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTAcggg                                                                                                 >  2:1092354/1‑148 (MQ=255)
                                                      cGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTAcgcg                                                                                                 <  1:1597703/150‑1 (MQ=255)
                                                        gCAGGCGCACCGACACCCCGTTGCCCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACgcggc                                                                                               <  2:535328/150‑1 (MQ=255)
                                                               cAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTAc                                                                                                    >  1:375495/1‑138 (MQ=255)
                                                               cAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTAc                                                                                                    <  2:375495/138‑1 (MQ=255)
                                                                 gCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCa                                                                                      <  2:2273123/150‑1 (MQ=255)
                                                                      acacCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCa                                                                                                                                                       >  1:1698481/1‑80 (MQ=255)
                                                                      acacCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCa                                                                                                                                                       <  2:1698481/80‑1 (MQ=255)
                                                                      acacCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGATGGCGCAGATCAGCGACACGTCGAGGTCCTGCTCGGCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTAcc                                                                                 >  2:2099417/1‑150 (MQ=255)
                                                                        acCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCg                                                                               >  1:1826205/1‑150 (MQ=255)
                                                                         ccGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGAt                                                                                                                                     >  2:1308705/1‑95 (MQ=255)
                                                                         ccGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGAt                                                                                                                                     <  1:1308705/95‑1 (MQ=255)
                                                                         ccGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGc                                                                              >  1:389848/1‑150 (MQ=255)
                                                                              gACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAg                                                                         <  2:676083/150‑1 (MQ=255)
                                                                                cccGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGc                                                                       <  2:389848/150‑1 (MQ=255)
                                                                                 ccGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCg                                                                      >  1:1317390/1‑150 (MQ=255)
                                                                                 ccGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCg                                                                      >  2:1335695/1‑150 (MQ=255)
                                                                                  cGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGt                                                                     <  2:552480/150‑1 (MQ=255)
                                                                                  cGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGt                                                                     <  1:1089154/150‑1 (MQ=255)
                                                                                    gTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGc                                                                   <  2:2386834/150‑1 (MQ=255)
                                                                                     tAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCt                                                                  <  2:667912/150‑1 (MQ=255)
                                                                                      aGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCtt                                                                 <  1:2096463/150‑1 (MQ=255)
                                                                                      aGAACAGGACTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCtt                                                                 <  2:1317390/150‑1 (MQ=255)
                                                                                       gAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCACTACCCGCAACAGGCGTGCTTg                                                                >  1:722649/1‑150 (MQ=255)
                                                                                        aaCAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTAcg                                                                                                   <  1:1287477/114‑1 (MQ=255)
                                                                                        aaCAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGAGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTAcg                                                                                                   >  2:1287477/1‑114 (MQ=255)
                                                                                         aCATGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGGCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCa                                                              <  1:312797/150‑1 (MQ=255)
                                                                                         aCAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCa                                                              >  1:1764920/1‑150 (MQ=255)
                                                                                         aCAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGCGAACGGCTGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCa                                                              <  1:24669/150‑1 (MQ=255)
                                                                                           aGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCTGATCTGCGCCa                                                                                                                                            >  2:975041/1‑70 (MQ=39)
                                                                                           aGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCTGATCTGCGCCa                                                                                                                                            <  1:975041/70‑1 (MQ=39)
                                                                                           aGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAgg                                                            >  2:313119/1‑150 (MQ=255)
                                                                                           aGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTAGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAgg                                                            >  2:139426/1‑150 (MQ=255)
                                                                                            ggTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAggg                                                           <  2:1539726/150‑1 (MQ=255)
                                                                                             gTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGc                                                                                                                   <  2:968467/93‑1 (MQ=255)
                                                                                             gTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGc                                                                                                                   >  1:968467/1‑93 (MQ=255)
                                                                                             gTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGc                                                          <  1:313119/150‑1 (MQ=255)
                                                                                             gTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGc                                                          >  1:866002/1‑150 (MQ=255)
                                                                                                ttCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGCGCTTGCAgcaccag                                                       >  2:1955507/1‑144 (MQ=255)
                                                                                                 tCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGc                                                      <  2:2269904/150‑1 (MQ=255)
                                                                                                   aCGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTTCACGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCgg                                                    >  1:145229/1‑150 (MQ=255)
                                                                                                   aCGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCgg                                                    >  1:1827686/1‑150 (MQ=255)
                                                                                                   aCGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCgg                                                    >  2:1139022/1‑150 (MQ=255)
                                                                                                    cGGAAGTGCTGCCCCGGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGGTCCGCGCCACGTCGCGGTCCTGCTCGTCACCCGCCCCCACGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGc                                                   <  2:153479/150‑1 (MQ=255)
                                                                                                         gTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTc                                              <  2:1764920/150‑1 (MQ=255)
                                                                                                           gctgctCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCTCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGc                                            >  2:806138/1‑150 (MQ=255)
                                                                                                            ctgctCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGcc                                           >  1:304789/1‑150 (MQ=255)
                                                                                                               ctCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTTGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGAt                                        <  2:1297682/150‑1 (MQ=255)
                                                                                                                 ccTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCCCGTCGATGTCCTGTGTGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGAttt                                      <  1:536484/150‑1 (MQ=255)
                                                                                                                      aCCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTtcgtcg                                 >  1:1689859/1‑150 (MQ=255)
                                                                                                                        cAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAg                               >  1:1429220/1‑150 (MQ=255)
                                                                                                                        cAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAg                               >  1:2042226/1‑150 (MQ=255)
                                                                                                                             gcTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTAGAGGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAGGAACa                         >  1:1136345/1‑150 (MQ=255)
                                                                                                                                 acactCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAGGAACAGCGt                     >  1:82681/5‑150 (MQ=255)
                                                                                                                                  aGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTGATTCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAGGAACAGCGTg                    >  1:413026/1‑150 (MQ=255)
                                                                                                                                   ggTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAGGAACAGCGTGc                   <  2:1826205/150‑1 (MQ=255)
                                                                                                                                       gcgGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTc                                              >  1:2132115/1‑120 (MQ=255)
                                                                                                                                       gcgGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTc                                              <  2:2132115/120‑1 (MQ=255)
                                                                                                                                          gtgggtggCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAGGAACAGCGTGCCACCATg            <  2:1827686/150‑1 (MQ=255)
                                                                                                                                           tgggtggCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAGGAACAGCGTGCCACCATgg           <  1:1335695/150‑1 (MQ=255)
                                                                                                                                             ggtggCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAGGAACAGCGTGCCACCATGGGc         <  1:276958/150‑1 (MQ=255)
                                                                                                                                              gtggCGCTGATCTGCGCCACGTCTATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCAGGTCGCCGATTTCGTCGAGGAACAGCGTTCCACCATGGGcc        >  2:1487676/1‑150 (MQ=255)
                                                                                                                                              gtggCGCTGATCTGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTAcgcg                                                                                                 >  2:790692/1‑62 (MQ=38)
                                                                                                                                              gtggCGCTGATCTGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTAcgcg                                                                                                 <  1:790692/62‑1 (MQ=38)
                                                                                                                                              gtggCGCTGATCTGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAgg                                                            >  1:2378983/1‑99 (MQ=255)
                                                                                                                                              gtggCGCTGATCTGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAgg                                                            <  2:2378983/99‑1 (MQ=255)
                                                                                                                                                 gcgcAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAGGAACAGCGTGCCACCATGGGCCTGc     >  2:2400329/1‑150 (MQ=255)
                                                                                                                                                   gcAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGAAGGGCGAGACGCATGCCGCTCATCTCCCGTGAGAACAGCGAACCGCCAGGGGCCTGCAt   >  2:623771/1‑150 (MQ=255)
                                                                                                                                                    cAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAGGAACAGCGTGCCACCATGGGCCTGCATc  <  1:572590/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTTGTCCAGCAGGCCCTGGATGATCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAGGAACAGCGTGCCACCATGGGCCTGCATC  >  NC_002947/3933882‑3934179

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: