Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,823,987 G→T G530C (GGT→TGT)  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,9870GT100.0% 261.2 / NA 68G530C (GGT→TGT) spoTbifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (35/33);  total (35/33)

CTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTGGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGCGTGCTGATCACCTTTGCGAAATGCTGCCGCCCTATTCCTG  >  NC_000913/3823844‑3824134
                                                                                                                                               |                                                                                                                                                   
cTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGc                                                                                                                                              >  1:1036195/1‑151 (MQ=255)
     aGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGa                                                                                                                                         <  1:1051731/151‑1 (MQ=255)
     aGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGa                                                                                                                                         <  2:11755/151‑1 (MQ=255)
       ccTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGc                                                                                                                                       <  2:939413/151‑1 (MQ=255)
       ccTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGc                                                                                                                                       <  1:126850/151‑1 (MQ=255)
          gggCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCgtg                                                                                                                                    <  2:596862/151‑1 (MQ=255)
           ggCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCgtgg                                                                                                                                   <  2:496926/151‑1 (MQ=255)
              cgtcgtCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGtggtgg                                                                                                                                <  2:299507/151‑1 (MQ=255)
               gtcgtcTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTggtggt                                                                                                                               >  2:269518/1‑151 (MQ=255)
               gtcgtcTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTggtggt                                                                                                                               <  2:274482/151‑1 (MQ=255)
                           aCCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCt                                                                                                                   >  2:790308/1‑151 (MQ=255)
                                   ttGGGTGGTAGCCGTAAGCTGAATGAAAGCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACa                                                                                                             <  1:74627/149‑1 (MQ=255)
                                    tGGGCGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCCACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATgg                                                                                                          >  1:524499/1‑151 (MQ=255)
                                            aGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGcc                                                                                                   <  1:256218/150‑1 (MQ=255)
                                            aGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCt                                                                                                  <  2:688644/151‑1 (MQ=255)
                                             gCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTc                                                                                                 <  2:1290020/151‑1 (MQ=255)
                                              ccGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTcc                                                                                                >  2:191073/1‑151 (MQ=255)
                                                   aGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTcc                                                                                           <  2:644274/151‑1 (MQ=255)
                                                    gCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCa                                                                                          <  1:444017/151‑1 (MQ=255)
                                                      tgaatgaaATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCAcc                                                                                        >  1:382554/1‑151 (MQ=255)
                                                      tgaatgaaATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCAcc                                                                                        >  1:148629/1‑151 (MQ=255)
                                                      tgaatgaaATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCAcc                                                                                        >  2:509200/1‑151 (MQ=255)
                                                      tgaatgaaATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCAcc                                                                                        >  1:678641/1‑151 (MQ=255)
                                                       gaatgaaATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCg                                                                                       <  2:607826/151‑1 (MQ=255)
                                                           gaaATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCa                                                                                    <  2:95762/150‑1 (MQ=255)
                                                            aaaTCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAAc                                                                                  >  2:350966/1‑151 (MQ=255)
                                                              aTCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAAccc                                                                                <  1:1229627/151‑1 (MQ=255)
                                                                 ccGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCaaa                                                                             <  1:85745/151‑1 (MQ=255)
                                                                 ccGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCaaa                                                                             <  2:711415/151‑1 (MQ=255)
                                                                  cGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCa                                                                               <  2:1076648/148‑1 (MQ=255)
                                                                    cAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGcc                                                                          >  2:570530/1‑151 (MQ=255)
                                                                     aGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGcc                                                                          >  1:373391/1‑150 (MQ=255)
                                                                       gAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACg                                                                       >  2:1171868/1‑151 (MQ=255)
                                                                           atatTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGAc                                                                    >  1:166884/1‑150 (MQ=255)
                                                                           atatTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGAc                                                                    >  1:1256934/1‑150 (MQ=255)
                                                                             atTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATc                                                                 >  2:1206814/1‑151 (MQ=255)
                                                                                    gcgAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCATTGAGCGTGGTGGACGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCa                                                           >  1:800718/1‑150 (MQ=255)
                                                                                      gAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTa                                                        <  2:388158/151‑1 (MQ=255)
                                                                                          tGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAgg                                                    >  1:797583/1‑151 (MQ=255)
                                                                                           ggATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAg                                                     <  1:431954/149‑1 (MQ=255)
                                                                                            gATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTg                                                  >  1:628764/1‑151 (MQ=255)
                                                                                            gATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTg                                                  <  1:431867/151‑1 (MQ=255)
                                                                                                 cATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGAt                                             >  2:47603/1‑151 (MQ=255)
                                                                                                        cTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGCGTGc                                      >  1:760352/1‑151 (MQ=255)
                                                                                                          ggCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGCGTGCTg                                    <  1:453056/151‑1 (MQ=255)
                                                                                                           gCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGCGTGCTg                                    >  2:621848/1‑150 (MQ=255)
                                                                                                            cAACGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGCGTGCTGAt                                  >  2:1110273/1‑151 (MQ=255)
                                                                                                             aaCGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAAGGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGCGTGCTGATc                                 <  1:245564/151‑1 (MQ=255)
                                                                                                              aCGCTTGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGCGTGCTGATc                                 >  2:359881/1‑150 (MQ=255)
                                                                                                                   tGACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGCGTGCTGATCACCtt                            <  2:893425/150‑1 (MQ=255)
                                                                                                                    gACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGCGTGCTGATCACCttt                           <  2:727261/150‑1 (MQ=255)
                                                                                                                    gACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGCGTGCTGATCACCttt                           <  2:775116/150‑1 (MQ=255)
                                                                                                                    gACGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGCGTGCTGATCACCttt                           <  1:1046255/150‑1 (MQ=255)
                                                                                                                     aCGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGCGTGCTGATCACCTTTGc                         <  2:1153305/151‑1 (MQ=255)
                                                                                                                      cGATCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGCGTGCTGATCACCTTTGCg                        >  2:932550/1‑151 (MQ=255)
                                                                                                                        aTCTGCTGGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGCGTGCTGATCACCTTTGCGa                       >  1:440303/1‑150 (MQ=255)
                                                                                                                          ctgctgGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGCGTGCTGATCACCTTTGCGAAAt                    >  2:127429/1‑151 (MQ=255)
                                                                                                                            gctgGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGCGTGCTGATCACCTTTGCGaaa                     <  2:1224327/148‑1 (MQ=255)
                                                                                                                             ctgGCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGCGTGCTGATCACCTTTGCGAAAtgct                 >  1:1082095/1‑151 (MQ=255)
                                                                                                                                gCAGAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGCGTGCTGATCACCTTTGCGAAATGCTgcc              >  1:466006/1‑151 (MQ=255)
                                                                                                                                   gAAATCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGCGTGCTGATCACCTTTGCGAAATGCTgccgcc           <  2:163713/151‑1 (MQ=255)
                                                                                                                                     aaTCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGCGTGCTGATCACCTTTGCGAAATGCTGCCGccc          >  1:539567/1‑150 (MQ=255)
                                                                                                                                      aTCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGCGTGCTGATCACCTTTGCGAAATGCTGCCGCCCt         >  2:1287529/1‑150 (MQ=255)
                                                                                                                                      aTCGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGCGTGCTGATCACCTTTGCGAAATGCTGCCGCCCt         >  1:1281724/1‑150 (MQ=255)
                                                                                                                                        cGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGCGTGCTGATCACCTTTGCGAAATGCTGCCGCCCTAtt      >  2:16895/1‑151 (MQ=255)
                                                                                                                                        cGGACTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGCGTGCTGATCACCTTTGCGAAATGCTGCCGCCCTAt       <  2:1036259/150‑1 (MQ=255)
                                                                                                                                           aCTTTGTAACGCAATGAGCGTGGTGGCCGCGAAAAATCAGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGCGTGCTGATCACCTTTGCGAAATGCTGCCGCCCTATTc     >  1:13612/1‑149 (MQ=255)
                                                                                                                                            cTTTGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGCGTGCTGATCACCTTTGCGAAATGCTGCCGCCCTATTCCTg  >  1:154095/1‑151 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
CTGTAAGCCTGGGCCGTCGTCTGCTCAACCATGCTTTGGGTGGTAGCCGTAAGCTGAATGAAATCCCGCAGGAAAATATTCAGCGCGAGCTGGATCGCATGAAGCTGGCAACGCTTGACGATCTGCTGGCAGAAATCGGACTTGGTAACGCAATGAGCGTGGTGGTCGCGAAAAATCTGCAACATGGGGACGCCTCCATTCCACCGGCAACCCAAAGCCACGGACATCTGCCCATTAAAGGTGCCGATGGCGTGCTGATCACCTTTGCGAAATGCTGCCGCCCTATTCCTG  >  NC_000913/3823844‑3824134

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