Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,693,728 C→G 7.3% P52R (CCT→CGT) ‡ PP_1489 → CheW domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,693,7280CG 7.3% 260.0 / 8.7 82P52R (CCT→CGT) ‡PP_1489CheW domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (39/37);  new base G (3/3);  total (42/40)
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.83e-01

TTGCAACTGCGCACCGCTCGGGCGGCGAACACCGAGGGGGCGTTGTACCTGCAGTTTCGCATCAGAGCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTG  >  NC_002947/1693583‑1693876
                                                                                                                                                  |                                                                                                                                                    
ttGCAACTGCGCACCGCTCGGGCGGCGAACACCGAGGGGGCGTTGTACCTGCAGTTTCGCATCAGAGCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCAcctgcc                                                                                                                                                  >  1:1008266/1‑150 (MQ=255)
    aaCTGCGCACCGCTCGGGCGGCGAACACCGAGGGGGCGTTGTACCTGCAGTTTCGCATCAGAGCGCAGCGTTTTGCCCTGAGCTTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTggg                                                                                                                                              >  1:668509/1‑150 (MQ=255)
    aaCTGCGCACCGCTCGGGCGGCGAACACCGAGGGGGCGTTGTACCTGCAGTTTCGCATCAGAGCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCAGGTGCCTggg                                                                                                                                              <  2:1202236/150‑1 (MQ=255)
     aCTGCGCACCGCTCGGGCGGCGAACACCGAGGGGGCGTTGTACCTGCAGTTTCGCATCAGAGCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCACCCAGGCAGGTGCCTgggc                                                                                                                                             <  2:1192083/150‑2 (MQ=255)
         cgcACCGCTCGGGCGGCGAACACCGAGGGGGCGTTGTACCTGCAGTTTCGCATCAGAGCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTcgcg                                                                                                                                         <  1:913243/150‑1 (MQ=255)
         cgcACCGCTCGGGCGGCGAACACCGAGGGGGCGTTGTACTTCAGATGTGTGATCATAGCGTATCGTTTTGCCCTTAGCTGGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCACCCAGGCAGGTGCCTgggcgatc                                                                                                                                         <  1:319204/150‑6 (MQ=255)
             ccGCTCGGGCGGCGAACACCGAGGGGGCGTTGTACCTGCAGTTTCGCATCAGAGCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCa                                                                                                                                     >  1:694099/1‑150 (MQ=255)
                    ggcggcGGACACCGAGGGGGCGTTGTAAATGAAGTTTCGCATCAGCGCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGTTCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGc                                                                                                                              <  2:1008266/150‑1 (MQ=255)
                         cgAACACCGAGGGGGCGTTGTACCTGCAGTTTCGCATCAGAGCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATc                                                                                                                         <  1:585594/150‑1 (MQ=255)
                          gAACACCGAGGGGGCGTTGTACCTGCAGTTTCGCATCAGAGCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCAGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCg                                                                                                                        >  2:665667/1‑150 (MQ=255)
                                    ggggCGTTGTACCTGCAGTTTCGCATCAGAGCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTg                                                                                                              >  1:1087145/1‑150 (MQ=255)
                                       gCGTTGTACCTGCAGTTTCGCATCAGAGCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGtt                                                                                                           >  1:477224/1‑150 (MQ=255)
                                        cGTTGTACCTGCAGTTTCGCATCAGAGCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCGGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTc                                                                                                          <  1:564999/150‑1 (MQ=255)
                                         gTTGTACCTGCAGTTTCGCATCAGAGCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTcc                                                                                                         >  2:1142597/1‑150 (MQ=255)
                                          ttGTACCTGCAGTTTCGCATCAGAGCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCg                                                                                                        >  2:782496/1‑150 (MQ=255)
                                               ccTGCAGTTTCGCATCAGAGCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCGGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCAt                                                                                                   >  1:948847/1‑150 (MQ=255)
                                                cTGCAGTTTCGCATCAGAGCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATc                                                                                                  <  2:263995/150‑1 (MQ=255)
                                                cTGCAGTTTCGCATCAGAGCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATc                                                                                                  <  2:477224/150‑1 (MQ=255)
                                                  gCAGTTTCGCATCAGAGCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCt                                                                                                               >  1:1364458/1‑135 (MQ=255)
                                                  gCAGTTTCGCATCAGAGCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCt                                                                                                               <  2:1364458/135‑1 (MQ=255)
                                                  gCAGTTTCGCATCAGAGCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCg                                                                                                        <  1:1262799/142‑1 (MQ=255)
                                                  gCAGTTTCGCATCAGAGCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCg                                                                                                        >  2:1262799/1‑142 (MQ=255)
                                                   cAGTTTCGCATCAGAGCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGAc                                                                                               <  2:934883/150‑1 (MQ=255)
                                                          cGCATCAGAGCGCAGCGTTTCGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCTACCTTgcc                                                                                         >  1:1318824/1‑150 (MQ=255)
                                                              tCAGAGCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCt                                                                                     <  2:301058/150‑1 (MQ=255)
                                                               cAGAGCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCtt                                                                                    <  2:668509/150‑1 (MQ=255)
                                                                agagCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCt                                                                                             <  1:1172787/140‑1 (MQ=255)
                                                                agagCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCt                                                                                             >  2:1172787/1‑140 (MQ=255)
                                                                 gagCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGa                                                                                  >  1:1395931/1‑150 (MQ=255)
                                                                  agCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAg                                                                                 >  1:514455/1‑150 (MQ=255)
                                                                  agCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAg                                                                                 <  1:1290779/150‑1 (MQ=255)
                                                                   gcgcAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGATCTTGCCGCCTTGAGc                                                                                <  1:759219/150‑1 (MQ=255)
                                                                       aGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCgg                                                                                                                                        >  1:1394992/1‑90 (MQ=255)
                                                                       aGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCgg                                                                                                                                        <  2:1394992/90‑1 (MQ=255)
                                                                       aGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCgg                                                                                                                                        >  2:421779/1‑90 (MQ=255)
                                                                       aGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCgg                                                                                                                                        <  1:421779/90‑1 (MQ=255)
                                                                       aGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTg                                                                            <  2:506090/150‑1 (MQ=255)
                                                                           ttttGCCCTGCGCGCGCACGAGGTGATCGAAGTGTTGCCGCGTCGGCCGCTCAAGCCGGTCGCCCAGGCACCTGCCCGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCa                                                                        <  1:810701/150‑1 (MQ=255)
                                                                                  cTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCcggcg                                                                 >  1:301128/1‑150 (MQ=255)
                                                                                   tGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCcg                                                                    >  1:853213/1‑146 (MQ=255)
                                                                                   tGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCcg                                                                    <  2:853213/146‑1 (MQ=255)
                                                                                      gCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGt                                                                                                                       <  2:456314/92‑1 (MQ=255)
                                                                                      gCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGt                                                                                                                       >  1:456314/1‑92 (MQ=255)
                                                                                               aGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTAGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCa                                                    <  2:408703/150‑1 (MQ=255)
                                                                                                  tGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCa                                                 <  1:1312362/150‑1 (MQ=255)
                                                                                                        aaGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGAc                                                                                               <  2:968894/98‑1 (MQ=255)
                                                                                                        aaGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGcc                                           >  1:709523/1‑150 (MQ=255)
                                                                                                        aaGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGACTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGAc                                                                                               >  1:968894/1‑98 (MQ=255)
                                                                                                             ttGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCctggtg                                      >  1:174386/1‑150 (MQ=255)
                                                                                                                ccGCGTCAGCCGCTCAAGCCGTTCGCCCAGGTATCTGCCTGGGTCTCGGGCATTCTCGCCTATCGTGTCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCAGGcggttgctg                                   <  2:98520/150‑1 (MQ=255)
                                                                                                                ccGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCctggtgctg                                   <  2:785291/150‑1 (MQ=255)
                                                                                                                ccGCGTCAGCCGCTCAAGACGATCGCCCAGGCATCGTGCTTGTTCGCGTGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCctggtgctg                                   <  2:514455/150‑1 (MQ=255)
                                                                                                                  gcgTCGGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCctggtgctggt                                 <  2:948847/150‑1 (MQ=255)
                                                                                                                     tCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCa                              >  1:48192/1‑150 (MQ=255)
                                                                                                                        gCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTa                           >  2:1248289/1‑150 (MQ=255)
                                                                                                                        gCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTa                           >  1:9907/1‑150 (MQ=255)
                                                                                                                         ccGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTAc                          >  1:847847/1‑150 (MQ=255)
                                                                                                                         ccGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTAc                          >  1:442681/1‑150 (MQ=255)
                                                                                                                           gCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCg                        >  1:256647/1‑150 (MQ=255)
                                                                                                                             tCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCTACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAg                      <  2:1318824/150‑1 (MQ=255)
                                                                                                                             gcAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAg                      <  1:911832/149‑1 (MQ=255)
                                                                                                                               aaGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCAc                                                      <  1:1316494/116‑1 (MQ=255)
                                                                                                                               aaGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCAc                                                      >  2:1316494/1‑116 (MQ=255)
                                                                                                                               aaGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAggc                    >  1:256437/1‑150 (MQ=255)
                                                                                                                               aaGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAggc                    <  2:577211/150‑1 (MQ=255)
                                                                                                                               aaGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAggc                    >  1:412642/1‑150 (MQ=255)
                                                                                                                                 gCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCtt                                                                                    <  1:66009/84‑1 (MQ=255)
                                                                                                                                 gCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCtt                                                                                    >  2:66009/1‑84 (MQ=255)
                                                                                                                                 gCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCAc                                                <  1:1239378/120‑1 (MQ=255)
                                                                                                                                 gCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCAc                                                >  2:1239378/1‑120 (MQ=255)
                                                                                                                                 gCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCAGGcgggtgctggt                                 <  2:4617/135‑1 (MQ=255)
                                                                                                                                 gCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCAGGcgggtgctggt                                 >  1:4617/1‑135 (MQ=255)
                                                                                                                                 gCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGAACCGACAGGGGGGCATGGTTCCAGCACCCGCCTGGTGCTGGTGCACTACCGAggcgg                  >  1:1028091/1‑150 (MQ=255)
                                                                                                                                  ccGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAggcggc                 >  1:523222/1‑150 (MQ=255)
                                                                                                                                    gATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGcca               >  2:135230/1‑149 (MQ=255)
                                                                                                                                    gATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCt               >  2:541680/1‑150 (MQ=255)
                                                                                                                                     aTCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTg              <  2:256437/150‑1 (MQ=255)
                                                                                                                                      tCGCCCAGGCACCTGCCTGGTTCTCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTGGAGCTGTGGCCAGTTTGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGc             <  1:1142597/150‑1 (MQ=255)
                                                                                                                                      tCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGc             <  2:48192/150‑1 (MQ=255)
                                                                                                                                        acccAGGCAGGTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAg           >  1:970665/2‑150 (MQ=255)
                                                                                                                                         cccAGGCAGGTGCCTGGGCCGCGGGCATTCTCGCCCATCGTGGCCAGCCGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGACAGCCGGCGGTGCAACACACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGCCGGCCTGTAGc          >  2:209867/1‑150 (MQ=255)
                                                                                                                                         cccAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGATGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGc          >  1:1327348/1‑150 (MQ=255)
                                                                                                                                          ccAGGCAGGTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGATGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCt         >  1:1085939/1‑150 (MQ=255)
                                                                                                                                                 cctgcctgGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTg  >  2:1333668/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TTGCAACTGCGCACCGCTCGGGCGGCGAACACCGAGGGGGCGTTGTACCTGCAGTTTCGCATCAGAGCGCAGCGTTTTGCCCTGAGCGTGCACGAGGTGATCGAAGTGTTGCCGCGTCAGCCGCTCAAGCCGATCGCCCAGGCACCTGCCTGGGTCGCGGGCATTCTCGCCCATCGTGGCCAGCTGGTTCCGGTCATCGACCTTGCCGCCTTGAGCTTTGGCCAGCCGGCGGTGCAACGCACCAGCACCCGCCTGGTGCTGGTGCACTACCGAGGCGGCCTGCAGCTTGGCCTG  >  NC_002947/1693583‑1693876

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: