Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,912,426 A→G F409S (TTT→TCT)  relA ← (p)ppGpp synthetase I/GTP pyrophosphokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,912,4260AG99.2% 427.5 / ‑4.7 122F409S (TTT→TCT) relA(p)ppGpp synthetase I/GTP pyrophosphokinase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (55/66);  total (56/66)
Fisher's exact test for biased strand distribution p-value = 4.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

Only 100 of 127 total aligned reads displayed.

CCCGATGCAGCGGTGTCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACGCAATCAGTTTACGC  >  NC_000913/2912329‑2912521
                                                                                                 |                                                                                               
cccGATGCAGCGGTGTCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTagag                                                                                               <  2:2052825/100‑1 (MQ=255)
cccGATGCAGCGGTGTCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTagag                                                                                               <  1:2010672/100‑1 (MQ=255)
cccGATGCAGCGGTGTCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTagag                                                                                               >  1:1856314/1‑100 (MQ=255)
 ccGATGCAGCGGTGTCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTagagn                                                                                              <  2:1356575/100‑2 (MQ=255)
 ccGATGCAGCGGTGTCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAgaga                                                                                              >  2:2286961/1‑100 (MQ=255)
     tGCAGCGGTGTCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTa                                                                                          <  2:2270051/100‑1 (MQ=255)
     tGCAGCGGTGTCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTa                                                                                          >  1:145201/1‑100 (MQ=255)
     tGCAGCGGTGTCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTa                                                                                          <  1:2104291/100‑1 (MQ=255)
     cgCAGCGGTGTCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTa                                                                                          <  2:2892231/99‑1 (MQ=255)
      gCAGCGGTGTCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTac                                                                                         >  1:1673695/1‑100 (MQ=255)
       cAGCGGTGTCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGCCGTaca                                                                                        >  2:2618312/1‑100 (MQ=255)
       cAGCGGTATCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTaca                                                                                        >  1:1619414/1‑100 (MQ=255)
        aGCGGTGTCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTaca                                                                                        <  2:2876089/99‑1 (MQ=255)
         gCGGTGTCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACAcc                                                                                      <  2:2354258/100‑1 (MQ=255)
          cGGTGTCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACAccc                                                                                     <  2:2483784/100‑1 (MQ=255)
             tgtCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGgt                                                                                  <  1:606116/100‑1 (MQ=255)
              gtCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGgtc                                                                                 >  1:1413659/1‑100 (MQ=255)
              gtCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGgtc                                                                                 >  2:454924/1‑100 (MQ=255)
              gtCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGgtc                                                                                 >  2:454964/1‑100 (MQ=255)
               tCCGACATCACTGTGGATGTGGTAAGCTAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGgtcg                                                                                <  2:2207991/100‑1 (MQ=255)
                 cGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGgtcgtc                                                                              >  2:2069221/1‑100 (MQ=255)
                 cGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGgtcgtc                                                                              >  1:1945833/1‑100 (MQ=255)
                 cGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGgtcgtc                                                                              <  1:2487651/100‑1 (MQ=255)
                  gACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCa                                                                             <  2:2586298/100‑1 (MQ=255)
                  gACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCa                                                                             <  2:2587170/100‑1 (MQ=255)
                   aCATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCaa                                                                            >  2:705172/1‑100 (MQ=255)
                   aCATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCaa                                                                            <  2:291115/100‑1 (MQ=255)
                     agcACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAg                                                                          <  1:1512161/98‑1 (MQ=255)
                      tCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGa                                                                         >  1:1653590/1‑100 (MQ=255)
                       cACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGAc                                                                        <  1:2127840/100‑1 (MQ=255)
                        aCTGTGGCTGTCGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTACGGCGTAGAGACGTACACCCGGTCGTCAAAGAcc                                                                       >  2:1620829/1‑100 (MQ=255)
                        aCTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGAcc                                                                       >  1:1539492/1‑100 (MQ=255)
                        aCTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGAcc                                                                       <  1:331269/100‑1 (MQ=255)
                          tgtgGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTg                                                                     >  1:1180984/1‑100 (MQ=255)
                          tgtgGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTg                                                                     <  2:2151654/100‑1 (MQ=255)
                              gATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTa                                                                 >  1:750596/1‑100 (MQ=255)
                               aTGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTAc                                                                <  2:1844284/100‑1 (MQ=255)
                                 gtgGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGt                                                              <  2:1945811/100‑1 (MQ=255)
                                  tgGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTa                                                             >  1:2946290/1‑100 (MQ=255)
                                  tgGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTa                                                             >  2:1334184/1‑100 (MQ=255)
                                     tAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACtt                                                          >  2:2293130/1‑100 (MQ=255)
                                       aGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTtcg                                                        >  1:42398/1‑100 (MQ=255)
                                       aGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTtcg                                                        >  1:92498/1‑100 (MQ=255)
                                       aGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTtcg                                                        <  2:47305/100‑1 (MQ=255)
                                         cGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTtcgtc                                                      >  1:1582784/1‑100 (MQ=255)
                                         cGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTtcgtc                                                      >  1:2614332/1‑100 (MQ=255)
                                          gAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACCCCCGGTCGTCAAAGACCTGACTACGTACTtcgtcg                                                     >  1:2448880/1‑100 (MQ=255)
                                          gAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTtcgtcg                                                     <  1:988256/100‑1 (MQ=255)
                                           aaGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGa                                                    <  2:2696077/100‑1 (MQ=255)
                                             gTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGc                                                  <  2:2588074/100‑1 (MQ=255)
                                              tCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCa                                                 >  1:783025/1‑100 (MQ=255)
                                              tCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCa                                                 <  1:1940936/100‑1 (MQ=255)
                                                    ggcgTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGc                                           >  2:2213884/1‑100 (MQ=255)
                                                    ggcgTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGc                                           <  2:521133/100‑1 (MQ=255)
                                                     gcgTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGcc                                          <  2:2668212/100‑1 (MQ=255)
                                                     gcgTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGcc                                          >  2:284168/1‑100 (MQ=255)
                                                      cgTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCg                                         <  2:1053676/100‑1 (MQ=255)
                                                      cgTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCg                                         >  1:1141048/1‑100 (MQ=255)
                                                       gTTGATCCCGCAGGCAAATCAACGACATCACCTTTTGGCGTAGAGACGTACACTCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCg                                         >  2:2450391/1‑99 (MQ=255)
                                                       gTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCgg                                        <  2:885284/100‑1 (MQ=255)
                                                          gATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAAt                                     >  2:991516/1‑100 (MQ=255)
                                                           aTCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATc                                    <  2:1309026/100‑1 (MQ=255)
                                                           aTCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATc                                    <  1:942135/100‑1 (MQ=255)
                                                             cccGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAg                                  >  1:2469748/1‑100 (MQ=255)
                                                             cccGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAg                                  >  2:736997/1‑100 (MQ=255)
                                                              ccGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAg                                  <  1:2470983/99‑1 (MQ=255)
                                                              ccGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGc                                 >  2:2618073/1‑100 (MQ=255)
                                                               cGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGcc                                <  2:1752333/100‑1 (MQ=255)
                                                               cGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGcc                                <  1:318046/100‑1 (MQ=255)
                                                               cGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGcc                                <  1:234501/100‑1 (MQ=255)
                                                               cGCAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGcc                                <  1:1354044/100‑1 (MQ=255)
                                                                 cAGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCAt                              <  1:3049194/100‑1 (MQ=255)
                                                                  aGGCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCAtc                             <  2:1327761/100‑1 (MQ=255)
                                                                    gCAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCAtctc                           <  2:1413647/100‑1 (MQ=255)
                                                                     cAAATCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATctct                          >  1:2224704/1‑100 (MQ=255)
                                                                      aaaTCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCtt                         >  1:2696790/1‑100 (MQ=255)
                                                                        aTCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTcc                       <  2:2188573/100‑1 (MQ=255)
                                                                        aTCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTcc                       >  1:2014604/1‑100 (MQ=255)
                                                                         tCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCt                      <  2:704871/100‑1 (MQ=255)
                                                                         tCAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCt                      <  2:2946260/100‑1 (MQ=255)
                                                                          cAACGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTg                     <  2:750591/100‑1 (MQ=255)
                                                                            aCGACCTCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGcc                   <  2:2925450/100‑1 (MQ=255)
                                                                            aCGACATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGcc                   <  2:2490869/100‑1 (MQ=255)
                                                                               aCATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACg                <  1:2723827/100‑1 (MQ=255)
                                                                                cATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACGc               <  1:3064036/100‑1 (MQ=255)
                                                                                cATCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACGc               <  2:874991/100‑1 (MQ=255)
                                                                                 aTCACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACGCa              >  2:2803817/1‑100 (MQ=255)
                                                                                   cACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACGCAAt            >  2:1476993/1‑100 (MQ=255)
                                                                                   cACCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACGCAAt            >  1:1706851/1‑100 (MQ=255)
                                                                                    aCCTTTCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGTCATCTCTTCCTGCCACGCAAt            >  1:1370619/1‑99 (MQ=255)
                                                                                     ccTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACGCAATCa          >  1:1099187/1‑100 (MQ=255)
                                                                                       tttCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACGCAATCAGt        <  1:705175/100‑1 (MQ=255)
                                                                                        ttCGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACGCAATCAGtt       <  2:1746259/100‑1 (MQ=255)
                                                                                          cGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACGCAATCAGttt      <  2:1482269/99‑1 (MQ=255)
                                                                                          cGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACGCAATCAGTTTa     <  1:454926/100‑1 (MQ=255)
                                                                                          cGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACGCAATCAGTTTa     <  1:454966/100‑1 (MQ=255)
                                                                                          cGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACGCAATCAGTTTa     >  1:1070109/1‑100 (MQ=255)
                                                                                          cGGCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACGCAATCAGTTTa     <  1:1334195/100‑1 (MQ=255)
                                                                                           ggCGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACGCAATCAGTTTAc    >  2:2464502/1‑100 (MQ=255)
                                                                                             cGTAGAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACGCAATCAGTTTACGc  >  1:2359556/1‑100 (MQ=255)
                                                                                                 |                                                                                               
CCCGATGCAGCGGTGTCCGACATCACTGTGGATGTGGTAAGCGAAGTCCAGCGGCGTTGATCCCGCAGGCAAATCAACGACATCACCTTTCGGCGTAAAGACGTACACCCGGTCGTCAAAGACCTGACTACGTACTTCGTCGAGCATTTCGCCGGAATCAGCCATCTCTTCCTGCCACGCAATCAGTTTACGC  >  NC_000913/2912329‑2912521

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: