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% 347.6 / NA 89R209H (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 (47/42);  total (47/42)

GTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCAT  >  NC_000913/3822878‑3823172
                                                                                                                                                   |                                                                                                                                                   
gTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATa                                                                                                                                                   >  1:1189920/1‑150 (MQ=255)
  gAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAag                                                                                                                                                <  1:1045608/151‑1 (MQ=255)
     aCTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAag                                                                                                                                                >  1:557046/1‑148 (MQ=255)
     aCTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAag                                                                                                                                                >  1:887446/1‑148 (MQ=255)
        ctcGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAgatg                                                                                                                                           >  1:733642/1‑150 (MQ=255)
        ctcGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAgat                                                                                                                                            <  1:379705/149‑1 (MQ=255)
         tcGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAgatgat                                                                                                                                         >  2:1140289/1‑151 (MQ=255)
          cGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAgatga                                                                                                                                          <  1:1174357/149‑1 (MQ=255)
            aaaTTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATcc                                                                                                                                       >  1:436314/1‑150 (MQ=255)
            aaaTTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCa                                                                                                                                      <  2:351228/151‑1 (MQ=255)
            aaaTTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCa                                                                                                                                      <  2:304421/151‑1 (MQ=255)
               tttATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCaga                                                                                                                                    >  1:532356/1‑150 (MQ=255)
                   taGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGAtt                                                                                                                               >  2:216928/1‑151 (MQ=255)
                    aGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGAtt                                                                                                                               <  2:18294/150‑1 (MQ=255)
                    aGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGAtt                                                                                                                               <  1:204410/150‑1 (MQ=255)
                          cTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTg                                                                                                                        >  2:371099/1‑151 (MQ=255)
                            ggCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGaa                                                                                                                      <  1:20771/151‑1 (MQ=255)
                             gcgcACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGaaa                                                                                                                     <  2:127372/151‑1 (MQ=255)
                              cgcACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAAt                                                                                                                    >  1:895308/1‑151 (MQ=255)
                               gcACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATc                                                                                                                   <  1:919957/151‑1 (MQ=255)
                                cACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCg                                                                                                                  <  1:505486/151‑1 (MQ=255)
                                      ttAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAggg                                                                                                             <  2:697851/150‑1 (MQ=255)
                                        aGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGt                                                                                                          <  1:469307/151‑1 (MQ=255)
                                         ggTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGtt                                                                                                         <  1:1006655/151‑1 (MQ=255)
                                          gTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGttt                                                                                                        >  1:159029/1‑151 (MQ=255)
                                           tATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTg                                                                                                       >  2:939558/1‑151 (MQ=255)
                                             tCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGc                                                                                                      >  1:1189766/1‑150 (MQ=255)
                                               caccacATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAgg                                                                                                   >  2:656271/1‑151 (MQ=255)
                                                 ccacATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGa                                                                                                  >  2:736692/1‑150 (MQ=255)
                                                    caTTAAAGCCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTAATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGc                                                                                               <  2:1201658/150‑1 (MQ=255)
                                                    caTTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCg                                                                                              >  2:436946/1‑151 (MQ=255)
                                                     aTTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCgg                                                                                             <  2:887402/151‑1 (MQ=255)
                                                       tAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCggg                                                                                            <  1:291525/150‑1 (MQ=255)
                                                         aaaCCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAAt                                                                                         <  1:1146717/151‑1 (MQ=255)
                                                          aaCCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAAt                                                                                         >  2:1045163/1‑150 (MQ=255)
                                                           aCCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATAc                                                                                       >  2:758933/1‑151 (MQ=255)
                                                              gAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCg                                                                                     >  1:919601/1‑150 (MQ=255)
                                                              gAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGt                                                                                    >  2:402082/1‑151 (MQ=255)
                                                                 cTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGc                                                                                  >  1:182098/1‑150 (MQ=255)
                                                                   cGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCcgc                                                                               >  1:71671/1‑151 (MQ=255)
                                                                   cGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCcgc                                                                               >  1:972623/1‑151 (MQ=255)
                                                                         gCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGt                                                                         >  1:535534/1‑151 (MQ=255)
                                                                            gggTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTgg                                                                       >  1:1225076/1‑150 (MQ=255)
                                                                             ggTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTc                                                                     <  1:1239588/151‑1 (MQ=255)
                                                                             ggTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTc                                                                     >  1:16462/1‑151 (MQ=255)
                                                                                tttGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTcgc                                                                   <  2:159251/150‑1 (MQ=255)
                                                                                  tGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCga                                                                 <  1:756806/150‑1 (MQ=255)
                                                                                   gAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGTAATACCGTGCCGCGTCAGTGGTCGCgaga                                                               <  2:68014/151‑1 (MQ=255)
                                                                                         cTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATc                                                         >  1:27933/1‑151 (MQ=255)
                                                                                           gTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCtt                                                       >  1:553926/1‑151 (MQ=255)
                                                                                           gTATCCCAACCGTTATCACGTAGTCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGGGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCtt                                                       <  2:328309/151‑1 (MQ=255)
                                                                                            tATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCtt                                                       <  2:774048/150‑1 (MQ=255)
                                                                                             aTCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCttt                                                      >  1:521871/1‑150 (MQ=255)
                                                                                             aTCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCttt                                                      <  1:842771/150‑1 (MQ=255)
                                                                                             aTCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTa                                                     <  1:1232206/151‑1 (MQ=255)
                                                                                               cccAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTAtg                                                   <  2:666207/151‑2 (MQ=255)
                                                                                               cccAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTAt                                                    <  1:1147823/150‑1 (MQ=255)
                                                                                               cccAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTAt                                                    <  1:420461/150‑1 (MQ=255)
                                                                                                  aaCCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCg                                                 >  1:146107/1‑150 (MQ=255)
                                                                                                   aCCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGAt                                               <  1:931118/151‑1 (MQ=255)
                                                                                                    ccGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGAtt                                              >  1:740417/1‑151 (MQ=255)
                                                                                                      gTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTa                                            <  2:1150026/151‑1 (MQ=255)
                                                                                                      gTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTa                                            <  1:1025571/151‑1 (MQ=255)
                                                                                                         aTCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACt                                          <  1:762727/150‑1 (MQ=255)
                                                                                                           cgcgTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGc                                        >  2:1212848/1‑150 (MQ=255)
                                                                                                           cgcgTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCa                                       <  2:239503/151‑1 (MQ=255)
                                                                                                              gTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCaaaa                                    <  2:281999/151‑1 (MQ=255)
                                                                                                              gTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCaa                                      >  1:698034/1‑149 (MQ=255)
                                                                                                               tAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCaaaa                                    >  2:973503/1‑150 (MQ=255)
                                                                                                                aaTCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATg                                  >  2:1070898/1‑151 (MQ=255)
                                                                                                                    aaaGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGc                              <  2:854546/151‑1 (MQ=255)
                                                                                                                    aaaGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGc                              >  1:896012/1‑151 (MQ=255)
                                                                                                                    aaaGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGc                              >  1:349027/1‑151 (MQ=255)
                                                                                                                         agTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCaa                          >  1:603503/1‑150 (MQ=255)
                                                                                                                          gtggtgAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAag                        <  1:130060/151‑1 (MQ=255)
                                                                                                                           tggtgAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGTAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCaa                          >  1:982022/1‑148 (MQ=255)
                                                                                                                           tggtgAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCaa                          >  2:551684/1‑148 (MQ=255)
                                                                                                                            ggtgAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAaga                       <  1:1017685/150‑1 (MQ=255)
                                                                                                                             gtgAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGagc                     <  2:952390/151‑1 (MQ=255)
                                                                                                                               gAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGagc                     >  2:29048/1‑149 (MQ=255)
                                                                                                                               gAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGagc                     >  1:288580/1‑149 (MQ=255)
                                                                                                                                 aaGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCg                 <  1:1082930/151‑1 (MQ=255)
                                                                                                                                       cgcgcgGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTc            >  1:321225/1‑150 (MQ=255)
                                                                                                                                         cgcgGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACt         <  1:963396/151‑1 (MQ=255)
                                                                                                                                          acgGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTc        >  1:870649/2‑151 (MQ=255)
                                                                                                                                             gCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGAt     >  1:974918/1‑151 (MQ=255)
                                                                                                                                             gCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGAt     <  1:974877/151‑1 (MQ=255)
                                                                                                                                                aCCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCAt  >  2:1045499/1‑151 (MQ=255)
                                                                                                                                                aCCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCAt  >  1:1192980/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCAT  >  NC_000913/3822878‑3823172

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: