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% 275.7 / NA 69R209H (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 (26/43);  total (26/43)

GTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATC  >  NC_000913/3822878‑3823170
                                                                                                                                                   |                                                                                                                                                 
gTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCTAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATaa                                                                                                                                                <  1:85890/151‑1 (MQ=255)
  gAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGTGCGGCAACCATAAag                                                                                                                                              <  2:819364/151‑1 (MQ=255)
        ctcGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAgat                                                                                                                                          <  1:531595/149‑1 (MQ=255)
          cGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAgatg                                                                                                                                         <  1:19620/148‑1 (MQ=255)
            aaaTTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCa                                                                                                                                    <  2:895749/151‑1 (MQ=255)
                 tataGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCagaag                                                                                                                                <  1:501562/150‑1 (MQ=255)
                 tataGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCaga                                                                                                                                  >  1:436579/1‑148 (MQ=255)
                    aGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGAtt                                                                                                                             <  2:873880/150‑1 (MQ=255)
                      cccGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGAttctt                                                                                                                          <  1:693611/151‑1 (MQ=255)
                         gCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGAttctttct                                                                                                                       <  1:500403/151‑1 (MQ=255)
                            ggCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGaa                                                                                                                    <  1:1103897/151‑1 (MQ=255)
                            ggCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGaa                                                                                                                    <  2:181420/151‑1 (MQ=255)
                            ggCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGaa                                                                                                                    <  1:613853/151‑1 (MQ=255)
                            ggCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGa                                                                                                                     <  2:1035830/150‑1 (MQ=255)
                             gcgcACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGaa                                                                                                                    >  2:1096714/1‑150 (MQ=255)
                               gcACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATc                                                                                                                 >  1:506526/1‑151 (MQ=255)
                                 aCCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGa                                                                                                               <  2:1035765/151‑1 (MQ=255)
                                  ccGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGa                                                                                                               >  1:903210/1‑150 (MQ=255)
                                   cGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAg                                                                                                             <  2:554677/151‑1 (MQ=255)
                                     tttAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAggg                                                                                                           >  1:1112485/1‑151 (MQ=255)
                                       tAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCg                                                                                                         >  2:883438/1‑151 (MQ=255)
                                       tAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCg                                                                                                         >  1:44667/1‑151 (MQ=255)
                                         ggTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGtt                                                                                                       <  2:328393/151‑1 (MQ=255)
                                            aTCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGc                                                                                                    >  1:1133272/1‑151 (MQ=255)
                                              ccaccaCATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAg                                                                                                  <  1:699873/151‑1 (MQ=255)
                                               caccacATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAgg                                                                                                 <  1:375183/151‑1 (MQ=255)
                                                 ccacATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGaa                                                                                               <  2:225114/151‑1 (MQ=255)
                                                  cacaTTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAg                                                                                              <  2:190289/151‑1 (MQ=255)
                                                     aTTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCgg                                                                                           <  2:254416/151‑1 (MQ=255)
                                                          aaCCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATa                                                                                      <  1:668922/151‑1 (MQ=255)
                                                              gAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGt                                                                                  >  2:751274/1‑151 (MQ=255)
                                                               aaCTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTg                                                                                 <  2:897423/151‑1 (MQ=255)
                                                                    gAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCcgcg                                                                            <  2:253482/151‑1 (MQ=255)
                                                                      agagCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTc                                                                          <  1:873037/151‑1 (MQ=255)
                                                                        agCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTc                                                                          >  1:95128/1‑149 (MQ=255)
                                                                           tGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTgg                                                                     >  1:802716/1‑151 (MQ=255)
                                                                                tttGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTcgc                                                                 <  1:173699/150‑1 (MQ=255)
                                                                                  tGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTcgcg                                                                >  2:383599/1‑149 (MQ=255)
                                                                                  tGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGACGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTcgcg                                                                >  1:516614/1‑149 (MQ=255)
                                                                                    aGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGaa                                                            <  2:673731/151‑1 (MQ=255)
                                                                                       cgcTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCa                                                         <  1:1042804/151‑1 (MQ=255)
                                                                                         cTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATc                                                       >  2:547972/1‑151 (MQ=255)
                                                                                            tATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCttt                                                    >  1:107291/1‑151 (MQ=255)
                                                                                             aTCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTa                                                   <  1:986510/151‑1 (MQ=255)
                                                                                               cccAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTAtt                                                 <  2:199639/151‑1 (MQ=255)
                                                                                                 cAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCg                                               <  1:648767/151‑1 (MQ=255)
                                                                                                 cAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCg                                               <  2:319543/151‑1 (MQ=255)
                                                                                                  aaCCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGa                                              <  1:449670/151‑1 (MQ=255)
                                                                                                      gTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGAttt                                           >  1:422803/1‑150 (MQ=255)
                                                                                                           cgcgTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGc                                      >  2:220759/1‑150 (MQ=255)
                                                                                                              gTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCaaaa                                  <  2:310369/151‑1 (MQ=255)
                                                                                                              gTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCaaaa                                  <  2:841202/151‑1 (MQ=255)
                                                                                                              gTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCaaa                                   >  2:253260/1‑150 (MQ=255)
                                                                                                                 aTCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATgg                               <  1:433805/151‑1 (MQ=255)
                                                                                                                 aTCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATgg                               >  2:666469/1‑151 (MQ=255)
                                                                                                                 aTCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATgg                               >  2:279423/1‑151 (MQ=255)
                                                                                                                    aaaGAAGTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGc                            >  1:220180/1‑151 (MQ=255)
                                                                                                                       gaagTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCa                         <  1:1040049/151‑1 (MQ=255)
                                                                                                                        aagTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCa                         >  1:712609/1‑150 (MQ=255)
                                                                                                                         agTGGTGAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCaaa                       <  1:788991/151‑1 (MQ=255)
                                                                                                                               gAAAGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGagcag                 <  2:604339/151‑1 (MQ=255)
                                                                                                                                 aaGCCGCGCGCGGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGagcag                 <  1:227021/149‑1 (MQ=255)
                                                                                                                                     cgcgcgcgGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGtttt           >  1:1144745/1‑151 (MQ=255)
                                                                                                                                       cgcgcgGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTc          >  1:1123301/1‑150 (MQ=255)
                                                                                                                                       cgcgcgGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCa         <  1:785642/151‑1 (MQ=255)
                                                                                                                                         cgcgGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACt       <  2:458993/151‑1 (MQ=255)
                                                                                                                                         cgcgGCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACt       >  1:108371/1‑151 (MQ=255)
                                                                                                                                             gCAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGGTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGa    <  1:941450/150‑1 (MQ=255)
                                                                                                                                              cAACCATAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATc  >  1:367038/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATC  >  NC_000913/3822878‑3823170

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