Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,866,854 A→C 53.7% V23G (GTT→GGT)  PP_3415 ← LacI family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,866,8540AC53.7% 25.6 / 78.8 67V23G (GTT→GGT) PP_3415LacI family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (17/14);  new base C (15/21);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 3.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

CGACTGGCTGCGTGCCGAGGCCAGCGCCCGGGCGGCCGGGTTGGCCACATAGCCAAGGTTCGCCGCTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAAT  >  NC_002947/3866709‑3867000
                                                                                                                                                 |                                                                                                                                                  
cGACTGGCTGCGTGCCGAGGCCAGCGCCCGGGCGGCCGGGTTGGCCACATAGCCAAGGTTCGCCGCTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAACcccc                                                                                                                                                >  2:1018468/1‑150 (MQ=255)
     ggCTGCGTGCCGAGGCCAGCGCCCGGGCGGCCGGGTTGGCCACATAGCCAAGGTTCGCCGCTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGaaa                                                                                                                                           >  1:706721/1‑150 (MQ=255)
     ggCTGCGTGCCGAGGCCAGCGCCCGGGCGGCAGGGTTGGCCACATAGCCAAGGTTCGCCGCTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAACCCCCGGaaa                                                                                                                                           >  1:273062/1‑150 (MQ=255)
      gCTGCGTGCCGAGGCCAGCGCCCGGGCGGCCGGGTTGGCCACATAGCCAAGGTTCGCCGCTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCCGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAg                                                                                                                                          >  2:1875627/1‑150 (MQ=255)
         gCGTGCCGAGGCCAGCGCCCGGGCGGCCGGGTTGGCCACATAGCCAAGGTTCGCCGCTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCt                                                                                                                                       <  2:1109445/150‑1 (MQ=255)
          cGTGCCGAGGCCAGCGCCCGGGCGGCCGGGTTGGCCACATAGCCAAGGTTCGCCGCTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAACCCCCGGAAAGCCTg                                                                                                                                      <  1:1895154/150‑1 (MQ=255)
                     cAGCGCCCGGGCGGCCGGGTTGGCCACATAGCCAAGGTTCGCCGCTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAACCCCCGGAAAGCCTGGCAACTTCgag                                                                                                                           >  1:418027/1‑148 (MQ=255)
                     cAGCGCCCGGGCGGCCGGGTTGGCCACATAGCCAAGGTTCGCCGCTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAACCCCCGGAAAGCCTGGCAACTTCGGc                                                                                                                           >  1:209702/1‑150 (MQ=255)
                     cAGCGCCCGGGCGGCCGGGTTGGCCACATAGCCAAGGTTCGCCGCTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGc                                                                                                                           <  2:1560711/150‑1 (MQ=255)
                     cAGCGCCCGGGCGGCCGGGTTGGCCACATAGCCAAGGTTCGCCGCTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGc                                                                                                                           >  1:1560711/1‑150 (MQ=255)
                     cAGCGCCCGGGCGGCCGGGTTGGCCACATAGCCAAGGTTCGCCGCTGCGGCGATGACCTTTTCCACCACGTTCTGTTCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAACCCCCGGAAAGCCTGGCAACTTCGGc                                                                                                                           <  1:723915/150‑1 (MQ=255)
                      aGCGCCCGGGCGGCCGGGTTGGCCACATAGCCAAGGTTCGCCGCTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAACCCCCGGAAAGCCTGGCAACTTCGGcc                                                                                                                          <  2:209702/150‑1 (MQ=255)
                         tcccGGGCGGCCGGGTTTGCCACATGGCCAGGGTTCGCCGCTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGTGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAgg                                                                                                                       <  2:141810/149‑1 (MQ=255)
                           ccGGGCGGCCGGGTTGGCCACATAGCCAAGGTTCGCCGCTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAACCCCCGGAAAGCCTGGCAACTTCGGCCAGGGt                                                                                                                     <  2:439818/150‑1 (MQ=255)
                            cGGGCGGCCGGGTTGGCCACATAGCCAAGGTTCGCCGCTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTg                                                                                                                    >  2:1406744/1‑150 (MQ=255)
                             gggCGGCCGGGTTGGCCACATAGCCAAGGTTCGCCGCTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAACCCCCg                                                                                                                                               <  2:930329/122‑1 (MQ=255)
                             gggCGGCCGGGTTGGCCACATAGCCAAGGTTCGCCGCTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAACCCCCg                                                                                                                                               >  1:930329/1‑122 (MQ=255)
                              ggcggcCGGGTTGGCCACATAGCCAAGGTTCGCCGCTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAACCCCCGGAAAGCCTGGCAACTTCGGCCAGGGTggg                                                                                                                  <  1:589529/150‑1 (MQ=255)
                                      ggTTGGCCACATAGCCAAGGTTCGCCGCTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAACCCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTg                                                                                                          >  1:498396/1‑150 (MQ=255)
                                       gTTGGCCACATAGCCAAGGTTCGCCGCTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAACCCCCGGAAAGCCt                                                                                                                                       <  2:631451/120‑1 (MQ=255)
                                       gTTGGCCACATAGCCAAGGTTCGCCGCTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAACCCCCGGAAAGCCt                                                                                                                                       >  1:631451/1‑120 (MQ=255)
                                         tGGCCACATAGCCAAGGTTCGCCGCTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAACCCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAg                                                                                                       <  1:1291607/150‑1 (MQ=255)
                                             cacaTAGCCAAGGTTCGCCGCTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAAAcc                                                                                                                                                  <  1:217292/103‑1 (MQ=255)
                                             cacaTAGCCAAGGTTCGCCGCTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAAAcc                                                                                                                                                  >  2:217292/1‑103 (MQ=255)
                                             cacaTAGCCAAGGTTCGCCGCTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCg                                                                                                   <  1:1591717/150‑1 (MQ=255)
                                              acaTAGCCAAGGTTCGCCGCTGCGGAGATTACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCgg                                                                                                  >  1:1509194/1‑150 (MQ=255)
                                                aTAGCCAAGGTTCGCCGCTGCGGCGATGACCTTTTCCGCCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTgg                                                                                                                                     <  2:833944/113‑1 (MQ=255)
                                                aTAGCCAAGGTTCGCCGCTGCGGCGATGACCTTTTCCGCCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTgg                                                                                                                                     >  1:833944/1‑113 (MQ=255)
                                                aTAGCCAAGGTTCGCCGCTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTgg                                                                                                                                     >  1:834116/1‑113 (MQ=255)
                                                aTAGCCAAGGTTCGCCGCTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTgg                                                                                                                                     <  2:834116/113‑1 (MQ=255)
                                                    ccAAGGTTCGCCGCTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCg                                                                                            >  2:1802467/1‑150 (MQ=255)
                                                    ccAAGGTTCGCCGCTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCg                                                                                            >  2:690297/1‑150 (MQ=255)
                                                        ggTTCGCCGCTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAACCCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGc                                                                                        <  1:1526465/150‑1 (MQ=255)
                                                             gccgcTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCAGTGCGGGAc                                                                                   >  1:1052971/1‑150 (MQ=255)
                                                              ccgcTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAACCCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACa                                                                                  <  1:1018468/150‑1 (MQ=255)
                                                               cgcTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAACCCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACAt                                                                                 <  2:1066667/150‑1 (MQ=255)
                                                                    cggcgATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTg                                                                                         <  2:589519/137‑1 (MQ=255)
                                                                    cggcgATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTg                                                                                         >  1:589519/1‑137 (MQ=255)
                                                                          tGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAACCCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTaa                                                                      >  1:1895078/1‑150 (MQ=255)
                                                                           gACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAACCCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAAt                                                                     >  2:64227/1‑150 (MQ=255)
                                                                                     accaGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACtt                                                           >  1:139379/1‑150 (MQ=255)
                                                                                         ggTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCaa                                                       <  2:706721/150‑1 (MQ=255)
                                                                                          gTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGGAATTTTACTACTTGCAAg                                                      >  1:1068829/1‑150 (MQ=255)
                                                                                              ggCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAACCCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCgg                                                                                      <  1:1793487/114‑1 (MQ=255)
                                                                                              ggCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAACCCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCgg                                                                                      >  2:1793487/1‑114 (MQ=255)
                                                                                               gcgcAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAACCCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGt                                                                                          <  2:1589549/109‑1 (MQ=255)
                                                                                               gcgcAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAACCCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGt                                                                                          >  1:1589549/1‑109 (MQ=255)
                                                                                               gcgcAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAACCCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGagag                                                 >  1:958158/1‑150 (MQ=255)
                                                                                                cgcAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAACCCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAgg                                                <  1:552610/150‑1 (MQ=255)
                                                                                                             aCCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTaa                                   <  1:386658/150‑1 (MQ=255)
                                                                                                               cccACGCAAGGCGCGCGAAGCGGTAATCGGGGAACCCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAAtg                                                                            <  2:722238/107‑1 (MQ=255)
                                                                                                               cccACGCAAGGCGCGCGAAGCGGTAATCGGGGAACCCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAAtg                                                                            >  1:722238/1‑107 (MQ=255)
                                                                                                                        ggCGCGCGAAGCGGTAATCGGGGAACCCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGAGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGt                        <  2:273062/150‑1 (MQ=255)
                                                                                                                             gcgAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCTGGTGCGGGACATGAATGTTGTAATTTTACTACTATCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAg                   >  1:990087/1‑150 (MQ=255)
                                                                                                                               gAAGCGGTAATCGGGGAACCCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAAc                                            <  2:1815625/123‑1 (MQ=255)
                                                                                                                               gAAGCGGTAATCGGGGAACCCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAAc                                            >  1:1815625/1‑123 (MQ=255)
                                                                                                                                    ggTAATCGGGGAAAACCCGGACAGCCTGGCACATTCTGCCATGGTGGGACGACCTGTAGTTCGGGAGCTGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAg            <  2:990087/150‑1 (MQ=255)
                                                                                                                                     gTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGcc                                                                                             >  2:631214/1‑68 (MQ=255)
                                                                                                                                     gTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGcc                                                                                             <  1:631214/68‑1 (MQ=255)
                                                                                                                                      tAATCGGGGAACCCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGc          >  1:222396/1‑150 (MQ=255)
                                                                                                                                       aaTCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCa         >  1:609050/1‑150 (MQ=255)
                                                                                                                                        aTCGGGGAACCCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCaa        <  1:1760335/150‑1 (MQ=255)
                                                                                                                                        aTCGGGGAACCCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCaa        <  2:108275/150‑1 (MQ=255)
                                                                                                                                         tCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAAt       <  1:690297/150‑1 (MQ=255)
                                                                                                                                          cGGGGAACCCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATg      <  2:1370635/150‑1 (MQ=255)
                                                                                                                                           ggggAACCCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGt     <  1:177898/150‑1 (MQ=255)
                                                                                                                                              gAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAAt  <  2:1813476/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CGACTGGCTGCGTGCCGAGGCCAGCGCCCGGGCGGCCGGGTTGGCCACATAGCCAAGGTTCGCCGCTGCGGCGATGACCTTTTCCACCAGGTCAGGCGCAACGGTGCTGACCCCACGCAAGGCGCGCGAAGCGGTAATCGGGGAAACCCCGGAAAGCCTGGCAACTTCGGCCAGGGTGGGACGACCTGTAGTTCGGGAGCCGGTGCGGGACATGAATGTTGTAATTTTACTACTTGCAAGGAGAGGGAACGGCCACTAAGGTAGCGCTGTCTCAGGATGCAGGCAATGTAAT  >  NC_002947/3866709‑3867000

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: