Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 334,605 A→T 11.3% L28* (TTG→TAG) ‡ PP_0276 ← Cro/CI family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947334,6050AT11.3% 219.0 / 13.9 78L28* (TTG→TAG) ‡PP_0276Cro/CI family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (38/31);  new base T (5/4);  total (43/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

CCAGCATGTCCATCATGGCGTCCAAGGTGAATTTCTCGATCTTGCCATTTGCAAGCTCGGATATGCGGGATTGGGTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCTCGGCAGAAATATCTAAAATTTGCGATAAATCGTAGTGTGGATTGGTTGGCGAGTCAA  >  NC_002947/334456‑334754
                                                                                                                                                      |                                                                                                                                                     
ccAGCATGTCCATCATGGCGTCCAAGGTGAATTTCTCGATCTTGCCATTTGCAAGCTCGGATATGCGGGATTGGGTTACGTTCAGCTTGGCTGCGGCCTCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACa                                                                                                                                                       <  1:1543467/150‑1 (MQ=255)
   gCATGTCCATCATGGCGTCCAAGGTGAATTTCTCGATCTTGCCATTTGCAAGCTCGGATATGCGGGATTGGGTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGt                                                                                                                                                    >  2:1499520/1‑150 (MQ=255)
   gCATGTCCATCATGGCGTCCAAGGTGAATTTCTCGATCTTGCCATTTGCAAGCTCGGATATGCGGGATTGGGTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGAAAAGt                                                                                                                                                    >  1:447322/1‑150 (MQ=255)
       gTCCATCATGGCGTCCAAGGTGAATTTCTCGATCTTGCCATTTGCAAGCTCGGATATGCGGGATTGGGTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCAt                                                                                                                                                >  1:1912175/1‑150 (MQ=255)
        tCCATCATGGCGTCCAAGGTGAATTTCTCGATCTTGCCATTTGCAAGCTCGGATATGCGGGATTGGGTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGATGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCAtt                                                                                                                                               <  1:1499520/150‑1 (MQ=255)
           atcatGGCGTCCAAGGTGAATTTCTCGATCTTGCCATTTGCAAGCTCGGATATGGGGGATTGGGTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTc                                                                                                                                            >  1:227920/1‑150 (MQ=255)
             catGGCGTCCAAGGTGAATTTCTCGATCTTGCCATTTGCAAGCTCGGATATGCGGGATTGGGTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCaa                                                                                                                                          <  1:1388825/150‑1 (MQ=255)
              atGGCGTCCAAGGTGAATTTCTCGATCTTGCCATTTGCAAGCTCGGATATGCGGGATTGGGTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAg                                                                                                                                         <  1:1462249/150‑1 (MQ=255)
               tGGCGTCCAAGGTGAATTTCTCGATCTTGCCATTTGCAAGCTCGGATATGCGGGATTGGGTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCACATGGTAAAATGGACAAGTCCATTTTCAAgg                                                                                                                                        <  1:1643921/150‑1 (MQ=255)
                        aGGTGAATTTCTCGATCTTGCCATTTGCAAGCTCGGATATGCGGGATTGGGTTACGTTCAGCTTGTCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGc                                                                                                                               <  2:42990/150‑1 (MQ=255)
                         ggTGAATTTCTCGATCTTGCCATTTGCAAGCTCGGATATGCGGGATTGGGTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCt                                                                                                                              >  2:2127059/1‑150 (MQ=255)
                           tGAATTTCTCGATCTTGCCATTTGCAAGCTCGGATATGCGGGATTGGGTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACTTGTCCAttt                                                                                                                                              >  2:441001/1‑132 (MQ=255)
                           tGAATTTCTCGATCTTGCCATTTGCAAGCTCGGATATGCGGGATTGGGTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACTTGTCCAttt                                                                                                                                              <  1:441001/132‑1 (MQ=255)
                            gAATTTCTCGATCTTGCCATTTGCAAGCTCGGATATGCGGGATTGGGTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGc                                                                                                                           >  1:668785/1‑150 (MQ=255)
                                  ctcGATCTTGCCATTTGCAAGCTCGGATATGCGGGATTGGGTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGTCATTTTTCTGGCGGTAtt                                                                                                                     >  2:972126/1‑149 (MQ=255)
                                  ctcGATCTTGCCATTTGCAAGCTCGGATATGCGGGATTGGGTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATc                                                                                                                     >  2:142101/1‑150 (MQ=255)
                                  ctcGATCTTGCCATTTGCAAGCTCGGATATGCGGGATTGGGTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATc                                                                                                                     >  2:1532864/1‑150 (MQ=255)
                                             cATTTGCAAGCTCGGATAGGCGGGATTGGGTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCgtggg                                                                                                          >  1:2012139/1‑150 (MQ=255)
                                              aTTTGCAAGCTCGGATATGCGGGATTGGGTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCgtgggt                                                                                                         >  1:870060/1‑150 (MQ=255)
                                               tttGCAAGCTCGGATATGCGGGATTGGGTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCgtgggtg                                                                                                        <  1:1532864/150‑1 (MQ=255)
                                                   cAAGCTCGGATATGCGGGATTGGGTTACGTTCAGCTTGGCTGCGGCTTCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTTCTGGCGGTATCTGCGTCGTGGGTGGCGt                                                                                                    >  2:1825339/1‑150 (MQ=255)
                                                   cAAGCTCGGATATGCGGGATTGGGTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGt                                                                                                    >  1:1648451/1‑150 (MQ=255)
                                                    aaGCTCGGATATGCGGGATTGGGTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACTTGTCCAttt                                                                                                                                              >  1:1624465/1‑107 (MQ=255)
                                                    aaGCTCGGATATGCGGGATTGGGTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACTTGTCCAttt                                                                                                                                              <  2:1624465/107‑1 (MQ=255)
                                                    aaGCTCGGATATGAGGGATTGGGTTACGTTCAGTTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGTAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTg                                                                                                   <  1:1148722/150‑1 (MQ=255)
                                                      gCTCGGATATGCGGGATTGGGTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGa                                                                                                 >  1:999858/1‑150 (MQ=255)
                                                        tCGGATATGCGGGATTGGGTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTTGTGCGTGGCGTGGaaa                                                                                               <  2:227920/150‑1 (MQ=255)
                                                        tCGGATATGCGGGATTGGGTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTAAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGaaa                                                                                               >  1:1717982/1‑150 (MQ=255)
                                                         cGGATATGCGGGATTGGGTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAg                                                                                              <  2:1912175/150‑1 (MQ=255)
                                                          ggATATGCGGGATTGGGTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAga                                                                                             <  2:1648451/150‑2 (MQ=255)
                                                                  gggATTGGGTTACGTTCAGCTTGGCTGCGGCTTCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTTCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGt                                                                                     <  1:1825339/150‑1 (MQ=255)
                                                                          gTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCTTCGTGGTTTGCGTGGAAAGGGCTTTCGTGGTATTGa                                                                             >  2:339793/1‑150 (MQ=255)
                                                                          gTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGaaa                                                                                               <  1:267085/131‑1 (MQ=255)
                                                                          gTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGaaa                                                                                               >  2:267085/1‑131 (MQ=255)
                                                                           ttACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGaa                                                                            >  1:2250861/1‑150 (MQ=255)
                                                                             aCGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAgc                                                                          <  2:2250861/150‑1 (MQ=255)
                                                                                 tCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCa                                                                      <  1:2127059/150‑1 (MQ=255)
                                                                                    gCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCt                                                                                                                    <  2:314708/101‑1 (MQ=255)
                                                                                    gCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCt                                                                                                                    >  1:314708/1‑101 (MQ=255)
                                                                                    gCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAgc                                                                          <  1:1178883/143‑1 (MQ=255)
                                                                                    gCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAgc                                                                          >  2:1178883/1‑143 (MQ=255)
                                                                                     cTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACaa                                                                                                                                                      <  2:1063128/66‑1 (MQ=255)
                                                                                     cTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACaa                                                                                                                                                      >  1:1063128/1‑66 (MQ=255)
                                                                                     cTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCgg                                                                  <  2:668785/150‑1 (MQ=255)
                                                                                      ttGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGc                                                                 >  1:5156/1‑150 (MQ=255)
                                                                                      ttGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGAAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGc                                                                 >  2:33059/1‑150 (MQ=255)
                                                                                                     tGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCTCGGCAGaaa                                                  >  2:2127516/1‑150 (MQ=255)
                                                                                                      gCTGTTTGAGGCCCATTTTTTTGATGCATTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTTTCGGTATCTGCGTCGTGGGTGGCTTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCTGCTTGCCTCGGCAGaact                                                 >  2:180546/1‑148 (MQ=255)
                                                                                                                gCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCt                                                                                                                    >  1:887762/1‑73 (MQ=255)
                                                                                                                gCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCt                                                                                                                    <  2:887762/73‑1 (MQ=255)
                                                                                                                 cccATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAggg                                                                                            <  2:76157/96‑1 (MQ=255)
                                                                                                                 cccATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAggg                                                                                            >  1:76157/1‑96 (MQ=255)
                                                                                                                  ccATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCTCGGCAGAAATATCTAGAATTTg                                     >  1:1176728/1‑150 (MQ=255)
                                                                                                                  ccATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCTCGGCAGAAATATCTAAAATTTg                                     <  2:447322/150‑1 (MQ=255)
                                                                                                                   cATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCTCGGCAGAAATATATAAAATTTGc                                    >  1:115433/1‑150 (MQ=255)
                                                                                                                     tgttgttGATGCAGTCAGTGATGAAAATGGAAAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCTCGGCAGAAATATCTAAAATTTGCGa                                  <  2:395662/145‑1 (MQ=255)
                                                                                                                          ttGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTGTCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATGGAAGCGCCATCGGGTTGCCGCGGCAGAAATAGCTGAAATTTGCGATAAAt                             >  2:177570/1‑150 (MQ=255)
                                                                                                                            gATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGc                                                                                                                  <  1:2308064/63‑1 (MQ=255)
                                                                                                                            gATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCTCGGCAGAAATATCTAAAATTTGCGATAAATCg                           >  1:1294611/1‑150 (MQ=255)
                                                                                                                            gATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGAGGTATATGc                                                                                                                  >  2:2308064/1‑63 (MQ=255)
                                                                                                                                   tcagtGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGGGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCAt                                                                     >  1:1499892/1‑101 (MQ=255)
                                                                                                                                   tcagtGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCAt                                                                     <  2:1499892/101‑1 (MQ=255)
                                                                                                                                       gtGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCTCGGCAGAAATATCTAAAATTTGCGATAAATCGTAGTGTGGAt                 <  1:33059/150‑1 (MQ=255)
                                                                                                                                        ggatgaAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCTCGGCAGAAATATCTAAAATTTGCGATAAATCGTAGTGTGGAtt                >  1:117631/2‑150 (MQ=255)
                                                                                                                                         cttgaAAATGGACTTGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGaaa                                                                                               >  1:2227357/3‑70 (MQ=37)
                                                                                                                                         cttgaAAATGGACTTGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGaaa                                                                                               <  2:2227357/68‑1 (MQ=37)
                                                                                                                                         cttgaAAATGGACTTGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCTCgg                                                        <  2:1442108/107‑1 (MQ=255)
                                                                                                                                         cttgaAAATGGACTTGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCTCgg                                                        >  1:1442108/3‑109 (MQ=255)
                                                                                                                                         cttgaAAATGGACAAGTCCATTTTCAAGGCCGTTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCTCGGCAGAAATATCTAAAATTTGCGa                                  >  2:2211763/3‑131 (MQ=255)
                                                                                                                                         cttgaAAATGGACAAGTCCATTTTCAAGGCCGTTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCTCGGCAGAAATATCTAAAATTTGCGa                                  <  1:2211763/129‑1 (MQ=255)
                                                                                                                                          atgaAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCTCGGCAGAAATATCTAAAATTTGCGATAAATCGTAGTGTGGAttgg              <  1:142101/150‑1 (MQ=255)
                                                                                                                                           tgaAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCTCGGCAGAAATATCTAAAATTTGCGATAAATCGTAGTGTGGAttggt             >  2:394620/1‑150 (MQ=255)
                                                                                                                                           tgaAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCTCGGCAGAAATATCTAAAATTTGCGATAAATCGTAGTGTGGAttg               >  1:1629590/1‑148 (MQ=255)
                                                                                                                                           tgaAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCTCGGCAGAAATATCTAAAATTTGCGATAAATCGTAGTGTGGAttg               <  2:1629590/148‑1 (MQ=255)
                                                                                                                                             aaaaTGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCTCGGCAGAAATATCTAAAATTTGCGATAAATCGTAGTGTGGAttggttg           <  2:117631/150‑1 (MQ=255)
                                                                                                                                              aaaTGGACTTGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCTCGGCAGAAATATCTAAAATTTGCGATAAATCGTAGTGTGGAttggttgg          >  1:116992/1‑150 (MQ=255)
                                                                                                                                              aaaTGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCTCGGCAGAAATATCTAAAATTTGCGATAAATCGTAGTGTGGAttggttgg          >  2:843148/1‑150 (MQ=255)
                                                                                                                                                 tgtaCAAGTCCATTTTCAAGGCCATTTTGTTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCTCGGCAGAAATATCTAAAATTTGCGATAAATCGTAGTGTGGATTGGTTGGCGa       <  2:999858/147‑1 (MQ=255)
                                                                                                                                                 tGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGAAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCTCGGCAGAAATATCTAAAATTTGCGATAAATCGTAGTGTGGATTGGTTGGCGa       <  1:972126/150‑1 (MQ=255)
                                                                                                                                                   gACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCTCGGCAGAAATATCTAAAATTTGCGATAAATCGTAGTGTGGATTGGTTGGCGAGt     <  2:1294611/150‑1 (MQ=255)
                                                                                                                                                    aCAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCTCGGCAGAAATATCTAAAATTTGCGATAAATCGTAGTGTGGATTGGTTGGCGAGTc    >  1:1057326/1‑150 (MQ=255)
                                                                                                                                                      aaGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCTCGGCAGAAATATCTAAAATTTGCGATAAATCGTAGTGTGGATTGGTTGGCGAGTCaa  >  2:46192/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CCAGCATGTCCATCATGGCGTCCAAGGTGAATTTCTCGATCTTGCCATTTGCAAGCTCGGATATGCGGGATTGGGTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCTCGGCAGAAATATCTAAAATTTGCGATAAATCGTAGTGTGGATTGGTTGGCGAGTCAA  >  NC_002947/334456‑334754

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