Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,823,249 T→C Y284H (TAC→CAC)  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,2490TC100.0% 251.2 / NA 67Y284H (TAC→CAC) spoTbifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (31/36);  total (31/36)

CATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGA  >  NC_000913/3823114‑3823382
                                                                                                                                       |                                                                                                                                     
cATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGcaca                                                                                                                                    >  2:1297297/1‑139 (MQ=255)
cATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGcaca                                                                                                                                    >  2:917722/1‑139 (MQ=255)
cATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGcaca                                                                                                                                    >  1:1417538/1‑139 (MQ=255)
   cTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGCACAAGc                                                                                                                                 <  2:1059395/139‑1 (MQ=255)
             tttACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCggg                                                                                                                       >  1:425829/1‑139 (MQ=255)
              ttACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGc                                                                                                                      <  2:467698/139‑1 (MQ=255)
                    cAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCCGCGTg                                                                                                                <  1:828347/139‑1 (MQ=255)
                    cAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCCGCGTg                                                                                                                <  2:756416/139‑1 (MQ=255)
                    cAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCCGCGTg                                                                                                                <  2:629470/139‑1 (MQ=255)
                            tGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCCGCGTGAAAGACta                                                                                                        <  2:20127/139‑1 (MQ=255)
                              cTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCCGCGTGAAAGACtata                                                                                                      >  1:1124039/1‑139 (MQ=255)
                                                      aTCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGa                                                                                    <  1:1012574/133‑1 (MQ=255)
                                                      aTCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGa                                                                                    >  2:1012574/1‑133 (MQ=255)
                                                         aTGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATc                                                                           <  1:1297297/139‑1 (MQ=255)
                                                           ggACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGCACAAGc                                                                                                                                 >  2:39814/1‑83 (MQ=255)
                                                           ggACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGCACAAGc                                                                                                                                 <  1:39814/83‑1 (MQ=255)
                                                                 cTACGCTTTCCGCGTGCTCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGc                                                                                                  >  1:716551/1‑108 (MQ=255)
                                                                 cTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGCACAAg                                                                                                                                  <  2:425733/76‑1 (MQ=255)
                                                                 cTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGCACAAg                                                                                                                                  >  1:425733/1‑76 (MQ=255)
                                                                 cTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCcgcg                                                                                                                  >  1:474213/1‑92 (MQ=255)
                                                                 cTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCcgcg                                                                                                                  <  2:474213/92‑1 (MQ=255)
                                                                 cTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGc                                                                                                  <  2:716551/108‑1 (MQ=255)
                                                                 cTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAAc                                                                                  <  1:277225/124‑1 (MQ=255)
                                                                 cTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAAc                                                                                  >  2:277225/1‑124 (MQ=255)
                                                                 cTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTg                                                                   >  1:383836/1‑139 (MQ=255)
                                                                 cTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTg                                                                   >  1:428841/1‑139 (MQ=255)
                                                                 cTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTg                                                                   >  1:318214/1‑139 (MQ=255)
                                                                 cTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTg                                                                   >  1:464916/1‑139 (MQ=255)
                                                                 cTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTg                                                                   >  2:611130/1‑139 (MQ=255)
                                                                 cTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTg                                                                   >  1:748365/1‑139 (MQ=255)
                                                                 cTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTg                                                                   >  1:975781/1‑139 (MQ=255)
                                                                 cTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTg                                                                   >  2:123044/1‑139 (MQ=255)
                                                                        ttCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCt                                                            <  2:1362663/139‑1 (MQ=255)
                                                                         tCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTc                                                           >  2:174997/1‑139 (MQ=255)
                                                                           cgcgTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCcgcg                                                                                                                  <  2:1192209/82‑1 (MQ=255)
                                                                           cgcgTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCcgcg                                                                                                                  >  1:1192209/1‑82 (MQ=255)
                                                                           cgcgTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCga                                                         >  2:893543/1‑139 (MQ=255)
                                                                            gcgTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCgat                                                        >  1:1018502/1‑139 (MQ=255)
                                                                              gTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCgatga                                                      >  1:876686/1‑139 (MQ=255)
                                                                              gTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCgatga                                                      <  1:174997/139‑1 (MQ=255)
                                                                                gATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATc                                                    <  1:123044/139‑1 (MQ=255)
                                                                                   cGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGc                                                 <  1:611130/139‑1 (MQ=255)
                                                                                   cGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGc                                                 >  2:703636/1‑139 (MQ=255)
                                                                                     tCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGccc                                               <  2:975781/139‑1 (MQ=255)
                                                                                        aTGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCa                                            <  2:1045949/139‑1 (MQ=255)
                                                                                          gATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACg                                          <  1:703636/139‑1 (MQ=255)
                                                                                            ttCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCCGCGTGAAAGa                                                                                                           >  1:337166/1‑72 (MQ=255)
                                                                                            ttCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCCGCGTGAAAGa                                                                                                           <  2:337166/72‑1 (MQ=255)
                                                                                                gACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGc                                    <  2:1018502/139‑1 (MQ=255)
                                                                                                    ccTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGt                                <  2:748365/139‑1 (MQ=255)
                                                                                                    ccTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGt                                >  1:97010/1‑139 (MQ=255)
                                                                                                      tGTTATCGCGTGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTg                              <  2:789542/139‑1 (MQ=255)
                                                                                                      tGTTATCGCGTGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTg                              <  1:917722/139‑1 (MQ=255)
                                                                                                           tCGCGTGCTGGGCCAGATGCACAGCCTCCACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTc                         <  2:414216/139‑1 (MQ=255)
                                                                                                              cgTGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTc                         >  2:199546/1‑136 (MQ=255)
                                                                                                              cgTGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTc                         <  1:199546/136‑1 (MQ=255)
                                                                                                               gTGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCCGCGTg                                                                                                                >  1:188484/1‑48 (MQ=255)
                                                                                                               gTGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCCGCGTg                                                                                                                <  2:188484/48‑1 (MQ=255)
                                                                                                                tGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGtgtg                                     <  1:233392/122‑1 (MQ=255)
                                                                                                                tGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGtgtg                                     >  2:233392/1‑122 (MQ=255)
                                                                                                                tGCTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGAt                    <  2:464916/139‑1 (MQ=255)
                                                                                                                  cTGGGCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCCCCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGCCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATcc                  >  2:105815/1‑139 (MQ=255)
                                                                                                                     ggCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAAc                                                                                  <  2:166181/72‑1 (MQ=255)
                                                                                                                     ggCCAGATGCACAGCCTGCACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAAc                                                                                  >  1:166181/1‑72 (MQ=255)
                                                                                                                       ccAGATGCACAGCCTGCACAAGCCGCGTCCGGACCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCGGTCTTTGCACACCTCGATGATCGGCCCGCa                                            >  2:176817/1‑108 (MQ=255)
                                                                                                                       ccAGATGCACAGCCTGCACAAGCCGCGTCCGGACCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCGGTCTTTGCACACCTCGATGATCGGCCCGCa                                            <  1:176817/108‑1 (MQ=255)
                                                                                                                         aGATGCACAGCCTGCACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGa           <  2:318214/139‑1 (MQ=255)
                                                                                                                         aGATGCACAGCCTGCACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGa           <  2:1124039/139‑1 (MQ=255)
                                                                                                                           aTGCACAGCCTGCACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAg         >  1:648790/1‑139 (MQ=255)
                                                                                                                                  gCCTGCACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGa  <  1:455527/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                     
CATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATCGCCATTCCAAAAGCGAACGGCTATCAGTCTTTGCACACCTCGATGATCGGCCCGCACGGTGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGA  >  NC_000913/3823114‑3823382

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