Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 160,736 A→T intergenic (‑132/+46) dksA ← / ← sfsA transcriptional regulator of rRNA transcription; DnaK suppressor protein/sugar fermentation stimulation protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913160,7360AT97.8% 446.4 / ‑4.3 139intergenic (‑132/+46)dksA/sfsAtranscriptional regulator of rRNA transcription; DnaK suppressor protein/sugar fermentation stimulation protein A
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base T (67/68);  total (69/70)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 140 total aligned reads displayed.

TGTTCGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAAT  >  NC_000913/160639‑160832
                                                                                                 |                                                                                                
  ttCGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTg                                                                                              <  2:1761536/100‑1 (MQ=255)
  ttCGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTg                                                                                              <  2:1868440/100‑1 (MQ=255)
    cGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTg                                                                                              <  2:3046205/98‑1 (MQ=255)
      ggggAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCa                                                                                          <  2:1705711/100‑1 (MQ=255)
      ggggAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCa                                                                                          >  2:2585218/1‑100 (MQ=255)
       gggAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTttttttt                                                                                            >  1:644407/1‑94 (MQ=255)
        ggAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCAcg                                                                                        <  2:1785420/100‑1 (MQ=255)
        ggAAAAATGAGGCCGCTATAAATAGCAGATGCTTTGCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCAcg                                                                                        <  2:1592306/100‑1 (MQ=255)
         gAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCAcgc                                                                                       <  1:1995984/100‑1 (MQ=255)
            aaaTGAGGCCGCTATAAATAGCAGATGCTTTTTTTTTTTTTTTTTATCTAAAAGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCAcgc                                                                                       >  1:1780961/1‑97 (MQ=255)
            aaaTGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGAc                                                                                    >  1:2558791/1‑100 (MQ=255)
            aaaTGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGAc                                                                                    <  1:1332264/100‑1 (MQ=255)
             aaTGAGGCCTCTATAAATAGCAGATGCTTTTCTGGATAGCAATTATCTAAACGTAACATTTGACAACTTTTTGAGGAAAAGCGTTTTTGCGCAcgcg                                                                                      <  1:2139412/97‑1 (MQ=255)
             aaTGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGAcc                                                                                   >  2:2905158/1‑100 (MQ=255)
             aaTGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGAca                                                                                   >  2:2247782/1‑99 (MQ=255)
               tGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAg                                                                                 <  1:3334113/100‑1 (MQ=255)
                 aGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGaa                                                                               <  2:1083739/100‑1 (MQ=255)
                   gCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATg                                                                             <  2:1901913/100‑1 (MQ=255)
                     cGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTa                                                                           <  2:694741/100‑1 (MQ=255)
                      gCTATAAATAGCAGATGCTTTTCTGGATAGCAATTATCTAAACGTAACGCTTGACAGCTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGt                                                                            <  1:2900274/98‑1 (MQ=255)
                         ataAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAAtt                                                                       <  1:2763039/100‑1 (MQ=255)
                           aaaTAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTaa                                                                     <  2:1300620/100‑1 (MQ=255)
                           aaaTAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTaa                                                                     >  1:3437720/1‑100 (MQ=255)
                           aaaTAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTaa                                                                     <  1:2955585/100‑1 (MQ=255)
                             aTAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAAcc                                                                   <  2:309520/100‑1 (MQ=255)
                             aTAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAAcc                                                                   <  1:242127/100‑1 (MQ=255)
                             aTAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAAcc                                                                   >  1:827056/1‑100 (MQ=255)
                              tAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCa                                                                  >  2:1820580/1‑100 (MQ=255)
                               aGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAg                                                                 <  1:87219/100‑1 (MQ=255)
                               aGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAg                                                                 >  2:3466113/1‑100 (MQ=255)
                                 cAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAgct                                                               >  1:2117739/1‑98 (MQ=255)
                                 cAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGtt                                                               <  1:1949193/100‑1 (MQ=255)
                                  aGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGtta                                                              <  2:357636/100‑1 (MQ=255)
                                   gATGCTTTTCCGTTTAGCAAAAAGCTAAACGTAACACTTGACAACTGTGTGAAGAAAAGCGTTTTTACGCACGCGACCAGAATGTAAATTAACCAGtt                                                               >  1:211241/1‑98 (MQ=255)
                                   gATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGttan                                                             <  1:96968/100‑2 (MQ=255)
                                        ttttCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTtacttact                                                        <  2:547640/100‑1 (MQ=255)
                                         tttCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTacttactt                                                       <  2:339111/100‑1 (MQ=255)
                                          ttCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACttacttt                                                      >  2:257284/1‑100 (MQ=255)
                                          ttCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACttacttt                                                      <  1:1872584/100‑1 (MQ=255)
                                            ccGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACttactttac                                                    >  1:128451/1‑100 (MQ=255)
                                            ccGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACttactttac                                                    >  2:67744/1‑100 (MQ=255)
                                            ccGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACttactttac                                                    <  2:2864738/100‑1 (MQ=255)
                                              ggttAGTAAGTATCTAAACGTAACGATTGACAACGTTGTTAGGCAAAGAGTTTTTGCGCACACCACCAGAATGTAAATTAACCAGTTACTTACTTtacta                                                  <  1:2810348/97‑1 (MQ=255)
                                              ggATAGTAAGTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACTCGACCAGAATGTAAATTAACCAGTTACTTACTTtacta                                                  <  1:2401741/100‑1 (MQ=255)
                                              ggATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTtacta                                                  <  2:1178198/100‑1 (MQ=255)
                                              ggATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAACTAACCAGTTACTTACTTtacta                                                  <  2:273935/100‑1 (MQ=255)
                                               gATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTtactac                                                 <  2:127523/100‑1 (MQ=255)
                                                 tcgCAATTATCTAGACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACaa                                               >  2:2358963/3‑100 (MQ=255)
                                                 tAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACaa                                               <  1:363971/100‑1 (MQ=255)
                                                  aGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAAt                                              <  2:1924014/100‑1 (MQ=255)
                                                   gCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCTCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACa                                                >  1:781488/1‑97 (MQ=255)
                                                   gCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGTCCAGAATGTAAATTAACCAGTTACTTACTTTACTACAAt                                              <  2:1934637/99‑1 (MQ=255)
                                                   gCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATg                                             >  2:1737289/1‑100 (MQ=255)
                                                     aaTTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTa                                           >  2:443035/1‑100 (MQ=255)
                                                     aaTTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTa                                           <  1:2273218/100‑1 (MQ=255)
                                                         aTCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCg                                       >  1:1890331/1‑100 (MQ=255)
                                                          tCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCgg                                      >  2:108552/1‑100 (MQ=255)
                                                          tCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCAACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAAcc                                        >  1:3144206/1‑98 (MQ=255)
                                                           cTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGc                                     >  1:921198/1‑100 (MQ=255)
                                                           cTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGc                                     <  1:6072/100‑1 (MQ=255)
                                                           cTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGc                                     >  2:1372956/1‑100 (MQ=255)
                                                            tAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCa                                    >  2:47931/1‑100 (MQ=255)
                                                            tAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCa                                    >  1:3162684/1‑100 (MQ=255)
                                                              aaCGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGt                                  >  1:3088854/1‑100 (MQ=255)
                                                              aaCGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACCTACTTTACTACAATGTAACCGGCAGt                                  <  2:2117098/100‑1 (MQ=255)
                                                                cGTAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGa                                >  2:795880/1‑100 (MQ=255)
                                                                  tAACACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGAtt                              >  1:1718829/1‑100 (MQ=255)
                                                                   aaCACTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGAttt                             >  2:3439504/1‑100 (MQ=255)
                                                                    acacTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGAtttt                            <  1:945297/100‑1 (MQ=255)
                                                                    aacattGACAACAATGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTTAATTAACTAGTTACGTACTTTACTACAATGTAACCGGCAGTGAtt                              >  1:1506701/5‑98 (MQ=255)
                                                                     cacTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGAttttt                           >  1:671006/1‑100 (MQ=255)
                                                                     cacTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGAttttt                           <  1:3228926/100‑1 (MQ=255)
                                                                     cacTTGACAACTGTGTGAGGAAAAGCGTGTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTATTTACTTTACTACAATGTAACCGGCAGTGAttttt                           <  1:3181538/100‑1 (MQ=255)
                                                                      acTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGAtttttt                          <  2:1357869/100‑1 (MQ=255)
                                                                       cTTGACAACTGTGTGAGGAAAAGCGTTTTTTTTTTTTTTACTAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTa                         >  1:3239884/1‑100 (MQ=255)
                                                                       cTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTa                         <  1:1107260/100‑1 (MQ=255)
                                                                       cTTGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTa                         >  1:1966168/1‑100 (MQ=255)
                                                                         tGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAag                       <  2:2265640/100‑1 (MQ=255)
                                                                         tGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAag                       <  1:1534501/100‑1 (MQ=255)
                                                                         tGACAACTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAag                       <  1:1673413/100‑1 (MQ=255)
                                                                               cTGTGTGAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCAt                 <  1:2828346/100‑1 (MQ=255)
                                                                                  tgtgAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGcc              <  1:1372936/100‑1 (MQ=255)
                                                                                   gtgAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCGGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCt             >  2:2820383/1‑100 (MQ=255)
                                                                                   gtgAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAagatcggaa                >  1:2570378/1‑91 (MQ=255)
                                                                                   gtgAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCt             >  2:2124704/1‑100 (MQ=255)
                                                                                     gAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACCATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCt             >  1:1997139/1‑98 (MQ=255)
                                                                                     gAGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAAACAGTTACTTACTTTACTACAATGTAACGGGCAGTGATTCTTTGAGAGCCATGCCt             >  1:354463/1‑98 (MQ=255)
                                                                                      aGGAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTca          <  2:1009854/100‑1 (MQ=255)
                                                                                       ggAAAAGCGTTTTTGCGCACGCGACCAGATTGTAAATTAACCAGTTACTTACTTTACTACAATGTATCCGGCAGTGATTTTTTAAGAGCCATGCCTTcag         <  1:2091878/100‑1 (MQ=255)
                                                                                       ggAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTcag         <  2:711964/100‑1 (MQ=255)
                                                                                        gAAAAGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTcagc        <  2:1888801/100‑1 (MQ=255)
                                                                                         aaaaGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTcagca       >  1:1533583/1‑100 (MQ=255)
                                                                                         aaaaGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTcagca       >  1:1030609/1‑100 (MQ=255)
                                                                                          aaaGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTcagcag      <  1:2034582/100‑1 (MQ=255)
                                                                                           aaGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGa     >  1:3366301/1‑100 (MQ=255)
                                                                                            aGCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGaa    >  2:3458415/1‑100 (MQ=255)
                                                                                             gCGTTTTTGTGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTATGATCCATGCATTCAGCAGa     >  1:306068/1‑98 (MQ=255)
                                                                                             gCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGaaa   >  2:3103597/1‑100 (MQ=255)
                                                                                             gCGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGaaa   <  2:2653079/100‑1 (MQ=255)
                                                                                              cGTTTTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAAt  >  1:2986615/1‑100 (MQ=255)
                                                                                                 |                                                                                                
TGTTCGGGGGAAAAATGAGGCCGCTATAAATAGCAGATGCTTTTCCGGATAGCAATTATCTAAACGTAACACTTGACAACTGTGTGAGGAAAAGCGTATTTGCGCACGCGACCAGAATGTAAATTAACCAGTTACTTACTTTACTACAATGTAACCGGCAGTGATTTTTTAAGAGCCATGCCTTCAGCAGAAAT  >  NC_000913/160639‑160832

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