Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 657,246 A→G intergenic (‑85/‑8) crcB ← / → ybeH hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913657,2460AG98.6% 658.2 / ‑1.3 210intergenic (‑85/‑8)crcB/ybeHhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (82/124);  total (83/127)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

Only 100 of 228 total aligned reads displayed.

TAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGG  >  NC_000913/657156‑657335
                                                                                          |                                                                                         
tAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACaggggg                                                                                         <  1:4251955/93‑1 (MQ=255)
tAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACaggggg                                                                                         >  2:3200937/1‑93 (MQ=255)
 aaCACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACggggggg                                                                                        >  2:1036627/1‑92 (MQ=255)
     cTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCa                                                                                    >  1:3030665/1‑93 (MQ=255)
      tGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCaa                                                                                   >  1:2954893/1‑93 (MQ=255)
      tGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCaa                                                                                   <  1:1375720/93‑1 (MQ=255)
      tGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCaa                                                                                   <  2:3053287/93‑1 (MQ=255)
       gCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAAt                                                                                  >  2:2761091/1‑93 (MQ=255)
       gCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAAt                                                                                  <  1:1910401/93‑1 (MQ=255)
         aaaaTCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGc                                                                                <  2:745329/93‑1 (MQ=255)
         aaaaTCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGc                                                                                <  2:406260/93‑1 (MQ=255)
          aaaTCCTGCTATTTTATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCt                                                                               <  1:1589410/93‑1 (MQ=255)
          aaaTCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCt                                                                               <  2:2300223/93‑1 (MQ=255)
          aaaTCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCt                                                                               >  1:3547568/1‑93 (MQ=255)
          aaaTCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTAGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCt                                                                               <  2:589223/93‑1 (MQ=255)
           aaTCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTg                                                                              <  2:4129455/93‑1 (MQ=255)
             tCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGt                                                                            <  2:2172551/93‑1 (MQ=255)
             tCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGt                                                                            >  1:1050156/1‑93 (MQ=255)
             tCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGt                                                                            >  1:2019239/1‑93 (MQ=255)
               cTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTg                                                                          >  1:1659862/1‑93 (MQ=255)
                tGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGc                                                                         >  2:2610743/1‑93 (MQ=255)
                tGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGc                                                                         <  1:2561899/93‑1 (MQ=255)
                tGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGc                                                                         >  1:4715015/1‑93 (MQ=255)
                tGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGCTGc                                                                         >  2:4493899/1‑93 (MQ=255)
                   tATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGc                                                                      >  2:192066/1‑93 (MQ=255)
                   tATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGc                                                                      <  2:2880335/93‑1 (MQ=255)
                    atttgatttgTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCa                                                                     <  2:2083818/93‑1 (MQ=255)
                    atttgatttgTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCa                                                                     >  2:1668927/1‑93 (MQ=255)
                    atttgatttgTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCa                                                                     <  2:874509/93‑1 (MQ=255)
                    atttgatttgTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCa                                                                     <  1:3487567/93‑1 (MQ=255)
                        gatttgTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGAc                                                                 <  1:3968639/93‑1 (MQ=255)
                        gatttgTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGAc                                                                 <  2:3107336/93‑1 (MQ=255)
                        gatttgTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTACCGCAGGAc                                                                 >  2:3719399/1‑93 (MQ=255)
                         atttgTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACa                                                                >  1:2205848/1‑93 (MQ=255)
                         atttgTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACa                                                                >  1:3618131/1‑93 (MQ=255)
                         atttgTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACa                                                                >  1:2443389/1‑93 (MQ=255)
                          tttgTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACat                                                               >  2:1451416/1‑92 (MQ=255)
                            tgTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGtt                                                             <  1:2348710/93‑1 (MQ=255)
                             gTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGttt                                                            <  2:2620864/93‑1 (MQ=255)
                             gTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGttt                                                            <  1:4504587/93‑1 (MQ=255)
                              tATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTg                                                           <  2:4228539/93‑1 (MQ=255)
                                tGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCt                                                         <  1:3107949/93‑1 (MQ=255)
                                tGAGCGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCt                                                         <  2:2470820/93‑1 (MQ=255)
                                  aGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGt                                                       <  2:3987161/93‑1 (MQ=255)
                                   gTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGtt                                                      >  1:1970572/1‑93 (MQ=255)
                                     gATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGTCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTAc                                                    >  1:2904497/1‑93 (MQ=255)
                                     gATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTAc                                                    <  2:3506047/93‑1 (MQ=255)
                                     gATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTAc                                                    >  2:4165625/1‑93 (MQ=255)
                                     gATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTAc                                                    >  1:3917321/1‑93 (MQ=255)
                                        aaGTGTAACGCCGAATAATCGTCGTTGGCTAATTTGACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATc                                                 <  1:159349/93‑1 (MQ=255)
                                        aaGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATc                                                 <  1:971209/93‑1 (MQ=255)
                                        aaGTGTAACGCCCAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATc                                                 <  2:2417950/93‑1 (MQ=255)
                                         aGTGTAACGCCGAATAATCGTCGTTGGTGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCt                                                <  1:3747434/93‑1 (MQ=255)
                                           tgtAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCtgt                                              <  1:2241802/93‑1 (MQ=255)
                                           tgtAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCtgt                                              <  1:2747296/93‑1 (MQ=255)
                                           tgtAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCtgt                                              <  2:1964605/93‑1 (MQ=255)
                                            gtAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCtgtg                                             <  2:1568959/93‑1 (MQ=255)
                                               aCGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTgg                                          >  2:2132035/1‑93 (MQ=255)
                                                  ccGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGaa                                       <  2:2232571/93‑1 (MQ=255)
                                                   cGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGaaa                                      >  2:3482670/1‑93 (MQ=255)
                                                    gAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGaaaa                                     <  1:1782867/93‑1 (MQ=255)
                                                    gAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGaaaa                                     >  2:2189187/1‑93 (MQ=255)
                                                       taatCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGTGAAAAGaa                                  <  2:1300996/93‑1 (MQ=255)
                                                       taatCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGaa                                  <  1:2477848/93‑1 (MQ=255)
                                                       taatCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGaa                                  <  1:3325830/93‑1 (MQ=255)
                                                        aatCGTCGTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAAc                                 <  1:3776192/93‑1 (MQ=255)
                                                          tcgtcgTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGTCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGc                               <  2:2454254/93‑1 (MQ=255)
                                                          tcgtcgTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGc                               >  1:1042909/1‑93 (MQ=255)
                                                          tcgtcgTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGc                               >  1:1498650/1‑93 (MQ=255)
                                                          tcgtcgTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGc                               >  2:4153690/1‑93 (MQ=255)
                                                          tcgtcgTTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGc                               >  2:4425648/1‑93 (MQ=255)
                                                           cgtcgtTGGCGAATTTTACGACTCTGACAGGGGGGGGCAATGCTGGTTGCCGCAGGGCAGGTTGCTGTTACATCTGTGTGGGAAAAGAACGCt                              >  2:4664967/1‑93 (MQ=255)
                                                             tcgtTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTga                            <  2:4270477/93‑1 (MQ=255)
                                                             tcgtTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTga                            <  2:3908403/93‑1 (MQ=255)
                                                             tcgtTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGAAAAAGAACGCTga                            <  2:2861818/93‑1 (MQ=255)
                                                              cgtTGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTgag                           <  1:6600/93‑1 (MQ=255)
                                                                ttGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGAt                         <  1:787466/93‑1 (MQ=255)
                                                                ttGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGAt                         <  1:4058708/93‑1 (MQ=255)
                                                                 tGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGAtt                        >  2:1840727/1‑93 (MQ=255)
                                                                 tGGCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGAtt                        >  2:1597616/1‑93 (MQ=255)
                                                                   gCGAATTTTGCGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGAAAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTtg                      <  1:3333306/93‑1 (MQ=255)
                                                                   gCGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTtg                      <  2:1093069/93‑1 (MQ=255)
                                                                    cGAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTtgt                     >  1:551329/1‑93 (MQ=255)
                                                                     gAATTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTtgtg                    <  1:4629810/93‑1 (MQ=255)
                                                                       aTTTTACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGcc                  <  1:3401110/93‑1 (MQ=255)
                                                                        ttttACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCt                 <  2:3587379/93‑1 (MQ=255)
                                                                         tttACGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTc                <  2:1364613/93‑1 (MQ=255)
                                                                           tACGACGCGGGCAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGt              <  2:3283725/93‑1 (MQ=255)
                                                                            aCGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGtt             <  1:327895/93‑1 (MQ=255)
                                                                            aCGACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGtt             <  1:4185684/93‑1 (MQ=255)
                                                                              gACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGa           >  1:3872044/1‑93 (MQ=255)
                                                                              gACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGa           <  1:1305935/93‑1 (MQ=255)
                                                                              gACTCTGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGa           <  1:1175831/93‑1 (MQ=255)
                                                                                  ctGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGgc       <  2:4404776/93‑1 (MQ=255)
                                                                                  ctGACAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGgc       <  2:519793/93‑1 (MQ=255)
                                                                                     aCAGGTGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCa    <  2:3652273/93‑1 (MQ=255)
                                                                                      cAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAg   >  2:4064948/1‑93 (MQ=255)
                                                                                      cAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAg   <  1:1489979/93‑1 (MQ=255)
                                                                                      cAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAg   <  2:4600481/93‑1 (MQ=255)
                                                                                      cAGGGGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAg   <  1:1685433/93‑1 (MQ=255)
                                                                                          |                                                                                         
TAACACTGCAAAATCCTGCTATTTGATTTGTATGAGTGATAAGTGTAACGCCGAATAATCGTCGTTGGCGAATTTTACGACTCTGACAGGAGGTGGCAATGCTGGTTGCCGCAGGACAGTTTGCTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGG  >  NC_000913/657156‑657335

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