Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,913,196 T→C K152K (AAA→AAG relA ← (p)ppGpp synthetase I/GTP pyrophosphokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,913,1960TC97.6% 304.5 / NA 83K152K (AAA→AAGrelA(p)ppGpp synthetase I/GTP pyrophosphokinase
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (41/40);  total (41/42)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

AGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGTTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATCAGTGTGCGTCGCTTTCAGCTGGCGGATCGCCGCCATATCACGCACGCCGTGAATAAGGTTAACG  >  NC_000913/2913048‑2913343
                                                                                                                                                    |                                                                                                                                                   
aGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGctt                                                                                                                                                   >  1:597952/1‑151 (MQ=255)
 gTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTtg                                                                                                                                                  <  1:561638/151‑1 (MQ=255)
     tCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTtgatg                                                                                                                                               >  1:669628/1‑150 (MQ=255)
         ccAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTa                                                                                                                                           <  1:545185/150‑1 (MQ=255)
         ccAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTAc                                                                                                                                          >  2:776538/1‑151 (MQ=255)
           aGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGc                                                                                                                                        >  1:1131969/1‑151 (MQ=255)
           aGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGc                                                                                                                                        >  1:604767/1‑151 (MQ=255)
            gTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCa                                                                                                                                       <  2:723037/151‑1 (MQ=255)
                ccATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCg                                                                                                                                    >  1:523764/1‑150 (MQ=255)
                ccATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGa                                                                                                                                   <  1:836557/151‑1 (MQ=255)
                  aTTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGaa                                                                                                                                  >  1:985890/1‑150 (MQ=255)
                   tttCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGaaa                                                                                                                                 >  1:809739/1‑150 (MQ=255)
                     tCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAatc                                                                                                                              >  2:295639/1‑151 (MQ=255)
                     tCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAat                                                                                                                               <  1:588638/150‑1 (MQ=255)
                     tCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAat                                                                                                                               <  1:517551/150‑1 (MQ=255)
                      cAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAatc                                                                                                                              >  1:212341/1‑150 (MQ=255)
                       aGTTGTCCGATTCCTAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAatcat                                                                                                                            <  2:998698/151‑1 (MQ=255)
                           gTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCg                                                                                                                          <  2:761878/149‑1 (MQ=255)
                           gTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGAc                                                                                                                        >  1:304342/1‑151 (MQ=255)
                           gTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGAc                                                                                                                        <  1:23276/151‑1 (MQ=255)
                             ccGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGAcc                                                                                                                       >  2:1250/1‑150 (MQ=255)
                                aTTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCg                                                                                                                   <  2:550232/151‑1 (MQ=255)
                                  tCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGcc                                                                                                                 >  2:829992/1‑151 (MQ=255)
                                     gagaCGTTTAGCCAGCGGTGCCTAGATGTGGGTTCACTCCTTTGCCGCCAGCACACGTTCATCTTCCGGCGCATCTTATACTTAGCGCAGATGAGCAATACGATAAGTCAGTTTGCTGAAAACACAGCGAAAAACATCGACCATCGc                                                                                                                  <  2:875134/147‑1 (MQ=255)
                                     gagaCGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAAt                                                                                                              >  1:160216/1‑151 (MQ=255)
                                      agaCGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATa                                                                                                             <  1:1163522/151‑1 (MQ=255)
                                       gaCGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATa                                                                                                             <  2:1041236/150‑1 (MQ=255)
                                         aggTTAGCCACCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATaa                                                                                                            >  2:719968/2‑149 (MQ=255)
                                          ggTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACAt                                                                                                         <  2:572002/151‑1 (MQ=255)
                                           gTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATc                                                                                                        <  2:307091/151‑1 (MQ=255)
                                           gTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATc                                                                                                        <  1:466502/151‑1 (MQ=255)
                                              aGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCgg                                                                                                     <  1:1187906/151‑1 (MQ=255)
                                              aGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCgg                                                                                                     <  2:1200021/151‑1 (MQ=255)
                                                 cAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCg                                                                                                   >  1:450165/1‑150 (MQ=255)
                                                   gCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAAc                                                                                                <  1:787911/151‑1 (MQ=255)
                                                        gCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTAt                                                                                           <  1:517176/151‑1 (MQ=255)
                                                        gCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTAt                                                                                           <  2:1131990/151‑1 (MQ=255)
                                                        gCGTAGATGTTGGTACACTCTTTTGCCGCCAGGACTCGTTCATCTTCCGGTGCATCTTTTACTTCGCGCGGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTAt                                                                                           <  2:636771/151‑1 (MQ=255)
                                                             gATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGAcc                                                                                      >  1:497999/1‑151 (MQ=255)
                                                              aTGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCt                                                                                     >  1:1239843/1‑151 (MQ=255)
                                                               tGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTg                                                                                    >  1:1110504/1‑151 (MQ=255)
                                                               tGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTg                                                                                    >  2:277330/1‑151 (MQ=255)
                                                                gTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTg                                                                                    <  2:1111777/150‑1 (MQ=255)
                                                                gTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTg                                                                                    <  2:1112230/150‑1 (MQ=255)
                                                                gTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTg                                                                                    <  2:879540/150‑1 (MQ=255)
                                                                 ttGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGtt                                                                                  >  1:754272/1‑151 (MQ=255)
                                                                  tGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGtt                                                                                  <  1:427409/150‑1 (MQ=255)
                                                                        acTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCggagg                                                                            >  2:386889/1‑150 (MQ=255)
                                                                           ctTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAAc                                                                        >  2:654675/1‑151 (MQ=255)
                                                                            ttttGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAAc                                                                        >  2:60736/1‑150 (MQ=255)
                                                                             tttGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAAc                                                                        >  2:219681/1‑149 (MQ=255)
                                                                               tGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGa                                                                     >  1:168522/1‑150 (MQ=255)
                                                                                gccgccAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAAt                                                                   <  1:715961/151‑1 (MQ=255)
                                                                                  cgccAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATc                                                                  >  2:847518/1‑150 (MQ=255)
                                                                                       gTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATCAGtgtg                                                            <  2:433393/151‑1 (MQ=255)
                                                                                                cATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATCAGTGTGCGTCGCttt                                                   >  2:907724/1‑151 (MQ=255)
                                                                                                cATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATCAGTGTGCGTCGCtt                                                    <  2:435160/150‑1 (MQ=255)
                                                                                                 aTCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATCAGTGTGCGTCGCTTTc                                                  >  1:511147/1‑151 (MQ=255)
                                                                                                  tCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATCAGTGTGCGTCGCTTTc                                                  >  2:880352/1‑150 (MQ=255)
                                                                                                     tCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATCAGTGTGCGTCGCTTTc                                                  <  1:424793/147‑1 (MQ=255)
                                                                                                     tCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATCAGTGTGCGTCGCTTTCAg                                                <  2:212353/149‑1 (MQ=255)
                                                                                                     tCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATCAGTGTGCGTCGCTTTCAGc                                               <  2:385056/150‑1 (MQ=255)
                                                                                                     tCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATCAGTGTGCGTCGCTTTCAGCt                                              >  2:269850/1‑151 (MQ=255)
                                                                                                     tCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATCAGTGTGCGTCGCTTTCAGCt                                              >  1:87879/1‑151 (MQ=255)
                                                                                                                 tttACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATCAGTGTGCGTCGCTTTCAGCTGGCGGATcgcc                                   <  1:994417/150‑1 (MQ=255)
                                                                                                                 tttACTTCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATCAGTGTGCGTCGCTTTCAGCTGGCGGATc                                      <  2:809732/147‑1 (MQ=255)
                                                                                                                       tCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATCAGTGTGCGTCGCTTTCAGCTGGCGGATCGCCGCCat                              <  1:1250/149‑1 (MQ=255)
                                                                                                                       tCGCGCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATCAGTGTGCGTCGCTTTCAGCTGGCGGATCGCCGCCat                              <  2:142433/149‑1 (MQ=255)
                                                                                                                        cgcgCAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATCAGTGTGCGTCGCTTTCAGCTGGCGGATCGCCGCCatat                            <  1:1088644/150‑1 (MQ=255)
                                                                                                                           gcAGATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATCAGTGTGCGTCGCTTTCAGCTGGCGGATCGCCGCCATATcacg                        >  2:1150160/1‑151 (MQ=255)
                                                                                                                              gATGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATCAGTGTGCGTCGCTTTCAGCTGGCGGATCGCCGCCATATcacgcac                     <  1:282908/151‑1 (MQ=255)
                                                                                                                               aTGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATCAGTGTGCGTCGCTTTCAGCTGGCGGATCGCCGCCATATcacgcacg                    >  1:113337/1‑151 (MQ=255)
                                                                                                                               aTGAGCAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATCAGTGTGCGTCGCTTTCAGCTGGCGGATCGCCGCCATATcacgcacg                    >  1:506531/1‑151 (MQ=255)
                                                                                                                                    cAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATCAGTGTGCGTCGCTTTCAGCTGGCGGATCGCCGCCATATCACGCACGCCGTg               >  2:1083442/1‑151 (MQ=255)
                                                                                                                                    cAATACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATCAGTGTGCGTCGCTTTCAGCTGGCGGATCGCCGCCATATCACGCACGCCGTg               >  2:1229564/1‑151 (MQ=255)
                                                                                                                                     aaTACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATCAGTGTGCGTCGCTTTCAGCTGGCGGATCGCCGCCATATCACGCACGCCGTg               <  1:1046870/150‑1 (MQ=255)
                                                                                                                                     aaTACGCTCCGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATCAGTGTGCGTCGCTTTCAGCTGGCGGATCGCCGCCATATCACGCACGCCGTGa              <  1:1059726/151‑1 (MQ=255)
                                                                                                                                         cGCTCCGCCAGTTTGATGACTCCGCAGCGAAAATCATCGACCATCCCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATCAGTGTGGGTCGCTTTCAGCTGGCGGATCGCCGCCATATcacaaatcggaagagcgt          <  2:251850/151‑16 (MQ=255)
                                                                                                                                             ccGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATCAGTGTGCGTCGCTTTCAGCTGGCGGATCGCCGCCATATCACGCACGCCGTGAATAAGGtt      >  1:122445/1‑151 (MQ=255)
                                                                                                                                             ccGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATCAGTGTGCGTCGCTTTCAGCTGGCGGATCGCCGCCATATCACGCACGCCGTGAATAAGGtt      >  1:674861/1‑151 (MQ=255)
                                                                                                                                             ccGCCAGCTTGATGACTACGCAACGAAAATCACCGACCTTCGCCAAAAACACCCGGCGAACGTTATCGGCCTGTTCGGAGAAAACAGAATCAGTGTGCGTCGTTTTCAGCTGGCGGTTCGCCGCCATATCACGCACGCCGTGAATAAGGt       <  2:669611/150‑1 (MQ=255)
                                                                                                                                              cGCCAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATCAGTGTGCGTCGCTTTCAGCTGGCGGATCGCCGCCATATCACGCACGCCGTGAATAAGGtt      >  1:844517/1‑150 (MQ=255)
                                                                                                                                                 cAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATCAGTGTGCGTCGCTTTCAGCTGGCGGATCGCCGCCATATCACGCACGCCGTGAATAAGGTTAAcg  >  1:933799/1‑151 (MQ=255)
                                                                                                                                                 cAGCTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATCAGTGTGCGTCGCTTTCAGCTGGCGGATCGCCGCCATATCACGCACGCCGTGAATAAGGTTAAcg  >  2:1202140/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AGTAATCTTCCAGTTCCCATTTCAGTTGTCCGATTCCGAGACGGTTAGCCAGCGGTGCGTAGATGTTGGTACACTCTTTTGCCGCCAGTACACGTTCATCTTCCGGCGCATCTTTTACTTCGCGCAGATGAGCAATACGCTCCGCCAGTTTGATGACTACGCAGCGAAAATCATCGACCATCGCCAATAACATCCGGCGAACGTTATCGACCTGTTCGGAGGAAACAGAATCAGTGTGCGTCGCTTTCAGCTGGCGGATCGCCGCCATATCACGCACGCCGTGAATAAGGTTAACG  >  NC_000913/2913048‑2913343

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