Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,197,853 T→C 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,8530TC98.6% 512.2 / ‑5.0 146N129D (AAT→GAT) ymfEe14 prophage; putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base C (65/79);  total (67/79)
Fisher's exact test for biased strand distribution p-value = 2.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 146 total aligned reads displayed.

GAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAGATA  >  NC_000913/1197756‑1197949
                                                                                                 |                                                                                                
gAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGc                                                                                                >  2:2383180/1‑100 (MQ=255)
gAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGc                                                                                                <  1:464313/100‑1 (MQ=255)
gAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGc                                                                                                <  1:464444/100‑1 (MQ=255)
gAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGc                                                                                                <  1:1465523/100‑1 (MQ=255)
  aaaTTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTc                                                                                              <  2:1519867/100‑1 (MQ=255)
     ttCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCAtt                                                                                           >  1:970898/1‑100 (MQ=255)
     tctaaaaaaTATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCAtt                                                                                           <  1:1919337/97‑1 (MQ=255)
      tCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTa                                                                                          <  1:338747/100‑1 (MQ=255)
       cAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTa                                                                                          >  2:2109722/1‑99 (MQ=255)
       cAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTAc                                                                                         <  1:2222361/100‑1 (MQ=255)
       cAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTAc                                                                                         <  1:1838144/100‑1 (MQ=255)
       cAAAAAATATCCAAAAATAGGATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTAc                                                                                         >  2:1763637/1‑100 (MQ=255)
        aaaaaaTATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTaga                                                                                        >  2:1919323/1‑98 (MQ=255)
        aaaaaaTATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACa                                                                                        >  2:487756/1‑100 (MQ=255)
           aaaTATCCAAAAATAGTATACTTTATTCCAGAAGATTTCAATATAATTTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACAtt                                                                                      <  2:2416581/99‑1 (MQ=255)
           aaaTATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCGTTACATtg                                                                                     <  1:1710006/100‑1 (MQ=255)
            aaTATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATtgt                                                                                    <  1:2067746/100‑1 (MQ=255)
            aaTATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATtgt                                                                                    <  1:618464/100‑1 (MQ=255)
             atatCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATtgtg                                                                                   <  2:2191892/100‑1 (MQ=255)
              tatCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGa                                                                                  <  1:2415181/100‑1 (MQ=255)
                  cAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGc                                                                                <  1:2252103/98‑1 (MQ=255)
                  cAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCtt                                                                              <  2:1042128/100‑1 (MQ=255)
                   aaaaaTAGTATACTTTATTCCAGAAGAGTTCAATATTATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCtt                                                                              <  1:883890/99‑1 (MQ=255)
                     aaaTAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCAt                                                                           <  1:1598199/100‑1 (MQ=255)
                      aaTAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATc                                                                          <  1:1805664/100‑1 (MQ=255)
                       aTAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCt                                                                         <  1:2324379/100‑1 (MQ=255)
                          gTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCttta                                                                      <  2:937603/100‑1 (MQ=255)
                           tataCTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCtttat                                                                     <  2:1938136/100‑1 (MQ=255)
                            ataCTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCtttatt                                                                    <  2:2493267/100‑1 (MQ=255)
                              aCTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAGTATAGATCGCTCATTACATTGTGAGCTTCATCtttattta                                                                  <  1:563751/100‑1 (MQ=255)
                                tttATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAAt                                                                <  2:940093/100‑1 (MQ=255)
                                tttATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAAt                                                                <  2:566212/100‑1 (MQ=255)
                                tttATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAAt                                                                <  1:207230/100‑1 (MQ=255)
                                 ttATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAAtt                                                               <  2:69287/100‑1 (MQ=255)
                                 ttATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAAtt                                                               >  2:2497824/1‑100 (MQ=255)
                                  tATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAAttt                                                              <  1:2574511/100‑1 (MQ=255)
                                    ttCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTc                                                            <  2:460903/100‑1 (MQ=255)
                                      ccAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTg                                                          >  1:2419335/1‑100 (MQ=255)
                                       cAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGt                                                         <  2:15402/100‑1 (MQ=255)
                                         gaagaGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTg                                                       <  1:841930/100‑1 (MQ=255)
                                           agagTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGAc                                                     <  1:208258/100‑1 (MQ=255)
                                           agagTTCAATATAAAGTTTGTCTTCAATTTTTCTTAATTATGAGTAATATAAATTAATCATTAAATAGTAAAATTCATCTTTTTTTAAttttattttt                                                       >  2:1206496/1‑92 (MQ=255)
                                            gagTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACt                                                    <  1:2153671/100‑1 (MQ=255)
                                              gTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTcc                                                  >  2:974910/1‑100 (MQ=255)
                                              gTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTcc                                                  >  1:2559505/1‑100 (MQ=255)
                                              gTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCCTTATTTAATTTTCTGTTGACTcc                                                  <  2:1785465/100‑1 (MQ=255)
                                                 cAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGc                                               <  2:1507009/100‑1 (MQ=255)
                                                 cAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGc                                               <  1:1839976/100‑1 (MQ=255)
                                                 cAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGc                                               >  2:1362433/1‑100 (MQ=255)
                                                  aaTATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGct                                              >  2:2817839/1‑100 (MQ=255)
                                                      taATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTcc                                          >  1:519072/1‑100 (MQ=255)
                                                      taATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTcc                                          >  1:1105644/1‑100 (MQ=255)
                                                        aTGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGt                                        <  2:278118/100‑1 (MQ=255)
                                                         tGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTg                                       <  2:340413/100‑1 (MQ=255)
                                                          gTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGa                                      >  2:1095930/1‑100 (MQ=255)
                                                              gTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAAc                                  >  2:2837112/1‑100 (MQ=255)
                                                              gTCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAAc                                  <  2:2615518/100‑1 (MQ=255)
                                                               tCTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACg                                 >  1:83694/1‑100 (MQ=255)
                                                                cTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACgg                                <  1:2656533/100‑1 (MQ=255)
                                                                cTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACgg                                >  1:2671805/1‑100 (MQ=255)
                                                                cTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACgg                                >  1:1774785/1‑100 (MQ=255)
                                                                cTTCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGATTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACg                                 >  2:805649/1‑99 (MQ=255)
                                                                 ttCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGt                               >  2:984802/1‑100 (MQ=255)
                                                                 ttCAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGt                               >  1:2045835/1‑100 (MQ=255)
                                                                   cAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGttt                             <  1:1362438/100‑1 (MQ=255)
                                                                   cAATTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGttt                             <  1:2383195/100‑1 (MQ=255)
                                                                    aaTTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGtttt                            <  1:917699/100‑1 (MQ=255)
                                                                     aTTTTTCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTta                           >  1:1270495/1‑100 (MQ=255)
                                                                        tttCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATaa                        <  1:322865/100‑1 (MQ=255)
                                                                        tttCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATaa                        >  1:2707313/1‑100 (MQ=255)
                                                                          tCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAAtt                      <  2:1624935/100‑1 (MQ=255)
                                                                          tCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAAtt                      <  2:2541423/100‑1 (MQ=255)
                                                                          tCTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTGTATAAtt                      <  1:1521857/100‑1 (MQ=255)
                                                                           cTTACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTa                     >  1:2095168/1‑100 (MQ=255)
                                                                             tACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGa                   <  2:2257794/100‑1 (MQ=255)
                                                                             tACTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGa                   >  1:2674071/1‑100 (MQ=255)
                                                                             tACCTCAGGGTAATCTAGATCGCTCATTACCTTGTGCGCTTCATCTTTATTTAATTTTCCGTTGACTCCAGCTCTCTGTGATAACGGTTTTATAAt                       >  2:172989/1‑96 (MQ=255)
                                                                               cTTCAGGGTAATATAGATCGCTTATTACATTATGAGCCTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATg                 <  2:1946821/100‑1 (MQ=255)
                                                                               cTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATg                 <  1:154272/100‑1 (MQ=255)
                                                                               cTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATg                 >  1:2728202/1‑100 (MQ=255)
                                                                               cTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATg                 <  1:1177437/100‑1 (MQ=255)
                                                                               cTTCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATg                 >  1:1061367/1‑100 (MQ=255)
                                                                                ttCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGc                <  2:890887/100‑1 (MQ=255)
                                                                                 tCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCt               >  1:1145984/1‑100 (MQ=255)
                                                                                 tCAGGGTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCt               <  1:1809412/100‑1 (MQ=255)
                                                                                 tCAGGGTAATATAGATCACTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCt               <  1:513897/100‑1 (MQ=255)
                                                                                    gggTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTAt            >  1:1981554/1‑100 (MQ=255)
                                                                                    cggTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTACTTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTAt            >  1:12546/2‑100 (MQ=255)
                                                                                      gTAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATcc          >  2:1217899/1‑100 (MQ=255)
                                                                                       tAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATccc         >  1:1600692/1‑100 (MQ=255)
                                                                                       tAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATccc         >  1:1032852/1‑100 (MQ=255)
                                                                                       tAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATccc         >  1:1600702/1‑100 (MQ=255)
                                                                                       tAATATAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATccc         >  1:2050581/1‑100 (MQ=255)
                                                                                         atatAGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCaa       <  2:2625878/100‑1 (MQ=255)
                                                                                          tataGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCaaa      >  2:1375434/1‑100 (MQ=255)
                                                                                          tataGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCaaa      >  2:251212/1‑100 (MQ=255)
                                                                                           ataGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAg     <  2:530556/100‑1 (MQ=255)
                                                                                           ataGATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAg     <  2:1972567/100‑1 (MQ=255)
                                                                                              gATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAGata  >  2:750961/1‑100 (MQ=255)
                                                                                              gATCGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAGata  >  2:439406/1‑100 (MQ=255)
                                                                                                 |                                                                                                
GAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTACTTCAGGGTAATATAGATTGCTCATTACATTGTGAGCTTCATCTTTATTTAATTTTCTGTTGACTCCAGCTCTCCGTGATAACGGTTTTATAATTAGATGCTTATCCCAAAGATA  >  NC_000913/1197756‑1197949

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: