Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,823,025 G→A R209H (CGT→CAT)  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,0250GA100.0% 259.0 / NA 65R209H (CGT→CAT) spoTbifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (30/35);  total (30/35)

CGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCAT  >  NC_000913/3822877‑3823172
                                                                                                                                                    |                                                                                                                                                   
cGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATa                                                                                                                                                   <  2:558068/151‑1 (MQ=255)
cGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATa                                                                                                                                                   <  1:305922/151‑1 (MQ=255)
         ctcGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAgat                                                                                                                                            <  2:338754/149‑1 (MQ=255)
           cGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAgatgat                                                                                                                                         <  1:162712/150‑1 (MQ=255)
                  tataGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCagaag                                                                                                                                  <  1:410595/150‑1 (MQ=255)
                   ataGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCagaaga                                                                                                                                 <  2:231395/150‑1 (MQ=255)
                      gCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGAttct                                                                                                                             >  2:773358/1‑151 (MQ=255)
                       cccGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGAttctt                                                                                                                            <  1:648131/151‑1 (MQ=255)
                       cccGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGAttctt                                                                                                                            >  2:69787/1‑151 (MQ=255)
                         cGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGAttctttc                                                                                                                          <  1:463905/151‑1 (MQ=255)
                            tGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTg                                                                                                                        >  2:225957/1‑150 (MQ=255)
                             ggCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGaa                                                                                                                      <  1:667438/151‑1 (MQ=255)
                                 cACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATc                                                                                                                   <  2:289464/150‑1 (MQ=255)
                                  aCCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCg                                                                                                                  >  2:325846/1‑150 (MQ=255)
                                       ttAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAggg                                                                                                             <  1:448478/150‑1 (MQ=255)
                                        tAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCg                                                                                                           >  2:313465/1‑151 (MQ=255)
                                         aGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGt                                                                                                          >  2:387590/1‑151 (MQ=255)
                                          ggTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGtt                                                                                                         <  2:559119/151‑1 (MQ=255)
                                          ggTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGtt                                                                                                         <  2:840526/151‑1 (MQ=255)
                                             aTCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTg                                                                                                       <  2:1086170/150‑1 (MQ=255)
                                             aTCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTg                                                                                                       <  1:1055332/150‑1 (MQ=255)
                                                      aTTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGAGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGGGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGACATCGAAGGGCGGTTGCAGGAAGCgg                                                                                             >  1:989543/1‑151 (MQ=255)
                                                           aaCCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATa                                                                                        <  1:184889/151‑1 (MQ=255)
                                                             ccGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATAcc                                                                                      >  2:1038790/1‑151 (MQ=255)
                                                             ccGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATAcc                                                                                      <  2:675259/151‑1 (MQ=255)
                                                                  cTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGcc                                                                                 >  1:64027/1‑151 (MQ=255)
                                                                    cGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCcgc                                                                               >  1:193117/1‑151 (MQ=255)
                                                                    cGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCcgc                                                                               >  2:88169/1‑151 (MQ=255)
                                                                     gAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCcgcg                                                                              >  1:1092395/1‑151 (MQ=255)
                                                                       agagCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTc                                                                            >  1:700549/1‑151 (MQ=255)
                                                                        gagCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCa                                                                           <  2:1039546/151‑1 (MQ=255)
                                                                         agCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTc                                                                            >  2:119659/1‑149 (MQ=255)
                                                                         agCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAg                                                                          <  1:46388/151‑1 (MQ=255)
                                                                             gggTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTgg                                                                       <  1:1071548/150‑1 (MQ=255)
                                                                                  ttGAGGTGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTcgcg                                                                  >  2:956523/1‑150 (MQ=255)
                                                                                    gAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCgaga                                                               <  1:190310/151‑1 (MQ=255)
                                                                                     aGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTTGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCgatag                                                              <  1:1159777/151‑4 (MQ=255)
                                                                                        cgcTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGc                                                            >  1:549479/1‑150 (MQ=255)
                                                                                         gcTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCAt                                                          <  2:459249/151‑1 (MQ=255)
                                                                                         gcTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCAt                                                          >  1:1174448/1‑151 (MQ=255)
                                                                                           tGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGTTGATCCAGAAGATTCTTCCTGAAATCGAAGGGGGTTTGCAGGAAGCGGGAATACCGTGCCGCGCCAGTGGTCGCGAGAAGCATCt                                                        >  2:933837/1‑151 (MQ=255)
                                                                                           tGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCt                                                        >  1:860392/1‑151 (MQ=255)
                                                                                            gTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCt                                                        <  1:384439/150‑1 (MQ=255)
                                                                                               tCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTAt                                                    >  2:1112066/1‑151 (MQ=255)
                                                                                               tCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTAt                                                    >  2:1120990/1‑151 (MQ=255)
                                                                                                   aaCCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCg                                                 >  2:1014652/1‑150 (MQ=255)
                                                                                                     ccGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGAtt                                              >  2:1073692/1‑151 (MQ=255)
                                                                                                     ccGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGAtt                                              <  2:507898/151‑1 (MQ=255)
                                                                                                      cGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGAttt                                             <  1:607606/151‑1 (MQ=255)
                                                                                                      cGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGAttt                                             <  2:64436/151‑1 (MQ=255)
                                                                                                       gTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTa                                            <  1:168401/151‑1 (MQ=255)
                                                                                                          aTCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACt                                          <  2:1172543/150‑1 (MQ=255)
                                                                                                          aTCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTg                                         >  1:705038/1‑151 (MQ=255)
                                                                                                            cgcgTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGc                                        >  1:1138741/1‑150 (MQ=255)
                                                                                                              cgTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCaaa                                     <  1:801236/151‑1 (MQ=255)
                                                                                                               gTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCaaaa                                    <  1:193107/151‑1 (MQ=255)
                                                                                                               gTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCaaaa                                    <  1:227178/151‑1 (MQ=255)
                                                                                                                    cAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGt                                <  2:1010572/150‑1 (MQ=255)
                                                                                                                    cAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGCGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGt                                <  2:743186/150‑1 (MQ=255)
                                                                                                                        gaagTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTc                            >  1:594841/1‑150 (MQ=255)
                                                                                                                         aagTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCa                           >  1:761698/1‑150 (MQ=255)
                                                                                                                               tgAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGagc                     >  1:1075891/1‑150 (MQ=255)
                                                                                                                                gAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGagca                    <  2:1066047/150‑1 (MQ=255)
                                                                                                                                   aGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGt                >  1:1010704/1‑151 (MQ=255)
                                                                                                                                                 aCCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCAt  >  2:1137508/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCAT  >  NC_000913/3822877‑3823172

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