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.

TGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑G‑AA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGC  >  NC_000913/1529970‑1530162
                                                                                                     |                                                                                              
    tATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑G‑aa‑aa                                                                                           >  1:3663804/1‑100 (MQ=255)
     aTTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGAACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑G‑aa‑aaa                                                                                          >  2:2656829/1‑100 (MQ=255)
       ttgttgTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAACCACCAAAGAGCGTTGAGATATATCTTTCT‑G‑aa‑aaaa                                                                                         >  2:141364/1‑99 (MQ=255)
        tgttgtTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑G‑aa‑aaaaaa                                                                                       >  2:1992623/1‑99 (MQ=255)
        tgttgtTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑G‑aa‑aaaaaa                                                                                       >  1:3313593/1‑99 (MQ=255)
          ttgttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑G‑aa‑aaaaaata                                                                                     <  2:1622887/100‑4 (MQ=255)
           tgttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑G‑aa‑aaaaaatat                                                                                    <  1:2669/100‑5 (MQ=255)
           tgttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑G‑aa‑aaaaaatat                                                                                    >  2:439905/1‑96 (MQ=255)
            gttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑G‑aa‑aaaaaatatt                                                                                   >  1:867548/1‑95 (MQ=255)
                aCCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATAttattt                                                                                >  1:3763221/1‑100 (MQ=255)
                 ccTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATAttattta                                                                               >  1:2389676/1‑100 (MQ=255)
                    tGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATaa                                                                            >  1:3244691/1‑100 (MQ=255)
                      ggTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATaaaa                                                                          <  1:2436654/100‑1 (MQ=255)
                      ggTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATaaaa                                                                          >  1:1809496/1‑100 (MQ=255)
                        tttGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATaaa                                                                           <  2:1986084/97‑1 (MQ=255)
                         ttGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAaataa                                                                       <  2:2156659/100‑1 (MQ=255)
                          tGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAaataat                                                                      >  2:904912/1‑100 (MQ=255)
                            aTAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAAtaatta                                                                    <  1:2943032/100‑1 (MQ=255)
                             ttccgatctaaaTTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAAtaattaa                                                                   <  1:1481607/91‑1 (MQ=255)
                             tAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAAtaattaa                                                                   >  1:3528272/1‑100 (MQ=255)
                              aaCAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAAtaattaat                                                                  <  1:315805/100‑1 (MQ=255)
                              aaCAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAAtaattaat                                                                  >  2:2893865/1‑100 (MQ=255)
                                 aaTTCAAATTCTGACATAGTATTTTATGGAAATAAAACACAAAAGAGCGTTGAGATATAAATATCT‑AAAA‑AAAAATATTATTTATAAAATAATTAATGa                                                                <  2:1445417/99‑1 (MQ=255)
                                  aTTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGAcc                                                              >  1:2404939/1‑100 (MQ=255)
                                   ttCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCa                                                             >  2:2409358/1‑100 (MQ=255)
                                   ttCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTGGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCa                                                             >  1:2914949/1‑100 (MQ=255)
                                     cAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACTaa                                                            >  2:4084079/1‑99 (MQ=255)
                                     cAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCaaa                                                           >  2:1078134/1‑100 (MQ=255)
                                     cAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCaaa                                                           >  2:270555/1‑100 (MQ=255)
                                     cAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCaaa                                                           >  2:1773451/1‑100 (MQ=255)
                                      aaaTTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAAtaattaaagatcggaa                                                          >  2:1481598/1‑91 (MQ=255)
                                       aaTTCTGTCATTTTATTTTATTTACAAAAAACAACAAAAAGTGTTAATATATAACTGCCT‑CAAA‑AAAAATATTATTTATAAAATAATTAATGAc                                                               <  2:834415/94‑1 (MQ=255)
                                          tCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAAtt                                                      <  1:159856/100‑1 (MQ=255)
                                          tCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAAtt                                                      <  1:421174/100‑1 (MQ=255)
                                           cTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATta                                                     >  2:1865466/1‑100 (MQ=255)
                                           cTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATta                                                     <  1:3198952/100‑1 (MQ=255)
                                            tGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATtag                                                    <  1:3496487/100‑1 (MQ=255)
                                            tGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATtag                                                    >  2:3055271/1‑100 (MQ=255)
                                            ctctttggtgttttATGTCCATAAAATACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATtag                                                    >  2:662259/11‑100 (MQ=25)
                                              aCATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCaaaaaaagatcg                                                  >  1:440224/1‑94 (MQ=255)
                                               catcatATTTTATGGACATAAAACACCAAAGAGCGTTTAGATATATCTCTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATtagtag                                                 >  2:847001/1‑100 (MQ=255)
                                                aTCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTaga                                                >  2:1639494/1‑100 (MQ=255)
                                                aTCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTaga                                                >  1:196817/1‑100 (MQ=255)
                                                aTCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACAAAAAAATTAATaga                                                >  2:505718/1‑100 (MQ=255)
                                                 tCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGGCCAAAAAATTAGTagag                                               >  1:766145/1‑100 (MQ=255)
                                                 tCGTATTTTATGGACAAAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTagag                                               <  1:904917/100‑1 (MQ=255)
                                                  ggtgttttATGTCCATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAgg                                              >  1:1007706/5‑100 (MQ=255)
                                                  cGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAgg                                              >  1:3686176/1‑100 (MQ=255)
                                                   gTATTTTGTGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAggg                                             >  2:2180257/1‑100 (MQ=255)
                                                   gTATTTAATGGAAATAAAACACCAAAGAGCGTTGAGATATATATTAAA‑G‑AATAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTagag                                               <  2:2303901/98‑1 (MQ=255)
                                                     aTTTTATTGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTagag                                               <  2:2984590/96‑1 (MQ=255)
                                                     aTTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGaa                                           >  1:719186/1‑100 (MQ=255)
                                                     aTTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGaa                                           >  2:2590433/1‑100 (MQ=255)
                                                      ttttATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAAt                                          >  1:451616/1‑100 (MQ=255)
                                                      ttttATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAAt                                          >  1:2027839/1‑100 (MQ=255)
                                                      ttttATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGAAGAGGGaa                                           >  2:2126995/1‑99 (MQ=255)
                                                       tatATGGACATAAAAAACCAAAGAGGGTTGAAATATAAAAAAAT‑AAAA‑AAAAATATTATTTATAAAATAATTAATGACAAAAAAATTAGTAGAGGGAAt                                          <  2:2837389/97‑1 (MQ=255)
                                                        ttATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATgg                                        <  1:1373122/100‑1 (MQ=255)
                                                          atgaaaaTAAAAAACAAAAGAGTGTTTAAATATAAATTTAA‑AAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATgg                                        <  2:959075/92‑1 (MQ=255)
                                                          aTGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTc                                      <  2:2390042/100‑1 (MQ=255)
                                                            ggACATAAAACACCAAAGAGCGTTGAGATATAAATTTCA‑GAAA‑AAAAATATTATTTATAAAATAAAAAAAAAAAAAAAAATTAGTAGAGGGAATGGTCAt                                    <  2:2495278/100‑1 (MQ=255)
                                                            ggACATAAAACACAAAAGAGCGTTGAGATATAACTTTAT‑GAAA‑AAAAATATTATTTATAAAATAAAAAAGAAAAAAAAAATTAGTAGAgg                                              <  2:2544407/90‑1 (MQ=255)
                                                             gACATAAAACACTAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtt                                   >  1:2733425/1‑100 (MQ=255)
                                                             gACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtt                                   >  1:1386906/1‑100 (MQ=255)
                                                             gACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtt                                   >  1:726630/1‑100 (MQ=255)
                                                             gACATAAAACAACAAAGAGTGTTAAGATATAAATTTATAA‑AA‑AAAAATATTATTTATAAAATAAAAAAAAAAAAAAAAATTAGTAGAGGGAATgg                                        <  2:847227/95‑1 (MQ=255)
                                                              aCATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAttt                                  <  1:355888/100‑1 (MQ=255)
                                                              aCATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtt                                   <  2:3129222/99‑1 (MQ=255)
                                                               cATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtttt                                 >  2:632783/1‑100 (MQ=255)
                                                                aTAAAACACCAAAGAGCGTTAAGATATAAAAAAAA‑AAAA‑AAAAATATTATTTATAAGATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtttt                                 <  2:46570/99‑1 (MQ=255)
                                                                 tAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtat                               >  1:1997244/1‑100 (MQ=255)
                                                                 tAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtat                               >  2:896912/1‑100 (MQ=255)
                                                                  aaaaCACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtata                              >  2:2395989/1‑100 (MQ=255)
                                                                  aaaaCACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtata                              >  2:3473241/1‑100 (MQ=255)
                                                                  aaaaCACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtata                              >  1:2249465/1‑100 (MQ=255)
                                                                  aaaaCACCAAAGAGCGTTAAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtata                              <  2:3027319/100‑1 (MQ=255)
                                                                    aaCACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAg                            <  1:2024873/100‑1 (MQ=255)
                                                                    aaCACCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAg                            <  1:3217234/100‑1 (MQ=255)
                                                                    aaCACCAAAGAGCGTCGAGATATATATTTAT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAg                            <  2:3972368/100‑1 (MQ=255)
                                                                     acacCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGt                           >  2:3915425/1‑100 (MQ=255)
                                                                      cacCAAAGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGta                          <  2:3143996/100‑1 (MQ=255)
                                                                          aaaGAGCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAAt                      >  2:1145950/1‑100 (MQ=255)
                                                                           atctgCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATg                     <  2:2785202/96‑1 (MQ=255)
                                                                             gagCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGt                   <  2:4079539/100‑1 (MQ=255)
                                                                               gCGTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATgagat                 >  1:2785218/1‑96 (MQ=255)
                                                                                 gTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTATTGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaat               >  2:350116/1‑100 (MQ=255)
                                                                                 gTTGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaat               <  1:3866061/100‑1 (MQ=255)
                                                                                  ttGAGATATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGTGAATGGTTATTTTATAAGTATAATGGTtaata              >  2:1905569/1‑100 (MQ=255)
                                                                                    gagaTATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGAACAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaata              <  2:1146365/98‑1 (MQ=255)
                                                                                     agaTATATCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAAtt           >  1:1916970/1‑100 (MQ=255)
                                                                                     agaTATATCTTTCT‑G‑AA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATTGTCATTTTATAAGTATATTGGTTaataatca          <  1:2808223/100‑3 (MQ=255)
                                                                                       atatatCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTAc         <  2:1052273/100‑1 (MQ=255)
                                                                                        tatatCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACa        >  1:2162148/1‑100 (MQ=255)
                                                                                         atatTTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAg       <  2:2242070/100‑1 (MQ=255)
                                                                                         atatCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACa        <  1:4038512/99‑1 (MQ=255)
                                                                                         atatCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAg       <  1:955248/100‑1 (MQ=255)
                                                                                         atatCTTTCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAg       <  1:1237477/100‑1 (MQ=255)
                                                                                          tatCTTTCT‑GAAA‑AAAAATATTATTTATAAAAT‑ATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGa     >  2:2796653/1‑100 (MQ=255)
                                                                                              tttCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGCAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGc  >  1:56337/1‑100 (MQ=255)
                                                                                               ttCT‑GAAA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTATTAGTAGGAATTGTTATTTTAAAAGTATAATAGTTAATAATTACAGGACg   >  2:3653887/1‑98 (MQ=255)
                                                                                                     |                                                                                              
TGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑G‑AA‑AAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGC  >  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: ‑