Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,530,076:1 (A)7→8 100% coding (155/483 nt) ydcD → putative immunity protein for RhsE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,530,0691.A95.8% 474.5 / 12.1 142coding (148/483 nt)ydcDputative immunity protein for RhsE
Reads supporting (aligned to +/- strand):  ref base . (2/4);  new base A (78/58);  total (80/62)
Fisher's exact test for biased strand distribution p-value = 4.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 159 total aligned reads displayed.

TGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGC  >  NC_000913/1529970‑1530162
                                                                                                     |                                                                                              
    tATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGaaaa                                                                                           >  1:3663804/1‑100 (MQ=255)
     aTTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGAACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGaaaaa                                                                                          >  2:2656829/1‑100 (MQ=255)
       ttgttgTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAACCACCAAAGAGCGTTGAGATATATCTTTCTGaaaaaa                                                                                         >  2:141364/1‑99 (MQ=255)
        tgttgtTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGaaaaaaaa                                                                                       >  2:1992623/1‑99 (MQ=255)
        tgttgtTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGaaaaaaaa                                                                                       >  1:3313593/1‑99 (MQ=255)
          ttgttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGaaaaaaaata                                                                                     <  2:1622887/100‑4 (MQ=255)
           tgttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGaaaaaaaatat                                                                                    <  1:2669/100‑5 (MQ=255)
           tgttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGaaaaaaaatat                                                                                    >  2:439905/1‑96 (MQ=255)
            gttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGaaaaaaaatatt                                                                                   >  1:867548/1‑95 (MQ=255)
                aCCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATAttattt                                                                                >  1:3763221/1‑100 (MQ=255)
                 ccTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATAttattta                                                                               >  1:2389676/1‑100 (MQ=255)
                    tGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATaa                                                                            >  1:3244691/1‑100 (MQ=255)
                      ggTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATaaaa                                                                          <  1:2436654/100‑1 (MQ=255)
                      ggTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATaaaa                                                                          >  1:1809496/1‑100 (MQ=255)
                        tttGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATaaa                                                                           <  2:1986084/97‑1 (MQ=255)
                         ttGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAaataa                                                                       <  2:2156659/100‑1 (MQ=255)
                          tGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAaataat                                                                      >  2:904912/1‑100 (MQ=255)
                            aTAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAAtaatta                                                                    <  1:2943032/100‑1 (MQ=255)
                             ttccgatctaaaTTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAAtaattaa                                                                   <  1:1481607/91‑1 (MQ=255)
                             tAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAAtaattaa                                                                   >  1:3528272/1‑100 (MQ=255)
                              aaCAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAAtaattaat                                                                  <  1:315805/100‑1 (MQ=255)
                              aaCAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAAtaattaat                                                                  >  2:2893865/1‑100 (MQ=255)
                                 aaTTCAAATTCTGACATAGTATTTTATGGAAATAAAACACAAAAGAGCGTTGAGATATAAATATCTAAAAAAAAATATTATTTATAAAATAATTAATGa                                                                <  2:1445417/99‑1 (MQ=255)
                                  aTTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGAcc                                                              >  1:2404939/1‑100 (MQ=255)
                                   ttCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCa                                                             >  2:2409358/1‑100 (MQ=255)
                                   ttCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTGGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCa                                                             >  1:2914949/1‑100 (MQ=255)
                                     cAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACTaa                                                            >  2:4084079/1‑99 (MQ=255)
                                     cAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCaaa                                                           >  2:1078134/1‑100 (MQ=255)
                                     cAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCaaa                                                           >  2:270555/1‑100 (MQ=255)
                                     cAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCaaa                                                           >  2:1773451/1‑100 (MQ=255)
                                      aaaTTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAAtaattaaagatcggaa                                                          >  2:1481598/1‑91 (MQ=255)
                                       aaTTCTGTCATTTTATTTTATTTACAAAAAACAACAAAAAGTGTTAATATATAACTGCCTCAAAAAAAATATTATTTATAAAATAATTAATGAc                                                               <  2:834415/94‑1 (MQ=255)
                                          tCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAAtt                                                      <  1:159856/100‑1 (MQ=255)
                                          tCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAAtt                                                      <  1:421174/100‑1 (MQ=255)
                                           cTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATta                                                     >  2:1865466/1‑100 (MQ=255)
                                           cTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATta                                                     <  1:3198952/100‑1 (MQ=255)
                                            tGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtag                                                    <  1:3496487/100‑1 (MQ=255)
                                            tGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtag                                                    >  2:3055271/1‑100 (MQ=255)
                                            ctctttggtgttttATGTCCATAAAATACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtag                                                    >  2:662259/11‑100 (MQ=25)
                                              aCATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCaaaaaaagatcg                                                  >  1:440224/1‑94 (MQ=255)
                                               catcatATTTTATGGACATAAAACACCAAAGAGCGTTTAGATATATCTCTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtagtag                                                 >  2:847001/1‑100 (MQ=255)
                                                aTCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTaga                                                >  2:1639494/1‑100 (MQ=255)
                                                aTCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTaga                                                >  1:196817/1‑100 (MQ=255)
                                                aTCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACAAAAAAATTAATaga                                                >  2:505718/1‑100 (MQ=255)
                                                 tCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGGCCAAAAAATTAGTagag                                               >  1:766145/1‑100 (MQ=255)
                                                 tCGTATTTTATGGACAAAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTagag                                               <  1:904917/100‑1 (MQ=255)
                                                  ggtgttttATGTCCATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAgg                                              >  1:1007706/5‑100 (MQ=255)
                                                  cGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAgg                                              >  1:3686176/1‑100 (MQ=255)
                                                   gTATTTTGTGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAggg                                             >  2:2180257/1‑100 (MQ=255)
                                                   gTATTTAATGGAAATAAAACACCAAAGAGCGTTGAGATATATATTAAAGAATAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTagag                                               <  2:2303901/98‑1 (MQ=255)
                                                     aTTTTATTGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTagag                                               <  2:2984590/96‑1 (MQ=255)
                                                     aTTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGaa                                           >  1:719186/1‑100 (MQ=255)
                                                     aTTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGaa                                           >  2:2590433/1‑100 (MQ=255)
                                                      ttttATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAAt                                          >  1:451616/1‑100 (MQ=255)
                                                      ttttATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAAt                                          >  1:2027839/1‑100 (MQ=255)
                                                      ttttATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGAAGAGGGaa                                           >  2:2126995/1‑99 (MQ=255)
                                                       tatATGGACATAAAAAACCAAAGAGGGTTGAAATATAAAAAAATAAAAAAAAATATTATTTATAAAATAATTAATGACAAAAAAATTAGTAGAGGGAAt                                          <  2:2837389/97‑1 (MQ=255)
                                                        ttATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATgg                                        <  1:1373122/100‑1 (MQ=255)
                                                          atgaaaaTAAAAAACAAAAGAGTGTTTAAATATAAATTTAAAAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATgg                                        <  2:959075/92‑1 (MQ=255)
                                                          aTGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTc                                      <  2:2390042/100‑1 (MQ=255)
                                                            ggACATAAAACACCAAAGAGCGTTGAGATATAAATTTCAGAAAAAAAATATTATTTATAAAATAAAAAAAAAAAAAAAAATTAGTAGAGGGAATGGTCAt                                    <  2:2495278/100‑1 (MQ=255)
                                                            ggACATAAAACACAAAAGAGCGTTGAGATATAACTTTATGAAAAAAAATATTATTTATAAAATAAAAAAGAAAAAAAAAATTAGTAGAgg                                              <  2:2544407/90‑1 (MQ=255)
                                                             gACATAAAACACTAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtt                                   >  1:2733425/1‑100 (MQ=255)
                                                             gACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtt                                   >  1:1386906/1‑100 (MQ=255)
                                                             gACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtt                                   >  1:726630/1‑100 (MQ=255)
                                                             gACATAAAACAACAAAGAGTGTTAAGATATAAATTTATAAAAAAAAATATTATTTATAAAATAAAAAAAAAAAAAAAAATTAGTAGAGGGAATgg                                        <  2:847227/95‑1 (MQ=255)
                                                              aCATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAttt                                  <  1:355888/100‑1 (MQ=255)
                                                              aCATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtt                                   <  2:3129222/99‑1 (MQ=255)
                                                               cATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtttt                                 >  2:632783/1‑100 (MQ=255)
                                                                aTAAAACACCAAAGAGCGTTAAGATATAAAAAAAAAAAAAAAAATATTATTTATAAGATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtttt                                 <  2:46570/99‑1 (MQ=255)
                                                                 tAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtat                               >  1:1997244/1‑100 (MQ=255)
                                                                 tAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtat                               >  2:896912/1‑100 (MQ=255)
                                                                  aaaaCACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtata                              >  2:2395989/1‑100 (MQ=255)
                                                                  aaaaCACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtata                              >  2:3473241/1‑100 (MQ=255)
                                                                  aaaaCACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtata                              >  1:2249465/1‑100 (MQ=255)
                                                                  aaaaCACCAAAGAGCGTTAAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtata                              <  2:3027319/100‑1 (MQ=255)
                                                                    aaCACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAg                            <  1:2024873/100‑1 (MQ=255)
                                                                    aaCACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAg                            <  1:3217234/100‑1 (MQ=255)
                                                                    aaCACCAAAGAGCGTCGAGATATATATTTATGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAg                            <  2:3972368/100‑1 (MQ=255)
                                                                     acacCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGt                           >  2:3915425/1‑100 (MQ=255)
                                                                      cacCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGta                          <  2:3143996/100‑1 (MQ=255)
                                                                          aaaGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAAt                      >  2:1145950/1‑100 (MQ=255)
                                                                           atctgCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATg                     <  2:2785202/96‑1 (MQ=255)
                                                                             gagCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGt                   <  2:4079539/100‑1 (MQ=255)
                                                                               gCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATgagat                 >  1:2785218/1‑96 (MQ=255)
                                                                                 gTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTATTGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaat               >  2:350116/1‑100 (MQ=255)
                                                                                 gTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaat               <  1:3866061/100‑1 (MQ=255)
                                                                                  ttGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGTGAATGGTTATTTTATAAGTATAATGGTtaata              >  2:1905569/1‑100 (MQ=255)
                                                                                    gagaTATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGAACAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaata              <  2:1146365/98‑1 (MQ=255)
                                                                                     agaTATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAAtt           >  1:1916970/1‑100 (MQ=255)
                                                                                     agaTATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATTGTCATTTTATAAGTATATTGGTTaataatca          <  1:2808223/100‑3 (MQ=255)
                                                                                       atatatCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTAc         <  2:1052273/100‑1 (MQ=255)
                                                                                        tatatCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACa        >  1:2162148/1‑100 (MQ=255)
                                                                                         atatTTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAg       <  2:2242070/100‑1 (MQ=255)
                                                                                         atatCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACa        <  1:4038512/99‑1 (MQ=255)
                                                                                         atatCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAg       <  1:955248/100‑1 (MQ=255)
                                                                                         atatCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAg       <  1:1237477/100‑1 (MQ=255)
                                                                                          tatCTTTCTGAAAAAAAATATTATTTATAAAATATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGa     >  2:2796653/1‑100 (MQ=255)
                                                                                              tttCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGCAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGc  >  1:56337/1‑100 (MQ=255)
                                                                                               ttCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTATTAGTAGGAATTGTTATTTTAAAAGTATAATAGTTAATAATTACAGGACg   >  2:3653887/1‑98 (MQ=255)
                                                                                                     |                                                                                              
TGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGC  >  NC_000913/1529970‑1530162

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: