Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,197,853 T→C 80.2% N129D (AAT→GAT)  ymfE ← e14 prophage; putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,197,8530TC80.2% 248.4 / 57.5 122N129D (AAT→GAT) ymfEe14 prophage; putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (13/11);  new base C (37/60);  total (51/71)
Fisher's exact test for biased strand distribution p-value = 1.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.83e-02

Only 100 of 123 total aligned reads displayed.

TCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAGATA  >  NC_000913/1197754‑1197949
                                                                                                   |                                                                                                
tCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGAtt                                                                                                  <  2:1525573/100‑1 (MQ=255)
  gAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGc                                                                                                >  1:2248641/1‑100 (MQ=255)
  gAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGc                                                                                                <  2:2246167/100‑1 (MQ=255)
  gAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAAAGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGc                                                                                                >  1:883295/1‑100 (MQ=255)
  gAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAATTCAATATAATATTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGc                                                                                                >  1:848035/1‑100 (MQ=255)
   aaaaTTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCt                                                                                               <  1:396743/100‑1 (MQ=255)
     aaTTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCa                                                                                             <  1:947705/100‑1 (MQ=255)
      aTTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCAt                                                                                            >  2:2140967/1‑100 (MQ=255)
       ttCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCAtt                                                                                           <  1:2546194/100‑1 (MQ=255)
       ttCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCAtt                                                                                           <  2:2628343/100‑1 (MQ=255)
        tCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTa                                                                                          <  1:521562/100‑1 (MQ=255)
         cAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTa                                                                                          >  2:595982/1‑99 (MQ=255)
          aaaaaaTATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACa                                                                                        >  1:2389320/1‑100 (MQ=255)
           aaaaaTATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACAt                                                                                       >  2:1284298/1‑100 (MQ=255)
           aaaaaTATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACAt                                                                                       <  2:420047/100‑1 (MQ=255)
            aaaaTATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACAtt                                                                                      <  2:1653809/100‑1 (MQ=255)
               atatCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATtgtg                                                                                   <  1:1948847/100‑1 (MQ=255)
               atatCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATtgtg                                                                                   <  2:2580727/100‑1 (MQ=255)
                tatCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGa                                                                                  <  2:1605802/100‑1 (MQ=255)
                   ccAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCt                                                                               <  2:1915801/100‑1 (MQ=255)
                   ccAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCt                                                                               <  2:2460394/100‑1 (MQ=255)
                    cAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCtt                                                                              >  2:690336/1‑100 (MQ=255)
                       aaaTAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCAt                                                                           <  2:684316/100‑1 (MQ=255)
                             tataCTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCtttat                                                                     <  1:83661/100‑1 (MQ=255)
                             tataCTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCtttat                                                                     <  2:883286/100‑1 (MQ=255)
                             tataCTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCtttat                                                                     >  2:1037023/1‑100 (MQ=255)
                              atgcTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCATGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCtttatt                                                                    <  2:970864/97‑1 (MQ=255)
                              ataCTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCtttatt                                                                    <  1:1184598/100‑1 (MQ=255)
                               taCTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCtttattt                                                                   <  1:429653/100‑1 (MQ=255)
                               taCTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGAATAATATAGATAGCTCATTACATTGTGACCTTCATCtttattt                                                                   >  1:713009/1‑100 (MQ=255)
                                aCTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCtttattta                                                                  >  2:227160/1‑100 (MQ=255)
                                aCTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCtttattaa                                                                  >  1:1648242/1‑98 (MQ=255)
                                       tCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCt                                                           <  2:553638/100‑1 (MQ=255)
                                        tcAGAAGAGTTCAATATAATCTTTGTCTTCAATTTTTCTTACTTCAGGTTAATATAGTTCGCTCATTACATTGTGAGCTTCTTCTTTATTTAATTGTCTg                                                          <  2:775449/99‑1 (MQ=255)
                                        ccAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTg                                                          <  2:2248614/100‑1 (MQ=255)
                                        ccAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTg                                                          >  2:1752270/1‑100 (MQ=255)
                                        ccAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTg                                                          <  1:2647/100‑1 (MQ=255)
                                         cAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGt                                                         >  1:813991/1‑100 (MQ=255)
                                         cAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGt                                                         >  1:1559878/1‑100 (MQ=255)
                                         cAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTAATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGt                                                         >  1:920443/1‑100 (MQ=255)
                                         cAGAAGAGTTCAATACAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTg                                                          <  2:1548101/99‑1 (MQ=255)
                                          agaagaGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGtt                                                        <  1:2492090/100‑1 (MQ=255)
                                          agaagaGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGtt                                                        <  1:784328/100‑1 (MQ=255)
                                          agaagaGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGtt                                                        >  1:837716/1‑100 (MQ=255)
                                           gaagaGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTg                                                       <  1:1575371/100‑1 (MQ=255)
                                            aagaGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCGCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTg                                                          <  1:2368032/96‑1 (MQ=255)
                                                 ttCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCa                                                 <  2:2040399/100‑1 (MQ=255)
                                                  tCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAg                                                >  2:1846492/1‑100 (MQ=255)
                                                   cAATATAATGTTTGTCTTTAATTTTTCTTACTTCAGTGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGc                                               <  1:2505360/100‑1 (MQ=255)
                                                   cAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGc                                               <  1:459317/100‑1 (MQ=255)
                                                   cAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGc                                               <  2:771890/100‑1 (MQ=255)
                                                   cAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGc                                               <  2:632347/100‑1 (MQ=255)
                                                     atatAATGTTTGTCTTCAATTTTTTTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGc                                               <  1:481102/98‑1 (MQ=255)
                                                      tataATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGctct                                            <  1:1305816/100‑1 (MQ=255)
                                                       ataATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCtctc                                           >  2:2436256/1‑100 (MQ=255)
                                                        taATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTcc                                          <  1:619235/100‑1 (MQ=255)
                                                          aTGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGt                                        >  1:2007848/1‑100 (MQ=255)
                                                              ttGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATa                                    <  1:2600645/100‑1 (MQ=255)
                                                                  cTTCAATTTTTTTTTTTTTTTGTTAATATAGAACGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATaa                                   >  2:531064/1‑97 (MQ=255)
                                                                  cTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACgg                                >  2:2524327/1‑100 (MQ=255)
                                                                  cTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACgg                                >  2:1157160/1‑100 (MQ=255)
                                                                   ttCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGt                               >  2:142109/1‑100 (MQ=255)
                                                                    tCAATTTTTCTTACTTCAGGGTAATAGAGATCGCTCATTACATTGTGATCTTCATCTTTATTTAATTTTCTGTTGACTTCAGCTCTCCGTGATAACGGtt                              >  2:2118264/1‑100 (MQ=255)
                                                                     cAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGttt                             <  2:955449/100‑1 (MQ=255)
                                                                     cAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGttt                             <  1:1037033/100‑1 (MQ=255)
                                                                     cAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGttt                             <  2:1242666/100‑1 (MQ=255)
                                                                          tttCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATaa                        <  2:2307562/100‑1 (MQ=255)
                                                                             cTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAAtaa                     >  2:1612968/1‑98 (MQ=255)
                                                                             cTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTa                     <  1:110858/100‑1 (MQ=255)
                                                                             cTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTa                     >  1:1956561/1‑100 (MQ=255)
                                                                             cTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTa                     >  1:2340098/1‑100 (MQ=255)
                                                                             cTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTa                     >  1:2468590/1‑100 (MQ=255)
                                                                              ttACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTTACTCCAGCTCTCCGTGATAACGGTTTTATAATTa                     >  1:1334800/1‑99 (MQ=255)
                                                                               tACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGa                   <  1:2317872/100‑1 (MQ=255)
                                                                               tACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGa                   >  2:305292/1‑100 (MQ=255)
                                                                                aCTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGAt                  >  2:2415658/1‑100 (MQ=255)
                                                                                 cTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATg                 <  1:1647714/100‑1 (MQ=255)
                                                                                 cTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATg                 >  2:1003795/1‑100 (MQ=255)
                                                                                  ttCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGc                <  1:1696511/100‑1 (MQ=255)
                                                                                   tCAGGGTAATATAGTTCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGAGATAACGGTTTTATAATTAGATGCt               <  1:2003969/100‑1 (MQ=255)
                                                                                   tCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCt               <  2:2389290/100‑1 (MQ=255)
                                                                                   tCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCt               >  1:1690400/1‑100 (MQ=255)
                                                                                   tCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCt               >  1:2416812/1‑100 (MQ=255)
                                                                                        gTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATcc          <  1:1406792/100‑1 (MQ=255)
                                                                                        gTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATcc          <  1:1414516/100‑1 (MQ=255)
                                                                                        gTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATcc          <  1:2456139/100‑1 (MQ=255)
                                                                                         tAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATccc         <  1:2311215/100‑1 (MQ=255)
                                                                                         tAATATAGATCGCTCATTACATTGTTAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATccc         <  2:1273234/100‑1 (MQ=255)
                                                                                         tAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCCGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATccc         >  2:843463/1‑100 (MQ=255)
                                                                                         tAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATccc         <  2:765049/100‑1 (MQ=255)
                                                                                          aaTATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTa             <  1:2056405/95‑1 (MQ=255)
                                                                                          aaTATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCa        <  1:2415688/100‑1 (MQ=255)
                                                                                             ataGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAg     <  2:929393/100‑1 (MQ=255)
                                                                                             ataGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAg     >  2:267722/1‑100 (MQ=255)
                                                                                              taGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAGa    <  2:272481/100‑1 (MQ=255)
                                                                                              taGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAGa    >  2:679114/1‑100 (MQ=255)
                                                                                              taGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAGa    <  1:979317/100‑1 (MQ=255)
                                                                                                gATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAGata  >  1:265835/1‑100 (MQ=255)
                                                                                                gATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAGata  >  1:1633726/1‑100 (MQ=255)
                                                                                                gATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAGata  >  2:1464460/1‑100 (MQ=255)
                                                                                                   |                                                                                                
TCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAGATA  >  NC_000913/1197754‑1197949

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