Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,530,076 (A)7→8 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.A96.0% 504.8 / 11.6 151coding (148/483 nt)ydcDputative immunity protein for RhsE
Reads supporting (aligned to +/- strand):  ref base . (1/5);  new base A (77/68);  total (78/73)
Fisher's exact test for biased strand distribution p-value = 1.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 167 total aligned reads displayed.

TGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGC  >  NC_000913/1529970‑1530162
                                                                                                      |                                                                                              
tGTTTATTTGTTGTTTACCTTGTGTTTGATATCAATTGAAATTCTGATATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTg                                                                                                 <  1:2349719/100‑1 (MQ=255)
tGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGCCATAAAACACCAAAGAGCGTTGAGATATATCTTTCTg                                                                                                 <  1:1852941/100‑1 (MQ=255)
tGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTg                                                                                                 <  2:565590/100‑1 (MQ=255)
   ttATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGaaa                                                                                             >  1:2699088/1‑100 (MQ=255)
   ttATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGaaa                                                                                             <  1:2361182/100‑1 (MQ=255)
    tgtttGTTGTTTAGCTTGGGTTTGATAACAAATCAAATTCTGGCATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGaaaa                                                                                            >  2:2343211/3‑100 (MQ=255)
     aTTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGaaaaa                                                                                           >  1:1688960/1‑100 (MQ=255)
       ttgttgTTTACCTTGGGTTTGATAACGATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGaaaaaaa                                                                                        >  2:2118884/1‑100 (MQ=255)
         gttgttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGaaaaaaaat                                                                                      <  2:94198/100‑3 (MQ=255)
            gttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGaaaaaaaatatt                                                                                   >  2:1065514/1‑95 (MQ=255)
              ttACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAAtattat                                                                                  <  1:503412/100‑1 (MQ=255)
              ttACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGAACATAAAACACCAAAGAGCGTTGATATATATCTTTCTGAAAAAAAAtattat                                                                                  <  1:963126/100‑1 (MQ=255)
               tACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATattatt                                                                                 <  1:2936474/100‑1 (MQ=255)
               tACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATattatt                                                                                 <  2:1106472/100‑1 (MQ=255)
                aCCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATAttattt                                                                                >  2:726955/1‑100 (MQ=255)
                 ccTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATAttattta                                                                               >  2:1363169/1‑100 (MQ=255)
                 ccTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATAttattta                                                                               >  2:54128/1‑100 (MQ=255)
                  cTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAATATTATTtata                                                                             <  2:2134714/100‑1 (MQ=255)
                   ttGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATAttattt                                                                                <  1:2683721/97‑1 (MQ=255)
                     gggTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATaaa                                                                           >  1:1262650/1‑100 (MQ=255)
                      ggTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATaaaa                                                                          >  1:452363/1‑100 (MQ=255)
                      ggTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATaaaa                                                                          >  2:922994/1‑100 (MQ=255)
                        tttGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAaata                                                                        >  2:1902072/1‑100 (MQ=255)
                          tGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAaataat                                                                      <  1:726946/100‑1 (MQ=255)
                           gATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAAtaatt                                                                     >  2:2405116/1‑100 (MQ=255)
                             tAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAAtaattaa                                                                   >  2:964378/1‑100 (MQ=255)
                              aaCAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAAtaattaat                                                                  >  2:23527/1‑100 (MQ=255)
                               aaaaatCAAATTCTGACATCTTATTTTATAAACATAAAACAAAAAAGATCTTTGAAATAAAAAAAAAAAAAAAAAAATATTATTTATAAAATAATTAATg                                                                 <  1:1466025/95‑1 (MQ=255)
                                 aaTTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTCGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGAc                                                               >  1:1453617/1‑100 (MQ=255)
                                     cAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCaaa                                                           >  2:896827/1‑100 (MQ=255)
                                       aaTTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCaaaaa                                                         >  1:1887809/1‑100 (MQ=255)
                                        aTTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCaaaaaa                                                        >  2:1278551/1‑100 (MQ=255)
                                        aTTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCaaaaaa                                                        >  1:65896/1‑100 (MQ=255)
                                          tCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAATTAATTAATGACCAAAAAAtt                                                      <  1:2118853/100‑1 (MQ=255)
                                          tCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAAtt                                                      <  1:279002/100‑1 (MQ=255)
                                           cTGAAAACGTATTTTATAAACATAAAACACCAAAGATCGTTTAGATAAAAAAAAAAAAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATta                                                     <  1:1585934/100‑1 (MQ=255)
                                            tGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtag                                                    <  1:1696957/100‑1 (MQ=255)
                                            tGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtag                                                    <  2:398724/100‑1 (MQ=255)
                                            tGAAATCGTATTTTATAGACATAAAACACAAAAGAGCGTTGAGATATATATATCTAAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtag                                                    <  1:1120635/100‑1 (MQ=255)
                                            ctctttggtgttttATGTCCATAAAATACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGAACAAAAAATtag                                                    >  2:3289352/11‑100 (MQ=25)
                                             gACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtagt                                                   <  2:855560/100‑1 (MQ=255)
                                             gACATATTATTTTATGGAAAAAAAACAACAAATAGCTTTGAGATATAACTTTCTAAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtag                                                    <  1:679345/99‑1 (MQ=255)
                                              aaataatataTTAAGGACATAAAACAACAAATAGCGTTGAAATATAACATACTGAAAAAAAATATTATTTATAAAATAAAAAAAAAAAAAAAAATtag                                                    <  1:347805/92‑1 (MQ=255)
                                                aTCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAAtcaat                                                   >  1:3167571/1‑93 (MQ=255)
                                                aTCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTaga                                                <  1:3157346/100‑1 (MQ=255)
                                                  ggtgttttATGTCCATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAgg                                              >  2:1740744/5‑100 (MQ=255)
                                                   ttATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAggg                                             <  2:3117522/99‑1 (MQ=255)
                                                   gTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAggg                                             <  2:3385436/100‑1 (MQ=255)
                                                   gTATTTTATGGAAATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAgg                                              <  1:1943022/99‑1 (MQ=255)
                                                    tATTTTATGGACATAAAACACCAAATATCGTTGAGATATAACAAACAGAAAAAAAATATTATTTATAAAATAATTACTGACCAAAAAATTAGTAGAggg                                             <  1:2017447/99‑1 (MQ=255)
                                                     aTTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGaa                                           >  1:22082/1‑100 (MQ=255)
                                                         tATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGt                                       <  1:3160259/100‑1 (MQ=255)
                                                           tGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCa                                     >  1:1968871/1‑100 (MQ=255)
                                                             gACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtt                                   <  2:2308465/100‑1 (MQ=255)
                                                             gACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtt                                   <  2:283203/100‑1 (MQ=255)
                                                             gACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAt                                    >  1:3313236/1‑99 (MQ=255)
                                                              aCATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAttt                                  >  2:456596/1‑100 (MQ=255)
                                                               cATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtttt                                 >  2:1912379/1‑100 (MQ=255)
                                                                aTAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtttt                                 <  1:964366/99‑1 (MQ=255)
                                                                 tAAAACACCAAAGGGCGTTGAGATATAACTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAGATTAGTAGAGGGAATGGTCATTTtan                               <  1:2644892/100‑2 (MQ=255)
                                                                 tAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtat                               >  2:2832154/1‑100 (MQ=255)
                                                                 tAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtat                               >  1:245404/1‑100 (MQ=255)
                                                                 tAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtat                               >  2:1272033/1‑100 (MQ=255)
                                                                 tAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtat                               >  2:1626519/1‑100 (MQ=255)
                                                                  aaaaCACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtata                              >  2:490492/1‑100 (MQ=255)
                                                                  aaaaCACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtata                              <  2:1705872/100‑1 (MQ=255)
                                                                  aaaaCACCAAAGAGCGTTAAAATATATCTTTCTAAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtata                              <  1:2489607/100‑1 (MQ=255)
                                                                   aaaCACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATaa                             >  2:1341678/1‑100 (MQ=255)
                                                                    aaaacCAAATAGCGTTGAAATATAACATAATGAAAAAAAATATTATTTATAAAATAAAAAAAAAAAAAAAAATTAGTAGAGGGAATGGTCATTTTATAAg                            <  1:2384266/97‑1 (MQ=255)
                                                                    aaCACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAg                            <  1:192690/100‑1 (MQ=255)
                                                                    aaCACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAg                            <  1:305707/100‑1 (MQ=255)
                                                                    aaCACCAAAGAGCGTTGAGATACATCTTTCTCAAAAACAATATTATTGATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAg                            <  1:2347992/100‑1 (MQ=255)
                                                                    aaCACCAAAGAGCGTAGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAg                            <  2:3417350/100‑1 (MQ=255)
                                                                       acCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGtat                         >  1:2855499/1‑100 (MQ=255)
                                                                           aaGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATg                     <  1:1765698/100‑1 (MQ=255)
                                                                            agagCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTTATTTTATAAGTATAATgg                    >  1:2425624/1‑100 (MQ=255)
                                                                            agagCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATgg                    <  2:3056391/100‑1 (MQ=255)
                                                                            agagCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATg                     >  1:336122/1‑99 (MQ=255)
                                                                                cGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaa                >  1:106147/1‑100 (MQ=255)
                                                                                 gTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATAAGTAGTGGGAATGGTCATTTTATAAGTATAATGGTtaa                >  1:3106659/1‑99 (MQ=255)
                                                                                  ttGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaata              >  2:1800552/1‑100 (MQ=255)
                                                                                   taagaTATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaat               <  1:465173/96‑1 (MQ=255)
                                                                                   tGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaataa             <  2:719168/100‑1 (MQ=255)
                                                                                      gaTATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTa          <  2:1936158/100‑1 (MQ=255)
                                                                                      gaTATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTa          >  1:2836985/1‑100 (MQ=255)
                                                                                      gaTATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTa          >  2:1087497/1‑100 (MQ=255)
                                                                                       atatatCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAAtt           <  1:687397/98‑1 (MQ=255)
                                                                                       atatatCTTTCGGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTa          <  2:3139354/99‑1 (MQ=255)
                                                                                        tatatCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATTACCAAAAAATTATTATAGGGAATGGTTATTTTATAAGTATAATGGTTAATAATTAc         >  1:333238/1‑99 (MQ=255)
                                                                                        tatatCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACa        >  1:772464/1‑100 (MQ=255)
                                                                                        tatatCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACa        <  2:452368/100‑1 (MQ=255)
                                                                                         atatCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAg       <  2:106148/100‑1 (MQ=255)
                                                                                         ataaCTATCTGAAAAAAAATATTATTTATAAAATAAAAAATAAAAAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAg       <  1:305634/100‑1 (MQ=255)
                                                                                           atCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGa     >  1:924370/1‑100 (MQ=255)
                                                                                            tCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGACTGGTCATTTTATAAGTATAATGGTTAATAATTACAGGAc    >  2:667223/1‑100 (MQ=255)
                                                                                              cTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGATCATATTATAAGTATAATGGTTAATAATTACAGGAc    >  1:3415688/1‑99 (MQ=255)
                                                                                               tttCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAGTGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGc  <  1:1907001/100‑1 (MQ=255)
                                                                                               tttCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGc  <  2:1215377/100‑1 (MQ=255)
                                                                                               tttCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGCGAATGATCATTTTACAAGAATAATGGTTAATAATTACAGAACGc  >  1:36242/1‑100 (MQ=255)
                                                                                                ttCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGc  <  1:236809/99‑1 (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: