Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,934,036 A→T 8.5% L474Q (CTG→CAG)  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,0360AT 8.5% 286.8 / 14.8 95L474Q (CTG→CAG) PP_3467sigma‑54 dependent transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (40/46);  new base T (4/4);  total (44/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 = 7.16e-01

AGCAGGCCCTGGATGATCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAGGAACAGCGTGCCACCATGGGCCTGCATCAGC  >  NC_002947/3933889‑3934182
                                                                                                                                                   |                                                                                                                                                   
agcagGCCCTGGATGATCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGATTAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAgc                                                                                                                                                   <  1:319156/150‑1 (MQ=255)
agcagGCCCTGGATGATCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAgc                                                                                                                                                   <  2:1558195/150‑1 (MQ=255)
agcagGCCCTGGATGATCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAgc                                                                                                                                                   <  1:318226/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)
                                                               acacCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGATGGCGCAGATCAGCGACACGTCGAGGTCCTGCTCGGCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTAcc                                                                                    >  2:2099417/1‑150 (MQ=255)
                                                                 acCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCg                                                                                  >  1:1826205/1‑150 (MQ=255)
                                                                  ccGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGAt                                                                                                                                        <  1:1308705/95‑1 (MQ=255)
                                                                  ccGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGAt                                                                                                                                        >  2:1308705/1‑95 (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                                                                                                                      >  1:968467/1‑93 (MQ=255)
                                                                                      gTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGc                                                                                                                      <  2:968467/93‑1 (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                                                       >  2:1139022/1‑150 (MQ=255)
                                                                                            aCGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCgg                                                       >  1:1827686/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                                                                                                    <  1:790692/62‑1 (MQ=38)
                                                                                                                                       gtggCGCTGATCTGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTAcgcg                                                                                                    >  2:790692/1‑62 (MQ=38)
                                                                                                                                       gtggCGCTGATCTGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAgg                                                               <  2:2378983/99‑1 (MQ=255)
                                                                                                                                       gtggCGCTGATCTGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAgg                                                               >  1:2378983/1‑99 (MQ=255)
                                                                                                                                          gcgcAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAGGAACAGCGTGCCACCATGGGCCTGc        >  2:2400329/1‑150 (MQ=255)
                                                                                                                                            gcAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGAAGGGCGAGACGCATGCCGCTCATCTCCCGTGAGAACAGCGAACCGCCAGGGGCCTGCAt      >  2:623771/1‑150 (MQ=255)
                                                                                                                                             cAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAGGAACAGCGTGCCACCATGGGCCTGCATc     <  1:572590/150‑1 (MQ=255)
                                                                                                                                                aTCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAGGAACAGCGTGCCACCATGGGCCTGCATCAgc  <  2:1781817/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AGCAGGCCCTGGATGATCGCCGCCAGGTCGTCGCGTTCACGCAGGGCCGGCAGGCGCAGCGACACACCGTTGACCCGGTAGAACAGGTCTTCACGGAAGTGCTGCTCCTGCACCAGGCGCTTGAGGTCGCGGTGGGTGGCGCAGATCAGCGCCACGTCGATGTCCTGCTCGTCACCCGCCCCCAAGGGTGCTACGCGGCGCTCCTGCAGTACCCGCAACAGGCGTGCTTGCAGGGCGAGCGGCATGTCGCCGATTTCGTCGAGGAACAGCGTGCCACCATGGGCCTGCATCAGC  >  NC_002947/3933889‑3934182

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: