Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,823,106 G→T R236L (CGC→CTC)  spoT → bifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,823,1060GT100.0% 241.0 / NA 63R236L (CGC→CTC) 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 (30/33);  total (30/33)

GAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACA  >  NC_000913/3822961‑3823252
                                                                                                                                                 |                                                                                                                                                  
gAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCgag                                                                                                                                                <  1:859720/150‑1 (MQ=255)
 aGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGaa                                                                                                                                              <  1:861668/151‑1 (MQ=255)
  agcgaTGTATCCCAACCGTTATCGAGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGAGGATGCAGAAGATGCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCgag                                                                                                                                                <  1:185948/147‑1 (MQ=255)
    cgcTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCa                                                                                                                                           <  1:366643/151‑1 (MQ=255)
     gcTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCAt                                                                                                                                          <  1:582358/151‑1 (MQ=255)
     gcTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCAt                                                                                                                                          >  2:267568/1‑151 (MQ=255)
       tGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCt                                                                                                                                        >  2:387522/1‑151 (MQ=255)
        gTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCtt                                                                                                                                       <  2:596782/151‑1 (MQ=255)
        gTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCtt                                                                                                                                       >  2:462520/1‑151 (MQ=255)
          aTCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTa                                                                                                                                     <  1:725894/151‑1 (MQ=255)
            cccAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTAt                                                                                                                                    <  1:57181/150‑1 (MQ=255)
              cAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCg                                                                                                                                 <  1:481859/151‑1 (MQ=255)
                aCCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGAt                                                                                                                               <  1:469220/151‑1 (MQ=255)
                  cGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGAttt                                                                                                                             <  2:600354/151‑1 (MQ=255)
                   gTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGAttt                                                                                                                             >  1:265155/1‑150 (MQ=255)
                       tCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGc                                                                                                                        >  1:107558/1‑151 (MQ=255)
                         gcgTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCaa                                                                                                                      <  2:425823/151‑1 (MQ=255)
                            tAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAAt                                                                                                                   >  1:940708/1‑151 (MQ=255)
                             aaTCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAAt                                                                                                                   <  2:625244/150‑1 (MQ=255)
                             aaTCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAAt                                                                                                                   <  2:937276/150‑1 (MQ=255)
                               tCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGt                                                                                                                >  1:97277/1‑151 (MQ=255)
                                     aagTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCa                                                                                                           >  1:780739/1‑150 (MQ=255)
                                            gAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGagca                                                                                                    <  2:431757/150‑1 (MQ=255)
                                            gAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGagca                                                                                                    <  1:221148/150‑1 (MQ=255)
                                             aaaGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGagcagc                                                                                                  >  2:313251/1‑151 (MQ=255)
                                             aaaGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGagcagc                                                                                                  >  2:790526/1‑151 (MQ=255)
                                                  cgcgcgcgGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGtttt                                                                                             >  2:377229/1‑151 (MQ=255)
                                                  cgcgcgcgGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGtt                                                                                               >  2:861840/1‑149 (MQ=255)
                                                    cgcgcgGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTc                                                                                            >  1:585762/1‑150 (MQ=255)
                                                    cgcgcgGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTc                                                                                            >  2:747549/1‑150 (MQ=255)
                                                      cgcgGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACt                                                                                         >  2:199261/1‑151 (MQ=255)
                                                         ggCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCg                                                                                       >  1:695620/1‑150 (MQ=255)
                                                            aaCCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCa                                                                                   <  2:703098/151‑1 (MQ=255)
                                                             aCCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCAt                                                                                  <  1:606680/151‑1 (MQ=255)
                                                              ccGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATg                                                                                 <  2:71578/151‑1 (MQ=255)
                                                                gTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATgg                                                                                >  2:874420/1‑150 (MQ=255)
                                                                   aaGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACAt                                                                            <  2:761827/151‑1 (MQ=255)
                                                                    agagATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATc                                                                           >  2:163109/1‑151 (MQ=255)
                                                                    agagATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATc                                                                           >  1:223952/1‑151 (MQ=255)
                                                                     gagaTGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCt                                                                          >  2:170182/1‑151 (MQ=255)
                                                                      agaTGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAACGAGCAGCGTTTTCACTCGATCATGGACATCt                                                                          >  1:356495/1‑150 (MQ=255)
                                                                        atgatCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACg                                                                       <  1:940124/151‑1 (MQ=255)
                                                                         tgatCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGc                                                                      >  1:783892/1‑151 (MQ=255)
                                                                             ccAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTc                                                                  >  1:616202/1‑151 (MQ=255)
                                                                                  agaTTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGt                                                             <  2:1015250/151‑1 (MQ=255)
                                                                                   gaTTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTg                                                            <  1:527572/151‑1 (MQ=255)
                                                                                   gaTTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTg                                                            >  2:64121/1‑151 (MQ=255)
                                                                                    aTTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGa                                                           <  2:585884/151‑1 (MQ=255)
                                                                                       ctttctGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCg                                                        >  2:630570/1‑151 (MQ=255)
                                                                                          tctGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTc                                                      <  2:18849/150‑1 (MQ=255)
                                                                                                  cGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCt                                             >  1:563388/1‑151 (MQ=255)
                                                                                                    aaGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGa                                           <  1:462453/151‑1 (MQ=255)
                                                                                                    aaGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGa                                           <  1:559517/151‑1 (MQ=255)
                                                                                                     aGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGac                                          <  2:949845/151‑1 (MQ=255)
                                                                                                          gTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTg                                     >  1:536862/1‑151 (MQ=255)
                                                                                                          gTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTg                                     <  1:393842/151‑1 (MQ=255)
                                                                                                                          aaTACCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGc                     <  1:220054/151‑1 (MQ=255)
                                                                                                                             aCCGTGCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGcc                    >  2:184245/1‑149 (MQ=255)
                                                                                                                                  gCCGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGc              <  1:79944/150‑1 (MQ=255)
                                                                                                                                   ccGCGTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGcac            <  1:546226/151‑1 (MQ=255)
                                                                                                                                     gcgTCAGTGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAg          <  1:267527/151‑1 (MQ=255)
                                                                                                                                            tGGTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTAc   >  1:1016429/1‑151 (MQ=255)
                                                                                                                                              gTCTCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACa  >  2:647280/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAGTGGTGAAAGCCGCGCGCGGCAACCGTAAAGAGATGATCCAGAAGATTCTTTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAGAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACA  >  NC_000913/3822961‑3823252

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