Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,823,664 T→C V422A (GTG→GCG)  spoT → bifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,823,6640TC98.4% 220.2 / ‑5.0 64V422A (GTG→GCG) spoTbifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (37/26);  total (38/26)
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.61e-01

TTTGAATTTATCGAGAGCGTTAAATCCGATCTCTTCCCGGATGAGATTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTC  >  NC_000913/3823516‑3823810
                                                                                                                                                    |                                                                                                                                                  
tttGAATTTATCGAGAGCGTTAAATCCGATCTCTTCCCGGATGAGATTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCgg                                                                                                                                                  >  1:589446/1‑151 (MQ=255)
tttGAATTTATCGAGAGCGTTAAATCCGATCTCTTCCCGGATGAGATTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCgg                                                                                                                                                  >  1:731956/1‑151 (MQ=255)
 ttGAATTTATCGAGAGCGTTAAATCCGATCTCTTCCCGGATGAGATTTACGTTTTCACACCGGACGGGCGCATTGTCGAGCTGCCAGCCGGTGCAACGCCCGTCGACTTCGCTTATGCCGCGCATACCGATATCGGTCCTGCCTGCGCggg                                                                                                                                                 >  1:468777/1‑151 (MQ=255)
   gAATTTATCGAGAGCGTTAAATCCGATCTCTTCCCGGATGAGATTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGgcg                                                                                                                                               >  1:182450/1‑151 (MQ=255)
    aaTTTATCGAGAGCGTTAAATCCGATCTCTTCCCGGATGAGATTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGgcgc                                                                                                                                              >  1:40796/1‑151 (MQ=255)
    aaTTTATCGAGAGCGTTAAATCCGATCTCTTCCCGGATGAGATTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGgcgc                                                                                                                                              >  2:136577/1‑151 (MQ=255)
    aaTTTATCGAGAGCGTTAAATCCGATCTCTTCCCGGATGAGATTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGgcgc                                                                                                                                              >  2:205953/1‑151 (MQ=255)
        tATCGAGAGCGTTAAATCCGATCTCTTCCCGGATGAGATTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCAcgc                                                                                                                                          >  1:422844/1‑151 (MQ=255)
         aTCGAGAGCGCTAAATCCGATCTCTTCCCGGATGAGATTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCAcgcg                                                                                                                                         >  1:231525/1‑151 (MQ=255)
           cGAGAGCGTTAAATCCGATCTCTTCCCGGATGAGATTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGtt                                                                                                                                       <  2:734490/151‑1 (MQ=255)
                    tAAATCCGATCTCTTCCCGGATGAGATTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGcc                                                                                                                                >  2:559288/1‑149 (MQ=255)
                     aaaTCCGATCTCTTCCCGGATGAGATTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGc                                                                                                                             >  1:578173/1‑151 (MQ=255)
                             tctcTTCCCGGATGAGATTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCg                                                                                                                     >  1:310644/1‑151 (MQ=255)
                             tctcTTCCCGGATGAGATTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCg                                                                                                                     >  1:899484/1‑151 (MQ=255)
                               tctTCCCGGATGAGATTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCt                                                                                                                   >  1:124258/1‑151 (MQ=255)
                                ctTCCCGGATGAGATTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTAccc                                                                                                                      <  2:910734/147‑1 (MQ=255)
                                ctTCCCGGATGAGATTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGc                                                                                                                    <  2:502305/149‑1 (MQ=255)
                                   cccGGATGAGATTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCg                                                                                                               <  1:762017/151‑1 (MQ=255)
                                       gATGAGATTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGc                                                                                                           >  2:580908/1‑151 (MQ=255)
                                             aTTTACGTTGTCACACCGGAGGGGCGCATGGTCGAGCGGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATCCCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTa                                                                                                     <  1:137558/151‑1 (MQ=255)
                                                   gTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCg                                                                                               >  1:973393/1‑151 (MQ=255)
                                                         acacCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCaaa                                                                                         <  2:899842/151‑1 (MQ=255)
                                                         acacCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCaa                                                                                          >  1:456553/1‑150 (MQ=255)
                                                          cacCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCaaa                                                                                         <  1:736374/150‑1 (MQ=255)
                                                          cacCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAAc                                                                                        >  2:152077/1‑151 (MQ=255)
                                                                aaGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTg                                                                                   >  1:1082371/1‑150 (MQ=255)
                                                                aaGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGa                                                                                  <  2:987375/151‑1 (MQ=255)
                                                                  gggCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGaaa                                                                                <  1:840939/151‑1 (MQ=255)
                                                                   ggCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCCTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAAt                                                                               <  1:929145/151‑1 (MQ=255)
                                                                       cATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCAtt                                                                           >  2:980359/1‑151 (MQ=255)
                                                                        aTTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCAtt                                                                           >  1:551419/1‑150 (MQ=255)
                                                                        aTTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTa                                                                          <  2:149710/151‑1 (MQ=255)
                                                                              gAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTc                                                                    <  2:87363/151‑1 (MQ=255)
                                                                               aGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTcc                                                                   >  1:397956/1‑151 (MQ=255)
                                                                                   gcctgccGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGGCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCgg                                                                 <  1:244606/149‑1 (MQ=255)
                                                                                              cAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGa                                                     >  2:666969/1‑150 (MQ=255)
                                                                                              cAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGa                                                     >  2:343651/1‑150 (MQ=255)
                                                                                                aCGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATg                                                  >  2:738736/1‑151 (MQ=255)
                                                                                                    ccGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGc                                              >  2:649523/1‑151 (MQ=255)
                                                                                                      gTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCtt                                            <  1:945728/151‑1 (MQ=255)
                                                                                                        cGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGCCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTg                                           >  1:144733/1‑150 (MQ=255)
                                                                                                          aCTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCt                                        <  1:902276/151‑1 (MQ=255)
                                                                                                           cTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTg                                       >  1:157471/1‑151 (MQ=255)
                                                                                                            ttCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTg                                       >  1:539431/1‑150 (MQ=255)
                                                                                                            ttCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTg                                       >  1:619831/1‑150 (MQ=255)
                                                                                                                 ttATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACttt                                 >  2:878167/1‑151 (MQ=255)
                                                                                                                  tATGCAGTGCCTAACGATATCGGTCATGCCTGCGCGGGCTGACGCGCTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTGACCAGCGGTCAAACCGTTGAAATCATTACCGCCCCCGGCCCGCGCCCGAATGCCGCTTGGCTGAACt                                   >  1:352816/1‑148 (MQ=255)
                                                                                                                  tATGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACttt                                 <  1:1065417/150‑1 (MQ=255)
                                                                                                                   aTGCAGTGCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGGCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTTCCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACttt                                 <  1:148044/149‑1 (MQ=255)
                                                                                                                          gCATACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAg                         <  2:874417/150‑1 (MQ=255)
                                                                                                                             tACCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCg                     >  2:279306/1‑151 (MQ=255)
                                                                                                                              aCCGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCg                     <  1:391800/150‑1 (MQ=255)
                                                                                                                                cGATATCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGaaa                  <  1:502929/151‑1 (MQ=255)
                                                                                                                                  atatCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAgc                >  1:477649/1‑151 (MQ=255)
                                                                                                                                    atCGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAgcgc              <  2:73629/151‑1 (MQ=255)
                                                                                                                                      cGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTCCCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAgcgcg             <  1:167203/150‑1 (MQ=255)
                                                                                                                                      cGGTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAgcgcgc            <  2:1008851/151‑1 (MQ=255)
                                                                                                                                       ggTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAgcgcgcg           >  2:11544/1‑151 (MQ=255)
                                                                                                                                       ggTCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAgcgcgcg           <  2:777292/151‑1 (MQ=255)
                                                                                                                                         tCATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGcc         >  2:1093282/1‑151 (MQ=255)
                                                                                                                                          cATGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGcc         >  2:1071326/1‑150 (MQ=255)
                                                                                                                                           aTGCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCaa       <  2:1077497/151‑1 (MQ=255)
                                                                                                                                             gCCTGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCaaa      <  1:104331/150‑1 (MQ=255)
                                                                                                                                                tGCGCGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTc  >  1:21879/1‑151 (MQ=255)
                                                                                                                                                 gcgcGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAAtt   <  1:1084998/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TTTGAATTTATCGAGAGCGTTAAATCCGATCTCTTCCCGGATGAGATTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGCTGCCTGCCGGTGCAACGCCCGTCGACTTCGCTTATGCAGTGCATACCGATATCGGTCATGCCTGCGTGGGCGCACGCGTTGACCGCCAGCCTTACCCGCTGTCGCAGCCGCTTACCAGCGGTCAAACCGTTGAAATCATTACCGCTCCGGGCGCTCGCCCGAATGCCGCTTGGCTGAACTTTGTCGTTAGCTCGAAAGCGCGCGCCAAAATTC  >  NC_000913/3823516‑3823810

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: