Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,823,139 T→A M247K (ATG→AAG)  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,1390TA99.0% 367.4 / NA 96M247K (ATG→AAG) spoTbifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (65/30);  minor base C (1/0);  total (66/30)
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.98e-01

ATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATC  >  NC_000913/3822991‑3823284
                                                                                                                                                    |                                                                                                                                                 
aTCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAgg                                                                                                                                                 <  1:313793/151‑1 (MQ=255)
 tCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGt                                                                                                                                                >  1:813849/1‑151 (MQ=255)
  cAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGt                                                                                                                                                <  2:406713/150‑1 (MQ=255)
  cAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTg                                                                                                                                               >  2:185303/1‑151 (MQ=255)
   aaaGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGc                                                                                                                                              >  1:954192/1‑151 (MQ=255)
    aagaagTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCt                                                                                                                                             >  1:1072856/1‑151 (MQ=255)
     agaagTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTc                                                                                                                                            >  2:753627/1‑151 (MQ=255)
      gaagTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTc                                                                                                                                            >  2:843967/1‑150 (MQ=255)
      gaagTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTc                                                                                                                                            >  1:805863/1‑150 (MQ=255)
       aagTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCa                                                                                                                                           >  2:1187165/1‑150 (MQ=255)
           ggtgAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAagag                                                                                                                                      <  2:370282/151‑1 (MQ=255)
               aaaGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGagcagc                                                                                                                                  >  1:1344206/1‑151 (MQ=255)
                   ccGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTAATGCAAAAAGGTGCTAAAAGAGCAGCGttt                                                                                                                              >  2:815361/1‑151 (MQ=255)
                    cgcgcgcgGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGtttt                                                                                                                             >  1:1235208/1‑151 (MQ=255)
                     ccgcgcgGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGtttt                                                                                                                             >  2:294822/2‑150 (MQ=255)
                              aaCCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCGTCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATc                                                                                                                    >  2:332737/1‑150 (MQ=255)
                              aaCCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATc                                                                                                                    >  1:95533/1‑150 (MQ=255)
                              aaCCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCa                                                                                                                   <  1:342500/151‑1 (MQ=255)
                               aCCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCAt                                                                                                                  >  1:46615/1‑151 (MQ=255)
                                 cGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATgg                                                                                                                >  1:256655/1‑151 (MQ=255)
                                 cGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCATGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATgg                                                                                                                <  1:1017751/151‑1 (MQ=255)
                                  gTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATgg                                                                                                                >  1:1044861/1‑150 (MQ=255)
                                    aaaGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTCTTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGAc                                                                                                              >  2:330424/1‑150 (MQ=255)
                                    aaaGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACa                                                                                                             <  1:86621/151‑1 (MQ=255)
                                     aaGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACAt                                                                                                            >  1:208169/1‑151 (MQ=255)
                                        agaTGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCt                                                                                                          >  1:491859/1‑150 (MQ=255)
                                        agaTGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCt                                                                                                          >  2:170540/1‑150 (MQ=255)
                                        agaTGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTa                                                                                                         <  1:217355/151‑1 (MQ=255)
                                          atgatCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACg                                                                                                       >  2:696291/1‑151 (MQ=255)
                                           tgatCCAGAAGATTCTTGTTGAAATCGAAGGGCGTTTGCAGGACGCGGGAATACCGTGCCGCGTCAGTGGGCGGGAGAAGCATCTTTATTCGATTTACTGCAAAACGGTGCTCAAAGAGCAGCGTTTTGCCTCGATCATGGATATCTAc                                                                                                        >  2:580939/1‑149 (MQ=255)
                                            gatCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCt                                                                                                     >  1:470973/1‑151 (MQ=255)
                                            gatCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCt                                                                                                     >  2:199241/1‑151 (MQ=255)
                                             ctCCAGAAGATTCTTTCTGAAATCGAAGGGCGAAGGCAGGAAGCAGGAATACGGTACTTCGTAATTGGTCGCGAGAAGCATCTATATTCGATTTACTGCAAACAGGAGCTCAAAGATCAGCGTCTTCACTCGATTTACGACCTCTACGcat                                                                                                    >  2:1173776/2‑149 (MQ=255)
                                             atCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCtt                                                                                                    >  1:550374/1‑151 (MQ=255)
                                             atCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCtt                                                                                                    >  1:41029/1‑151 (MQ=255)
                                              tCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCttt                                                                                                   >  2:642854/1‑151 (MQ=255)
                                              tCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCttt                                                                                                   >  1:413337/1‑151 (MQ=255)
                                               ccAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTc                                                                                                  >  2:1356545/1‑151 (MQ=255)
                                               ccAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAAAACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCttt                                                                                                   <  1:705835/150‑1 (MQ=255)
                                                cAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTcc                                                                                                 >  2:574690/1‑151 (MQ=255)
                                                cAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTcc                                                                                                 >  2:407381/1‑151 (MQ=255)
                                                 agaagaTTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCcg                                                                                                >  2:18519/1‑151 (MQ=255)
                                                         ctttctGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCg                                                                                        >  2:1249488/1‑151 (MQ=255)
                                                           ttctGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATc                                                                                         <  1:987828/148‑1 (MQ=255)
                                                           ttctGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGt                                                                                       <  2:872427/150‑1 (MQ=255)
                                                               gAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCa                                                                                     <  2:654305/148‑1 (MQ=255)
                                                                aaaTCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGa                                                                                 <  2:453599/151‑1 (MQ=255)
                                                                 caTCGAAGCGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGAt                                                                                >  2:419983/2‑151 (MQ=255)
                                                                   tCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTc                                                                              >  1:1303776/1‑151 (MQ=255)
                                                                    cGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCt                                                                             >  2:832856/1‑151 (MQ=255)
                                                                     gAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTg                                                                            >  1:357616/1‑151 (MQ=255)
                                                                      aaGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGa                                                                           <  1:1106513/151‑1 (MQ=255)
                                                                        gggCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGaca                                                                         <  1:922264/151‑1 (MQ=255)
                                                                          gCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACAcc                                                                       >  1:943639/1‑151 (MQ=255)
                                                                           cGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCt                                                                      >  1:547007/1‑151 (MQ=255)
                                                                            gTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTg                                                                     <  2:1072189/151‑1 (MQ=255)
                                                                            gTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTg                                                                     >  1:478720/1‑151 (MQ=255)
                                                                                        cGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCt                                                         <  1:844553/151‑1 (MQ=255)
                                                                                         gggAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTg                                                        <  2:83332/151‑1 (MQ=255)
                                                                                           gAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTggg                                                      <  1:623730/151‑1 (MQ=255)
                                                                                              tACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGcc                                                    >  2:1284114/1‑150 (MQ=255)
                                                                                                          tCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCt                                       >  2:308859/1‑151 (MQ=255)
                                                                                                           cAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTg                                      <  1:185355/151‑1 (MQ=255)
                                                                                                           cAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTg                                      <  2:1310732/151‑1 (MQ=255)
                                                                                                            aGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGt                                     >  1:1149886/1‑151 (MQ=255)
                                                                                                             gTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTa                                    <  1:757337/151‑1 (MQ=255)
                                                                                                               ggTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTAc                                   >  2:1075261/1‑150 (MQ=255)
                                                                                                               ggTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTAc                                   >  2:233214/1‑150 (MQ=255)
                                                                                                                 tCGCGAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACa                                  >  1:549534/1‑149 (MQ=255)
                                                                                                                  cgcgAGAAGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAg                                <  1:170546/150‑1 (MQ=255)
                                                                                                                  cgcgAGAAGCATCTTTAGTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAg                                >  2:93125/1‑150 (MQ=255)
                                                                                                                        aaGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTc                         >  2:387830/1‑151 (MQ=255)
                                                                                                                        aaGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTc                         >  1:947807/1‑151 (MQ=255)
                                                                                                                         aGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTcc                        >  2:316690/1‑151 (MQ=255)
                                                                                                                         aGCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTcc                        <  1:361872/151‑1 (MQ=255)
                                                                                                                          gCATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCg                       >  1:568251/1‑151 (MQ=255)
                                                                                                                           cATCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCgg                      >  2:388400/1‑151 (MQ=255)
                                                                                                                            aTCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCggg                     >  2:1051378/1‑151 (MQ=255)
                                                                                                                            aTCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCg                       >  2:1244840/1‑149 (MQ=255)
                                                                                                                             tCTTTATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCgg                      <  2:1002428/149‑1 (MQ=255)
                                                                                                                               tttATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCcg                  >  2:1005047/1‑151 (MQ=255)
                                                                                                                               tttATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCcg                  >  2:623334/1‑151 (MQ=255)
                                                                                                                               tttATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCcg                  >  1:804880/1‑151 (MQ=255)
                                                                                                                                ttATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGcc                   <  2:813278/149‑1 (MQ=255)
                                                                                                                                ttATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGcc                   <  2:805296/149‑1 (MQ=255)
                                                                                                                                ttATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCcgc                 >  1:1083126/1‑151 (MQ=255)
                                                                                                                                ttATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCcgc                 >  2:445766/1‑151 (MQ=255)
                                                                                                                                 tATTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCcgcg                >  2:1328475/1‑151 (MQ=255)
                                                                                                                                  aTTCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCcgcg                <  2:273932/150‑1 (MQ=255)
                                                                                                                                    tCGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTg              <  2:1189316/150‑1 (MQ=255)
                                                                                                                                     cGATTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGa             <  2:491748/150‑1 (MQ=255)
                                                                                                                                       aTTTACTGCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGaa            >  1:987244/1‑149 (MQ=255)
                                                                                                                                              gCAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCTGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTatat   <  2:321219/151‑1 (MQ=255)
                                                                                                                                               cAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATc  >  2:1119327/1‑151 (MQ=255)
                                                                                                                                               cAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATc  >  1:1098128/1‑151 (MQ=255)
                                                                                                                                               cAAAAAGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATc  >  2:1174670/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
ATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATC  >  NC_000913/3822991‑3823284

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