Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 334,606 A→T 11.5% L28M (TTG→ATG) ‡ 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,6060AT11.5% 223.4 / 17.2 78L28M (TTG→ATG) ‡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 = 8.87e-01

GCATGTCCATCATGGCGTCCAAGGTGAATTTCTCGATCTTGCCATTTGCAAGCTCGGATATGCGGGATTGGGTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCTCGGCAGAAATATCTAAAATTTGCGATAAATCGTAGTGTGGATTGGTTGGCGAGTCAAA  >  NC_002947/334459‑334755
                                                                                                                                                    |                                                                                                                                                     
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                                                                                                                     <  2:887762/73‑1 (MQ=255)
                                                                                                             gCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCt                                                                                                                     >  1:887762/1‑73 (MQ=255)
                                                                                                              cccATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAggg                                                                                             >  1:76157/1‑96 (MQ=255)
                                                                                                              cccATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAggg                                                                                             <  2:76157/96‑1 (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                                                                                                <  2:2227357/68‑1 (MQ=37)
                                                                                                                                      cttgaAAATGGACTTGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGaaa                                                                                                >  1:2227357/3‑70 (MQ=37)
                                                                                                                                      cttgaAAATGGACTTGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCTCgg                                                         >  1:1442108/3‑109 (MQ=255)
                                                                                                                                      cttgaAAATGGACTTGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCTCgg                                                         <  2:1442108/107‑1 (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)
                                                                                                                                                    aGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCTCGGCAGAAATATCTAAAATTTGCGATAAATCGTAGTGTGGATTGGTTGGCGAGTCaaa  <  2:2012139/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCATGTCCATCATGGCGTCCAAGGTGAATTTCTCGATCTTGCCATTTGCAAGCTCGGATATGCGGGATTGGGTTACGTTCAGCTTGGCTGCGGCATCATGCTGTTTGAGGCCCATTTTTTTGATGCAGTCAGTGATGAAAATGGACAAGTCCATTTTCAAGGCCATTTTGCTGGCGGTATCTGCGTCGTGGGTGGCGTGGAAAGGGCTTTCGTGGAATTGAAGCGCCATCGGCTTGCCTCGGCAGAAATATCTAAAATTTGCGATAAATCGTAGTGTGGATTGGTTGGCGAGTCAAA  >  NC_002947/334459‑334755

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: