Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,823,106 G→T R236L (CGC→CTC)  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,1060GT98.7% 291.1 / NA 79R236L (CGC→CTC) spoTbifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (41/37);  total (41/38)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTAC  >  NC_000913/3822960‑3823251
                                                                                                                                                  |                                                                                                                                                 
tGAGGCGCTGTATCCCAACCGTTATCGCGTAAGCAAAGAAGTGGTGAAAGCCGGGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGGGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCga                                                                                                                                                <  1:443859/150‑1 (MQ=255)
   ggCGGTGGATCCCAACCGTTATCGCGTAATCAAAGAAGTGATGAAAGCCGCGCGCGGCAACCGTAAAGAGATGAGCCAGAAGATTCTTTCTGAAATCGAAGGGCGGTTGCAGGGAGCGGGAATACCGTGCAGCGTCAGTGGGCTCGAGaa                                                                                                                                             <  1:122661/150‑1 (MQ=255)
      gcTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCAt                                                                                                                                         >  1:1243868/1‑151 (MQ=255)
       cTGGATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCAt                                                                                                                                         <  1:49423/150‑1 (MQ=255)
        tGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCt                                                                                                                                       >  1:61207/1‑151 (MQ=255)
         gTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCt                                                                                                                                       <  1:381732/150‑1 (MQ=255)
           aTCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCttt                                                                                                                                     <  1:1193506/150‑1 (MQ=255)
             cccAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTAtt                                                                                                                                  <  2:63948/151‑1 (MQ=255)
              ccAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTc                                                                                                                                 <  1:377026/151‑1 (MQ=255)
               cAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCg                                                                                                                                >  1:1232274/1‑151 (MQ=255)
                  ccGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGAtt                                                                                                                             >  2:722703/1‑151 (MQ=255)
                   cGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGAttt                                                                                                                            <  2:1068028/151‑1 (MQ=255)
                    tttagagcgTAATCAACGAAATGGTTAAAGCCGCGAGCGGCAACGGTAAAGAGATGATCCAGAAGATTCTTTATGAAATCGAAGCGCGTTTGCAGGAAGCGGAAATACCGTGCCGCGTCGGTGGTCTCGAGCAGCATGATTATTCGAttt                                                                                                                            >  1:103811/7‑150 (MQ=255)
                    gggagagaagaagaaaaGAAGGAGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTa                                                                                                                           <  1:1196217/137‑1 (MQ=255)
                       aTCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACt                                                                                                                         <  1:984828/150‑1 (MQ=255)
                       aTCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTg                                                                                                                        >  1:795890/1‑151 (MQ=255)
                            gTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCaaa                                                                                                                    >  2:551022/1‑150 (MQ=255)
                              aaTCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGTCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAAt                                                                                                                  <  1:893792/150‑1 (MQ=255)
                              aaTCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATg                                                                                                                 >  2:550267/1‑151 (MQ=255)
                                 cAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTg                                                                                                              <  2:656677/151‑1 (MQ=255)
                                 cAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTg                                                                                                              >  2:1066960/1‑151 (MQ=255)
                                        gtggtgAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCaa                                                                                                         >  2:868519/1‑149 (MQ=255)
                                         tggtgAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCaa                                                                                                         >  1:29751/1‑148 (MQ=255)
                                           gtgAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGagc                                                                                                    <  1:53824/151‑1 (MQ=255)
                                             gAAAGCCGCGCGCGGCCACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCTTGAATACCGTGCCACGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGagca                                                                                                   <  1:890336/150‑1 (MQ=255)
                                              aaaGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGagcagc                                                                                                 <  1:451842/151‑1 (MQ=255)
                                              aaaGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGagcagc                                                                                                 >  1:592079/1‑151 (MQ=255)
                                               aaGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCg                                                                                                >  1:279867/1‑151 (MQ=255)
                                                aGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGt                                                                                               >  1:356278/1‑151 (MQ=255)
                                                       cgcgGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACt                                                                                        >  2:560430/1‑151 (MQ=255)
                                                        gcgGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTc                                                                                       >  1:107720/1‑151 (MQ=255)
                                                         cgGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCg                                                                                      <  1:496076/151‑1 (MQ=255)
                                                         cgGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCg                                                                                      >  2:779648/1‑151 (MQ=255)
                                                               ccGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATg                                                                                <  2:1225849/151‑1 (MQ=255)
                                                                   aaaGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGAc                                                                             >  2:728640/1‑150 (MQ=255)
                                                                      gagaTGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACTGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCt                                                                         <  2:37260/151‑1 (MQ=255)
                                                                          tgatCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGc                                                                     >  2:1073906/1‑151 (MQ=255)
                                                                           gatCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCt                                                                    >  2:116377/1‑151 (MQ=255)
                                                                           gatCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCt                                                                    <  2:1179249/151‑1 (MQ=255)
                                                                            atCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCtt                                                                   <  2:771239/151‑1 (MQ=255)
                                                                              ccAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTc                                                                 >  2:620002/1‑151 (MQ=255)
                                                                               aagaagaTTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGCAATACCGTGGCGCGTCAGTGGTATCGAGCATCATCTTTATGCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATAATGGACATCTAAGCtctc                                                                 >  1:759151/2‑150 (MQ=255)
                                                                                      ttctttctGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGAt                                                         >  2:871891/1‑151 (MQ=255)
                                                                                       tctttctGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGAt                                                         <  2:906817/150‑1 (MQ=255)
                                                                                         tttctGAAATCGAAGGGCGTTTGCAGGAGGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGt                                                      >  2:612929/1‑151 (MQ=255)
                                                                                          ttctGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATc                                                        <  2:909366/148‑1 (MQ=255)
                                                                                           tctGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTc                                                     >  2:173049/1‑150 (MQ=255)
                                                                                              gAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCa                                                    <  1:154104/148‑1 (MQ=255)
                                                                                              gAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATg                                                 <  2:1318636/151‑1 (MQ=255)
                                                                                              gAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGGGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTTCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACTCTTTCCGCGTGATCGTCAATg                                                 <  1:1255608/151‑1 (MQ=255)
                                                                                               aaaTCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATg                                                 >  1:903138/1‑150 (MQ=255)
                                                                                               aaaTCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATg                                                 >  1:642272/1‑150 (MQ=255)
                                                                                                aaTCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGAt                                               <  2:380778/151‑1 (MQ=255)
                                                                                                 aTCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGAtt                                              >  2:1110115/1‑151 (MQ=255)
                                                                                                  tCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTc                                             >  1:1031677/1‑151 (MQ=255)
                                                                                                    gAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTg                                           <  1:753216/151‑1 (MQ=255)
                                                                                                     aaGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTg                                           >  1:365604/1‑150 (MQ=255)
                                                                                                       gggCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGaca                                        <  2:796533/151‑1 (MQ=255)
                                                                                                       gggCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGac                                         >  2:1139245/1‑150 (MQ=255)
                                                                                                           gTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTg                                    <  2:346762/151‑1 (MQ=255)
                                                                                                             ttGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGtt                                  >  2:1087041/1‑151 (MQ=255)
                                                                                                                 aGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATcg                              >  1:897666/1‑151 (MQ=255)
                                                                                                                 aGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATcg                              >  2:719267/1‑151 (MQ=255)
                                                                                                                 aGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATcg                              <  2:606380/151‑1 (MQ=255)
                                                                                                                        gggAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTg                       >  2:836599/1‑151 (MQ=255)
                                                                                                                         ggAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTgg                      <  1:899603/151‑1 (MQ=255)
                                                                                                                            aTACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGcc                   <  2:391278/151‑1 (MQ=255)
                                                                                                                            aTACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGcc                   >  1:1234885/1‑151 (MQ=255)
                                                                                                                             tACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGcc                   >  2:889758/1‑150 (MQ=255)
                                                                                                                              aCCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGcc                   >  1:888342/1‑149 (MQ=255)
                                                                                                                              aCCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCa                  <  1:244438/150‑1 (MQ=255)
                                                                                                                              aCCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCa                  <  2:567245/150‑1 (MQ=255)
                                                                                                                              aCCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAg                 <  1:550098/151‑1 (MQ=255)
                                                                                                                                cGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGcc                   >  1:213488/1‑147 (MQ=255)
                                                                                                                                 gTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATg              >  2:819134/1‑151 (MQ=255)
                                                                                                                                            gTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTa   <  1:1310930/151‑1 (MQ=255)
                                                                                                                                              ggTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTAc  <  1:376615/150‑1 (MQ=255)
                                                                                                                                              ggTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTAc  >  1:1125146/1‑150 (MQ=255)
                                                                                                                                               gTCTCGAGAAGCATCTTTATTCGTTTTGCTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATTGAGATCTACGCTGTCCGCCTGATCGTCAATGATTCTATCACCTGTTATCGCGTTCTGGGCCAGATGCACAGCATGt    <  1:458939/147‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTAC  >  NC_000913/3822960‑3823251

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