Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,590,854 C→A 5.9% E5* (GAA→TAA)  ypfM ← stress‑induced small enterobacterial protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,590,8540CA 5.9% 226.1 / 2.2 68E5* (GAA→TAA) ypfMstress‑induced small enterobacterial protein
Reads supporting (aligned to +/- strand):  ref base C (21/43);  new base A (2/2);  total (23/45)
Fisher's exact test for biased strand distribution p-value = 5.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

AAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGACAAATCTTATGTAAATACTATGGTCCTGCGG  >  NC_000913/2590718‑2590991
                                                                                                                                        |                                                                                                                                         
aaGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACg                                                                                                                                        <  1:1146643/140‑1 (MQ=255)
    cGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCa                                                                                                                                    <  1:5277/140‑1 (MQ=255)
    cGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCa                                                                                                                                    <  1:5309/140‑1 (MQ=255)
            cAATGTGACGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATaac                                                                                                                            >  1:1543763/1‑140 (MQ=255)
                  gtCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACctc                                                                                                                      >  1:460095/1‑140 (MQ=255)
                     gCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACctctct                                                                                                                   >  1:859002/1‑140 (MQ=255)
                     gCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACctctct                                                                                                                   >  1:141009/1‑140 (MQ=255)
                     gCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAAAAACctctct                                                                                                                   >  1:858962/1‑140 (MQ=255)
                      ccTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCtt                                                                                                                  <  2:13462/140‑1 (MQ=255)
                      ccTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCtt                                                                                                                  <  2:96465/140‑1 (MQ=255)
                       cTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTta                                                                                                                 <  1:96621/140‑1 (MQ=255)
                       cTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTta                                                                                                                 <  1:453321/140‑1 (MQ=255)
                       cTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTta                                                                                                                 <  1:1662287/140‑1 (MQ=255)
                                cTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATgg                                                                                                        <  1:1196758/140‑1 (MQ=255)
                                   tCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCGTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGtg                                                                                                     >  1:1097617/1‑139 (MQ=255)
                                           aGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACgg                                                                                             <  1:785883/140‑1 (MQ=255)
                                           aGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACgg                                                                                             <  2:1265980/140‑1 (MQ=255)
                                           aGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACgg                                                                                             <  2:769390/140‑1 (MQ=255)
                                           aGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCAAGTTCCCCAGTTCACGTTCAATCATAACAACCTCCCTTATAATTATGGGTATTCTACgg                                                                                             <  2:353414/140‑1 (MQ=255)
                                            gAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTAGATAAAGTCTTTCCAGTTCCCCATTTCACGTTCAATCATAAAAACCTCTCTTATAATTATGGGTATTCTACGGa                                                                                            <  2:677674/140‑1 (MQ=255)
                                              aCCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGaaa                                                                                          <  2:17667/140‑1 (MQ=255)
                                                 tGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACaa                                                                                       >  2:1602892/1‑140 (MQ=255)
                                                     cAGCGGGTTCTTGCTGTTTGTACTTTGTCTAAGGAATTACTTACGAATCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAAtata                                                                                   <  1:1457263/140‑1 (MQ=255)
                                                       gCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATAcc                                                                                 <  1:1419312/140‑1 (MQ=255)
                                                            ttCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCGAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCgtggt                                                                            <  2:1389680/140‑1 (MQ=255)
                                                             tCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCgtggtg                                                                           >  1:201253/1‑140 (MQ=255)
                                                             tCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATAATACGGAAACAATATACCttggtg                                                                           >  1:201141/1‑140 (MQ=255)
                                                                      ttGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAAtt                                                                  <  2:1162063/140‑1 (MQ=255)
                                                                           cTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCGTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTc                                                             >  1:1639746/1‑140 (MQ=255)
                                                                             ttGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCTGTACACGTTCAATCATAACAACCTCTCTTAGAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGa                                                           >  2:1186126/1‑140 (MQ=255)
                                                                             ttGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGa                                                           >  1:1049999/1‑140 (MQ=255)
                                                                              tGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCACTCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGAt                                                          <  1:638900/140‑1 (MQ=255)
                                                                                     ggAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCg                                                   <  2:221989/140‑1 (MQ=255)
                                                                                          tacttacGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGAc                                              <  2:1386630/140‑1 (MQ=255)
                                                                                           acttacGAAGGATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCACTCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATAGACCGTGGTGAAGCTAATTTACTCGATTGCTGCGAGGACt                                             <  2:384331/140‑1 (MQ=255)
                                                                                             ttacGACGCATAACTTCGATAAAGTCTTTCCAGTACCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTAc                                           <  2:16518/140‑1 (MQ=255)
                                                                                              tacGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTAACGTTCAATCATAACAACCTCTCTTATAAt                                                                                                             >  2:1484754/1‑73 (MQ=255)
                                                                                              tacGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTAACGTTCAATCATAACAACCTCTCTTATAAt                                                                                                             <  1:1074585/73‑1 (MQ=255)
                                                                                              tacGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTAACGTTCAATCATAACAACCTCTCTTATAAt                                                                                                             >  2:1074585/1‑73 (MQ=255)
                                                                                              tacGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTAACGTTCAATCATAACAACCTCTCTTATAAt                                                                                                             <  1:1484754/73‑1 (MQ=255)
                                                                                                     cATAAATTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCgggg                                   <  2:954386/140‑1 (MQ=255)
                                                                                                       tAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGa                                 <  2:1533076/140‑1 (MQ=255)
                                                                                                       tAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGa                                 <  2:1226898/140‑1 (MQ=255)
                                                                                                        aaCTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTTAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGAc                                <  1:1360234/140‑1 (MQ=255)
                                                                                                              gATAAAGTCTTTCCAGTTCCCCAGTGCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGACAAATCt                          <  1:1206591/140‑1 (MQ=255)
                                                                                                              gATAAAGTCACTCCCGTTCCCCCGTTCACGTTCAATCATAACAACCTCTCTTATCATTATGGGTATTCTACGGAAACAATATACCGTGGTGATGCTAATTTACACGATTGCTGCGATGACGACCTCCGGGTGACAAATCt                          <  2:855196/140‑1 (MQ=255)
                                                                                                               aTAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGACAAATCtt                         >  2:732782/1‑140 (MQ=255)
                                                                                                               aTAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGACAAATCtt                         >  2:652114/1‑140 (MQ=255)
                                                                                                               aTAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGACAAATCtt                         >  1:1454802/1‑140 (MQ=255)
                                                                                                                ttaaGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGACAAATCTTa                        <  2:933442/138‑1 (MQ=255)
                                                                                                                tAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTATTTTACTCGATTGCTGCGATGACTACCTCCGGGGGACAAATCTTa                        <  2:935560/140‑1 (MQ=255)
                                                                                                                   aGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGACAAATCTTATGt                     <  2:648580/140‑1 (MQ=255)
                                                                                                                   aGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGACAAATCTTATGt                     <  2:841742/140‑1 (MQ=255)
                                                                                                                   aGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGACAAATCTTATGt                     <  1:187217/140‑1 (MQ=255)
                                                                                                                   aGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGACAAATCTTATGt                     <  2:1463376/140‑1 (MQ=255)
                                                                                                                   aGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGACAAATCTTATGt                     <  1:761554/140‑1 (MQ=255)
                                                                                                                    gTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGACAAATCTTATGTa                    >  2:496558/1‑140 (MQ=255)
                                                                                                                     tCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCTATGACTACCTCCGGGGGACAAATCTTATGTaa                   >  2:962905/1‑140 (MQ=255)
                                                                                                                     tCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGACAAATCTTATGTaa                   >  2:672492/1‑140 (MQ=255)
                                                                                                                     tCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGACAAATCTTATGTaa                   >  2:731171/1‑140 (MQ=255)
                                                                                                                     tCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTATACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGACAAATCTTATGTaa                   <  2:743159/140‑1 (MQ=255)
                                                                                                                      cTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGACAAATCTTATGTaaa                  <  1:1501953/140‑1 (MQ=255)
                                                                                                                      cTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGACAAATCTTATGTaaa                  <  1:1501981/140‑1 (MQ=255)
                                                                                                                      cTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGACAAATCTTATGTaaa                  <  1:1625594/140‑1 (MQ=255)
                                                                                                                               tCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGACAAATCTTATGTAAATACTATGGt         >  2:1135041/1‑140 (MQ=255)
                                                                                                                                   cAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGACAAATCTTATGTAAATACTATGGTCCTg     <  2:1521995/140‑1 (MQ=255)
                                                                                                                                     gTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGACAAATCTTATGTAAATACTATGGTCCTGCg   <  1:165794/140‑1 (MQ=255)
                                                                                                                                      ttCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGACAAATCTTATGTAAATACTATGGTCCTGCgg  >  1:1343355/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                         
AAGGCGACATGCCAATGTGTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTATAATTATGGGTATTCTACGGAAACAATATACCGTGGTGAAGCTAATTTACTCGATTGCTGCGATGACTACCTCCGGGGGACAAATCTTATGTAAATACTATGGTCCTGCGG  >  NC_000913/2590718‑2590991

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