Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,499,479 T→G 10.6% intergenic (+49/‑117) trpS → / → zapE tryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,499,4790TG10.6% 242.7 / 2.1 89intergenic (+49/‑117)trpS/zapEtryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (31/46);  new base G (4/8);  total (35/54)
Fisher's exact test for biased strand distribution p-value = 7.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAG  >  NC_002947/1499332‑1499623
                                                                                                                                                   |                                                                                                                                                 
ggCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGc                                                                                                                                                 <  2:4150391/150‑1 (MQ=255)
ggCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGc                                                                                                                                                 <  2:1004264/150‑1 (MQ=255)
  cGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCAt                                                                                                                                               <  2:1862551/150‑1 (MQ=255)
   gAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATg                                                                                                                                              <  1:4139072/150‑1 (MQ=255)
   gAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATg                                                                                                                                              <  1:3534254/150‑1 (MQ=255)
   gAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATg                                                                                                                                              >  2:3320591/1‑150 (MQ=255)
     gCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGc                                                                                                                                            <  1:400850/150‑1 (MQ=255)
       aGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCAc                                                                                                                                          <  2:2804715/150‑1 (MQ=255)
         gTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCg                                                                                                                                       <  1:851590/150‑1 (MQ=255)
         gTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCg                                                                                                                                       <  1:851578/150‑1 (MQ=255)
         gTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCg                                                                                                                                       >  2:1018412/1‑150 (MQ=255)
          gtGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCgg                                                                                                                                      <  1:3562265/149‑1 (MQ=255)
           tGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGc                                                                                                                                     <  1:2581556/150‑1 (MQ=255)
           tGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGc                                                                                                                                     <  1:2581846/150‑1 (MQ=255)
             ccGACGGGCCGCAGTTCGACCACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGACGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGCGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCCATGCACCgcggc                                                                                                                                   <  1:3054544/150‑5 (MQ=255)
              cGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCttt                                                                                                                                  <  2:3944663/150‑1 (MQ=255)
              cGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATAGCACCGGCttt                                                                                                                                  <  1:80813/150‑1 (MQ=255)
                aCGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGc                                                                                                                                <  1:2598742/150‑1 (MQ=255)
                aCGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGc                                                                                                                                <  1:582846/150‑1 (MQ=255)
                  gggCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCg                                                                                                                              <  2:3963422/150‑1 (MQ=255)
                  gggCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCg                                                                                                                              <  2:2390604/150‑1 (MQ=255)
                  gggCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCg                                                                                                                              <  2:28204/150‑1 (MQ=255)
                      cGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGt                                                                                                                          >  1:3673257‑M2/1‑149 (MQ=255)
                      cGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGGGCATGGCACCCGCGTTGCCGGTGt                                                                                                                          <  1:723818‑M2/150‑2 (MQ=255)
                          gTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgc                                                                                                                      >  2:1306689‑M2/1‑145 (MQ=255)
                           ttCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcg                                                                                                                     >  1:261301‑M2/1‑144 (MQ=255)
                           ttCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCCCGAAACAGGCGCCGCGGTGCATGGCACCCGCTTTGCCGGTGttcgcg                                                                                                                     <  2:2433534‑M2/150‑7 (MQ=255)
                            tCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcgg                                                                                                                    >  2:4521239‑M2/1‑143 (MQ=255)
                            tCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcgg                                                                                                                    >  1:1520773‑M2/1‑143 (MQ=255)
                            tCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcgg                                                                                                                    >  1:4321799‑M2/1‑143 (MQ=255)
                                  aCGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcgggcacgc                                                                                                              >  1:2621711‑M2/1‑137 (MQ=255)
                                   cGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcgggcacgcc                                                                                                             >  1:4453065‑M2/1‑136 (MQ=255)
                                    gCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCCCGCAAAAGGCGCCGCGGGGCCTGGCCCCCGCGTTGGCCGTGgtggcgggccccccc                                                                                                            <  2:471670‑M2/150‑16 (MQ=255)
                                          gCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGGGCATGGCACCGGCGTTGCCGGTGgtcgcgggcacgcccgctccc                                                                                                      <  1:3767934‑M2/150‑22 (MQ=255)
                                          gCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACACGCGCCGCGGGGCATGGCACCCGCGTTTCCGGTGttggcgggcacgcccccgccc                                                                                                      <  1:3767886‑M2/150‑22 (MQ=255)
                                                     ggCCATTCAGAGCCTGAACGTAGCCTTGCCACCACACCAGGACTGCGGCTGCCAACGGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCTGGGCACCGGCTTTGCCGGTGTGCGCGGGCACGCCCGCTCCCACAGACTGAtt                                                                                           <  1:1078798/150‑1 (MQ=25)
                                                     ggCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCCCGAAACAGGCGCCCCGGGGCAGGGCACCGGCGTTGCCGGTGTTCGCCGGCACGCCCGCACCCACAGACAGAtt                                                                                           >  2:2252251/1‑150 (MQ=255)
                                                      gCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCCCGAAACAGGCCCCGCGGGGGATGGCACCGGCTTTTCCGGTGTTAGCGGGCACCCCCGCACCCCCAGACTGATTg                                                                                          >  2:2811648/1‑150 (MQ=255)
                                                        cATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGCGGCTGCCAACCGGTGGGTGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGcc                                                                                        <  2:1824517/150‑1 (MQ=255)
                                                               agCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGc                                                                                 >  2:1422206/1‑150 (MQ=255)
                                                                gCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCg                                                                                <  2:3317108/150‑1 (MQ=255)
                                                                   ttaaCGAAGCCTTGGCACCACAGCTGGACTGTGGCTGCAAAAGGGTGGGAGGGGGCGGGCCCGCTAAACTGGCGCGGCGGGGCATGGCACCGGCTTTGCCGGTGTCCGGGGCCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGc                                                                             <  1:1820193/148‑1 (MQ=25)
                                                                   ggaccGACGCGTTGCCACCCCAGCCGTAGTGCGGGTGCCAACGGGGGGGCGGGGGCGGGCCCGCCAAACCGGGGCGGGGGTGCATGGCCCGGGCTTGGCGGGTTTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGc                                                                             <  2:405704/146‑1 (MQ=255)
                                                                       cGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCCCCCAACACGCGCCGCCGGGCCAGGCACCCGCGTTTCCCGGGGTCGCGGGGCCCCCCCCGCCCACCGACGGGTTGCCAGTAAGCggggcgat                                                                         >  2:1568448/1‑145 (MQ=255)
                                                                                ggCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCCCCAAACACGCGCCCCCGGGCATGGCACCCGCGTTGCCGGTGTTCGCGGGCCCGCCCCCGCCCACAGACCGGTTGGCCTTAAGCGGGCCGATCGCGAccgg                                                                >  2:4213838/1‑147 (MQ=255)
                                                                                ggCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGGGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACACACTGATTGCCATTAAGCGGGCCTATCGCTACAGt                                                                >  2:2611401/1‑150 (MQ=255)
                                                                                    ccacAGCAGGCCCGCGGCTGCCAAACGGTGGGAGCGGGCGGGCCCGCGAAACAGGCGCCCCGGGGCATGGCACCGGCTTTGCCGGGGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTACGCGGGCCTATCGCTACAGTGACg                                                            <  2:3010703/150‑1 (MQ=25)
                                                                                              aCTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGttttt                                                  <  2:2180273/150‑1 (MQ=255)
                                                                                                tGAGGCTGCCAATCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttg                                                <  2:3271013/150‑1 (MQ=255)
                                                                                                   ggCTGCCAACCGGTGGGCGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCAGGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgttg                                             <  2:2565420/150‑1 (MQ=255)
                                                                                                       gCCAACCGGGGGGCGGGGGGGTGCCCGCGAACCCGGCGCCGCGGGGCTTGCCCCCGGCTGGGCCGGTGCTCGCGGGCACCCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTc                                         <  1:3107816/150‑1 (MQ=255)
                                                                                                                ggggAGCGGGCGTGCCCGCGAAACCGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGaa                                <  2:3473818/149‑1 (MQ=255)
                                                                                                                 gggTGCGGGCGGGCCCGCGAAACAGGCGCCGGGGTGCAGGGCACGGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGCCTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAAt                               <  2:480281/150‑1 (MQ=38)
                                                                                                                  ggcggggggGTGCCCGCGAAACTGGCGCCGCGGTGCCTGGCACCGGCTTTGCGGGTGTCCGCGGCCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATc                              <  1:3608968/145‑1 (MQ=38)
                                                                                                                  ggcgggggCGTGCCCGCGACACTGGCCCCGCGGTGCATGCCACCGGCTTGGCCGGTGTCCGCGGGCACCCCCGCTCCCACAGCCTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATc                              <  1:3608792/145‑1 (MQ=255)
                                                                                                                  ggAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATc                              <  2:538671/150‑1 (MQ=255)
                                                                                                                   gcgggggcgggcCCGCCAACCCGGCGCGGGGGTGCTGGCCCCGGGCTTGGCGGGTGTGCGCGGCCCCCCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                             <  1:485408/146‑1 (MQ=255)
                                                                                                                    ggggggcgggcCCGCTAACCGGGCGCGGGGGGGCTGGCCCCGGGGTTCGGGGGTGTCCGCGGCCCCCCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                            <  1:1889055/147‑1 (MQ=255)
                                                                                                                         gCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTc                       <  2:2879644/150‑1 (MQ=255)
                                                                                                                          cGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTcc                      >  2:2006347/1‑150 (MQ=255)
                                                                                                                             gCCCGCTAAACTGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCtt                   <  1:3320591/150‑1 (MQ=255)
                                                                                                                               ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcgggcacgcccgctcccaccgtttggcagcctcagtcctgctgtggtgccaaggcttcgttcaggctctgaatggccgcatcacgggcttcagcgttgtcga                 <  1:144304‑M2/150‑107 (MQ=255)
                                                                                                                               ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcgggcacgcccgctcccaccgtttggcagcctcagtcctgctgtggtgccaaggcttcgttcaggctctgaatggccgcatcacgggcttcagcgttgtcga                 <  2:3435588‑M2/150‑107 (MQ=255)
                                                                                                                               ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcgggcacgcccgctcccaccgtttggcagcctcagtcctgctgtggtgccaaggcttcgttcaggctctgaatggccgcatcacgggcttcagcgttgtcga                 <  2:4345845‑M2/150‑107 (MQ=255)
                                                                                                                               ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcgggcacgcccgctcccaccgtttggcagcctcagtcctgctgtggtgccaaggcttcgttcaggctctgaatggccgcatcacgggcttcagcgttgtcga                 <  1:2353039‑M2/150‑107 (MQ=255)
                                                                                                                               ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcgggcacgcccgctcccaccgtttggcagcctcagtcctgctgtggtgccaaggcttcgttcaggctctgaatggccgcatcacgggcttcagcgttgtcga                 <  1:433864‑M2/150‑107 (MQ=255)
                                                                                                                               ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcgggcacgcccgctcccaccgtttggcagcctcagtcctgctgtggtgccaaggcttcgttcaggctctgaatggccgcatcacgggcttcagcgttgtcga                 <  1:4464755‑M2/150‑107 (MQ=255)
                                                                                                                               ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcgggcacgcccgctcccaccgtttggcagcctcagtcctgctgtggtgccaaggctccgttcaggctctgaatggccgcatcacgggcttcagcgttgtcga                 <  2:1722502‑M2/150‑107 (MQ=255)
                                                                                                                               ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAg                 >  1:264722/1‑150 (MQ=255)
                                                                                                                               ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAg                 >  2:2313881/1‑150 (MQ=255)
                                                                                                                               ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAg                 >  2:3882713/1‑150 (MQ=255)
                                                                                                                               ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTGGTTGTTGCCTCGCTAACGAATCATGACTCCCTTAg                 >  2:3326351/1‑150 (MQ=255)
                                                                                                                               ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCCCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGACTCATGACTCCCTTAg                 >  1:2877646/1‑150 (MQ=255)
                                                                                                                                cgcgAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcgggcacgcccgctcccaccgtttggcagcctcagtcctgctgtggtgccaaggcttcgttcaggctctgaatggccgcatcacgggcttcagcgttgtcgaa                <  2:3806531‑M2/150‑108 (MQ=255)
                                                                                                                                cgcgAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGa                >  1:2183149/1‑150 (MQ=255)
                                                                                                                                cgcgAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGa                <  2:1032410/150‑1 (MQ=255)
                                                                                                                                cgcgAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTGTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGa                >  2:2882727/1‑150 (MQ=255)
                                                                                                                                cgcgAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTGGTTGTTGCCTCGCTCACGACTCATGACTCCCGTAGa                >  2:2882684/1‑150 (MQ=37)
                                                                                                                                 gcgAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGaa               >  2:4415279/1‑150 (MQ=255)
                                                                                                                                    aaaCAGGAGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCAGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCGTAGAACGc            >  2:3184219/1‑150 (MQ=255)
                                                                                                                                      aCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcgggcacgcccgctcccaccgtttggcagcctcagtcctgctgtggtgccaaggcttcgttcaggctctgaatggccgcatcacgggcttcagcgttgtcgaactgcgg          >  1:4397623‑M2/1‑37 (MQ=255)
                                                                                                                                      aCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTa          >  1:3265965/1‑150 (MQ=255)
                                                                                                                                        aGGCGCCGCGGTGCATGGCACCGGCTTTGCGGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc        <  1:717280/150‑1 (MQ=255)
                                                                                                                                           cgccgcGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAg     >  2:1109845/1‑150 (MQ=255)
                                                                                                                                           cgccgcGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAg     >  1:2659381/1‑150 (MQ=255)
                                                                                                                                           aaagccgGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAg     >  1:1500358/6‑150 (MQ=255)
                                                                                                                                            gccgcGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCGTTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGGTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGc    >  2:2584521/1‑150 (MQ=255)
                                                                                                                                            gccgcGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGc    >  2:2401843/1‑150 (MQ=255)
                                                                                                                                              cgcgGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAg  >  2:302310/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GGCGAGCAGGTTGCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAG  >  NC_002947/1499332‑1499623

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