Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,197,739 A→G 100% S167P (TCC→CCC)  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,7390AG97.8% 343.3 / ‑3.1 105S167P (TCC→CCC) ymfEe14 prophage; putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base G (56/44);  total (58/47)
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.99e-01

Only 100 of 117 total aligned reads displayed.

ACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTAC  >  NC_000913/1197642‑1197835
                                                                                                 |                                                                                                
 ctGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAAt                                                                                               >  2:222365/1‑100 (MQ=255)
 ctGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAAt                                                                                               >  2:2033606/1‑100 (MQ=255)
  tGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAAtt                                                                                              >  2:2276852/1‑100 (MQ=255)
   gagaATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTa                                                                                             <  2:1336052/100‑1 (MQ=255)
   gagaATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTa                                                                                             >  1:558198/1‑100 (MQ=255)
     gaATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGAGAATTAAt                                                                                           >  1:2050532/1‑100 (MQ=255)
      aaTGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATg                                                                                          >  2:3310075/1‑100 (MQ=255)
       aTGATAAAATTACAAAAAAGAATAAAGGCAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGc                                                                                         <  2:2865582/100‑1 (MQ=255)
       aTGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGc                                                                                         <  1:2482380/100‑1 (MQ=255)
        tGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCa                                                                                        <  2:427414/100‑1 (MQ=255)
        tGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCa                                                                                        <  1:1693128/100‑1 (MQ=255)
        tGAAAAAATTACAAAAAAGAATAAAGAAAGCAAAGAAAATTAAAAATAAAAAATATCCATTTGTTTGAAAAAAAAGATTATTAGGAAGGGAATTAATGCa                                                                                        <  1:1674467/100‑1 (MQ=255)
         gATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCaa                                                                                       >  2:2511032/1‑100 (MQ=255)
         gATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCAATTTGTTTGCAAAAAAAGATTATTAGGAAGGAAATTAATGCaa                                                                                       >  1:68652/1‑100 (MQ=255)
          aTAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAAt                                                                                      >  2:387462/1‑100 (MQ=255)
          aTAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAAt                                                                                      <  2:267828/100‑1 (MQ=255)
           tAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAAtt                                                                                     >  2:443822/1‑100 (MQ=255)
             aaaTTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTAt                                                                                   <  2:1143136/100‑1 (MQ=255)
             aaaTTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTAt                                                                                   >  2:147364/1‑100 (MQ=255)
              aaTTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATc                                                                                  <  2:632294/100‑1 (MQ=255)
               aTTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGCAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCg                                                                                 >  1:1854392/1‑100 (MQ=255)
               aTTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCg                                                                                 <  2:2210493/100‑1 (MQ=255)
               aTTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCg                                                                                 <  2:3121823/100‑1 (MQ=255)
                ttACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGa                                                                                <  2:2918943/100‑1 (MQ=255)
                 tACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGaa                                                                               >  2:45361/1‑100 (MQ=255)
                   cAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGaaaa                                                                             <  2:1464444/100‑1 (MQ=255)
                   cAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGaaaa                                                                             >  2:3094512/1‑100 (MQ=255)
                    aaaaaaGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGaaaa                                                                             >  1:1676876/1‑99 (MQ=255)
                    aaaaaaGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGaaaa                                                                             >  1:1423303/1‑99 (MQ=255)
                    aaaaaaGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAAt                                                                            >  1:2404847/1‑100 (MQ=255)
                            aTAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCaaaaaa                                                                    >  1:1158881/1‑100 (MQ=255)
                            aTAAAGACAAAAAAGAAACAAAAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTATGAAGGGAATTAATTCAATTATAaaaaaaaaaaaaaa                                                                    <  1:2834600/100‑1 (MQ=255)
                             tAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAat                                                                   <  2:2211306/100‑1 (MQ=255)
                               aaGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAata                                                                  >  1:1481598/1‑99 (MQ=255)
                                aGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAGAATTCAAAAAATATc                                                                >  2:3339298/1‑100 (MQ=255)
                                aGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATc                                                                <  2:3008891/100‑1 (MQ=255)
                                aGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGAATAATAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAata                                                                  >  1:3167007/1‑98 (MQ=255)
                                  aCAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCa                                                              <  2:2150974/100‑1 (MQ=255)
                                   cAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCaa                                                             >  1:1259003/1‑100 (MQ=255)
                                     gCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCaaaa                                                           >  1:1161083/1‑100 (MQ=255)
                                     gCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCaaaa                                                           >  2:2085153/1‑100 (MQ=255)
                                     gCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCaaaa                                                           >  1:82434/1‑100 (MQ=255)
                                     gCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCaaaa                                                           >  2:2564477/1‑100 (MQ=255)
                                       aaaGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAAt                                                         >  2:728008/1‑100 (MQ=255)
                                         aGAACCTTCAACATGAAAAATATCCAAATGAAAAAAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAg                                                       <  1:2559478/100‑1 (MQ=255)
                                             ccTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGtata                                                   >  2:901646/1‑100 (MQ=255)
                                             ccTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGtata                                                   >  2:3275335/1‑100 (MQ=255)
                                              ctacaacaaAAAAAATATCCATTTAAATGAAAAAAAAGATTATTAGGAAGTGAATTAATGCAATTATAGAAAATTCAAAAAATATCCAAAAATAGTATAc                                                  <  1:2439813/97‑1 (MQ=255)
                                              cTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATAc                                                  >  1:1008946/1‑100 (MQ=255)
                                              cTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATAc                                                  <  2:2093338/100‑1 (MQ=255)
                                               ttaaaCATGAAAAATATCCATTTAATAGAAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATAc                                                  <  1:918467/96‑1 (MQ=255)
                                                tCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACtt                                                >  1:1731172/1‑100 (MQ=255)
                                                    cATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGAAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTAt                                             >  2:804323/1‑99 (MQ=255)
                                                      tGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTcc                                          >  2:2423254/1‑100 (MQ=255)
                                                      tGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTcc                                          <  1:1181675/100‑1 (MQ=255)
                                                      tGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTcc                                          <  2:403120/100‑1 (MQ=255)
                                                      tGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTc                                           >  2:3032044/1‑99 (MQ=255)
                                                      tGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACGTTATTcc                                          <  2:2520625/100‑1 (MQ=255)
                                                       gAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCa                                         <  2:2931465/100‑1 (MQ=255)
                                                       gAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGCAAATTAATCCAATTATCAAAAATTCAAAAAATATCCAAAAATAGTATACTATATTCCa                                         >  1:261201/1‑100 (MQ=255)
                                                        aaaaaTATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAAAATAAAAAAATAGTATACTTTATTCCa                                         <  1:935535/99‑1 (MQ=255)
                                                         aaaaTATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCaga                                       >  1:1793828/1‑100 (MQ=255)
                                                         aaaaTATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCaga                                       >  1:994944/1‑100 (MQ=255)
                                                             tatCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCagaagan                                   <  2:2702422/100‑2 (MQ=255)
                                                             tatCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAagag                                   <  2:250609/100‑1 (MQ=255)
                                                              atCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGt                                  >  2:2147582/1‑100 (MQ=255)
                                                               tCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGtt                                 >  1:1325860/1‑100 (MQ=255)
                                                                ccATTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTc                                <  2:3035391/100‑1 (MQ=255)
                                                                  aTTTGTTTGCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCaa                              >  2:3386595/1‑100 (MQ=255)
                                                                    ttgtttgCAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAata                            <  2:2513061/100‑1 (MQ=255)
                                                                     tgtttgCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAatat                           <  2:1424010/100‑1 (MQ=255)
                                                                       tttgCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATaa                         <  2:106079/100‑1 (MQ=255)
                                                                       tttgCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATaa                         >  1:2401300/1‑100 (MQ=255)
                                                                        ttgCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAAt                        >  2:3390276/1‑100 (MQ=255)
                                                                         tgCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAAt                        >  1:1775359/1‑99 (MQ=255)
                                                                         tgCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATg                       >  2:648679/1‑100 (MQ=255)
                                                                         taaaaaaaaaGATTATTAGGAAGGTAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAat                             <  1:2804761/91‑1 (MQ=255)
                                                                          gCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGt                      <  1:133438/100‑1 (MQ=255)
                                                                          gCAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGt                      <  1:2722212/100‑1 (MQ=255)
                                                                           cAAAAAAAGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGtt                     >  1:2114345/1‑100 (MQ=255)
                                                                            aaaaaaaGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAAt                        <  1:988378/96‑1 (MQ=255)
                                                                             aaaaaaGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTg                   <  2:2748556/100‑1 (MQ=255)
                                                                             aaaaaaGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTg                   <  1:730580/100‑1 (MQ=255)
                                                                               aaaaGATTATTAGGAAGGGAATTAATCCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTc                 >  2:2060753/1‑100 (MQ=255)
                                                                                 aaGATTATTAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCtt               <  2:536056/100‑1 (MQ=255)
                                                                                    attattAGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCaa            >  2:3231748/1‑100 (MQ=255)
                                                                                     ttattaGGAAGGTAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAAt           <  1:1598517/100‑1 (MQ=255)
                                                                                     ttattaGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAAt           >  1:3067813/1‑100 (MQ=255)
                                                                                       attaGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAAttt         >  1:336196/1‑100 (MQ=255)
                                                                                         taGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAAttttt       >  1:2415452/1‑100 (MQ=255)
                                                                                         taGGAAGGGAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAAttttt       <  2:2132010/100‑1 (MQ=255)
                                                                                         taGGAAGGGAATTAATGCAACTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAAttttt       >  1:2074419/1‑100 (MQ=255)
                                                                                          agtaaggaaATTAATGCAATTATCTAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTc      <  1:1396019/97‑1 (MQ=255)
                                                                                          aggaagggAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTc      <  1:1919475/100‑1 (MQ=255)
                                                                                          aggaaaggaaTTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTc      <  1:2218234/100‑1 (MQ=255)
                                                                                           ggaagggaATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCt     >  2:2001565/1‑100 (MQ=255)
                                                                                           ggaagggaATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCt     >  1:1260248/1‑100 (MQ=255)
                                                                                            gaagggaATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCt     <  2:130945/99‑1 (MQ=255)
                                                                                             aagggaATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTAATTCCAGAAGAATTCAATATAATGTTTGTATTCAAttttt       >  1:111241/1‑96 (MQ=255)
                                                                                              agggaATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTAc  >  2:1869410/1‑100 (MQ=255)
                                                                                                 |                                                                                                
ACTGAGAATGATAAAATTACAAAAAAGAATAAAGACAGCAAAGAACCTTCAACATGAAAAATATCCATTTGTTTGCAAAAAAAGATTATTAGGAAGGAAATTAATGCAATTATCGAAAATTCAAAAAATATCCAAAAATAGTATACTTTATTCCAGAAGAGTTCAATATAATGTTTGTCTTCAATTTTTCTTAC  >  NC_000913/1197642‑1197835

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: ‑