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.A93.6% 355.0 / 16.9 109coding (148/483 nt)ydcDputative immunity protein for RhsE
Reads supporting (aligned to +/- strand):  ref base . (4/3);  new base A (52/50);  total (56/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 123 total aligned reads displayed.

TGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGC  >  NC_000913/1529970‑1530162
                                                                                                     |                                                                                               
tGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTGTTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTg                                                                                                  <  2:1418266/100‑1 (MQ=255)
tGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTg                                                                                                  >  2:2429154/1‑100 (MQ=255)
tGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTg                                                                                                  <  1:606764/100‑1 (MQ=255)
 gTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑a                                                                                               <  2:111113/100‑1 (MQ=255)
 gTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑a                                                                                               <  1:1260427/100‑1 (MQ=255)
  tttATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑aa                                                                                              >  1:1518616/1‑100 (MQ=255)
    tATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑aaaa                                                                                            >  1:1078774/1‑100 (MQ=255)
    tATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑aaaa                                                                                            >  1:1078846/1‑100 (MQ=255)
        tgttgtTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑aaaaaaaa                                                                                       >  1:1061786/1‑99 (MQ=255)
         gttgttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑aaaaaaaat                                                                                      >  2:2097848/1‑98 (MQ=255)
            ntttACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑aaaaaaaatatt                                                                                   >  2:2159388/2‑95 (MQ=255)
            gttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑aaaaaaaatatt                                                                                   >  2:2554929/1‑95 (MQ=255)
             tttACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAAtatta                                                                                   <  2:1295218/100‑1 (MQ=255)
             tttACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAAtatta                                                                                   <  1:1629781/100‑1 (MQ=255)
               tACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATattatt                                                                                 >  1:682806/1‑100 (MQ=255)
               tACCTTGGGGTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGAGTTTAGATATATCTTTCTGAAAAAAAATattatt                                                                                 >  2:1113000/1‑100 (MQ=255)
                aCCTTGGGTTTGATAACAATGCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATAttattt                                                                                <  2:545139/100‑1 (MQ=255)
                 ccTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATAttattta                                                                               <  1:282832/100‑1 (MQ=255)
                 ccTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATAttattta                                                                               >  1:1907638/1‑100 (MQ=255)
                   ttGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGGGATATATCTTTCTGAAAAAAAATAatatttat                                                                              >  2:629685/1‑99 (MQ=255)
                    tGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATaa                                                                            <  1:2256350/100‑1 (MQ=255)
                     gggTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATaaa                                                                           >  1:1524075/1‑100 (MQ=255)
                        tttGATAACAATTCAAATTCTGACATCGTATATTATGGACATAAAACACCAAAGAGCGTTGATATATATCTTTCTG‑‑AAAAAAATATTATTTATaaaa                                                                          <  2:1433118/97‑1 (MQ=255)
                         ttaaTAACAAATCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAAAAAAAAGAAAAAAAATATTATTTATAAaata                                                                        <  2:2434482/96‑1 (MQ=255)
                         ttGATAACAATTGAAATTATGACATTGTGTTTTATGGACATAAAAAACAAAACAGATTTGAAATAAATATTTCAGAAAAAAAATATTATTTATaaaa                                                                          <  1:1231658/97‑1 (MQ=255)
                         ttGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAaataa                                                                       >  1:2188018/1‑100 (MQ=255)
                          tGATAACAATTCAAATTCTGACATCGAATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAaataat                                                                      >  1:1561980/1‑100 (MQ=255)
                            aTAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAAtaatta                                                                    >  1:717630/1‑100 (MQ=255)
                            aTAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAAtaatta                                                                    >  2:458917/1‑100 (MQ=255)
                            aTAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAAtaatta                                                                    >  2:1637844/1‑100 (MQ=255)
                             tAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAAtaattaa                                                                   >  1:1657000/1‑100 (MQ=255)
                             tAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAAtaattaa                                                                   >  1:2060668/1‑100 (MQ=255)
                              aaCAAGTCAAATTCAGATATCGTATTTTATGAACATAAAACATCAAAGAACTTTGAGATATATCTTTCTG‑‑AAAAAAATATTATTTATAAAATAATTAATg                                                                 <  2:2157379/100‑1 (MQ=255)
                                 aaTTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGAc                                                               <  1:1986815/100‑1 (MQ=255)
                                  aTTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGAcc                                                              >  1:1135890/1‑100 (MQ=255)
                                   ttCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCa                                                             >  1:1610509/1‑100 (MQ=255)
                                   ttCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCa                                                             >  2:1193114/1‑100 (MQ=255)
                                     cAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCaaa                                                           >  1:181584/1‑100 (MQ=255)
                                     cAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCaaa                                                           >  1:1334260/1‑100 (MQ=255)
                                     cAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCCTTATG‑‑AAAACAAATATTATTTATAAAATAATTAATGACCaaa                                                           >  2:1798696/1‑100 (MQ=255)
                                        aTTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCaaaaaa                                                        >  2:663323/1‑100 (MQ=255)
                                         ttCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAAt                                                       <  2:2060654/100‑1 (MQ=255)
                                          tCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAAtt                                                      <  1:147796/100‑1 (MQ=255)
                                            tGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtag                                                    <  1:180219/100‑1 (MQ=255)
                                            tGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtag                                                    <  1:2239155/100‑1 (MQ=255)
                                            tGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTCATAAAATAATTAATGACCAAAAAATtag                                                    <  1:1414324/100‑1 (MQ=255)
                                             gACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtagt                                                   <  1:2389043/100‑1 (MQ=255)
                                             gACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtagt                                                   >  1:2109207/1‑100 (MQ=255)
                                               cATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtagta                                                  <  2:505266/99‑1 (MQ=255)
                                               cATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtagta                                                  >  2:151813/1‑98 (MQ=255)
                                                 tCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTagag                                               <  1:135255/100‑1 (MQ=255)
                                                 tCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAAAATTATTTATAAAATAATTAATGAACAAAAAATTAGTagag                                               >  2:2315547/1‑100 (MQ=255)
                                                  cGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAgg                                              >  2:151168/1‑100 (MQ=255)
                                                   gTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTagagag                                             >  1:1322554/1‑98 (MQ=255)
                                                   gTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAggg                                             >  1:407594/1‑100 (MQ=255)
                                                    tATTTTATGGAGATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGa                                            >  2:1967054/1‑100 (MQ=255)
                                                       tttATGGACATAAAACACAAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATg                                         <  2:2109193/100‑1 (MQ=255)
                                                         tATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGt                                       <  1:1112992/100‑1 (MQ=255)
                                                         tATGGACATAAAACACCAAAGAGCGTTGAGATATAAATTTCTGAAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATgg                                        <  2:1078766/100‑1 (MQ=255)
                                                         tATGGACATAAAACACCAAAGAGCATTGAGATATATCTTTCTGAAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATgg                                        <  2:1078838/100‑1 (MQ=255)
                                                             gACATAAAACACCAAAGAGCTTTGAGATATATCTATATGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtt                                   <  2:546498/100‑1 (MQ=255)
                                                              aCATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAttt                                  <  1:782889/100‑1 (MQ=255)
                                                              aCATAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAttt                                  <  1:2390335/100‑1 (MQ=255)
                                                                aTAAAACAGCAAAGCGAGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtttt                                 <  2:2225775/99‑1 (MQ=255)
                                                                aTAAAACACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTta                                <  1:1165322/100‑1 (MQ=255)
                                                                  aaaaCACCAAAGAGCGTTGATATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtata                              <  2:346781/100‑1 (MQ=255)
                                                                  aaaaCACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATTTTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtata                              <  1:1116755/100‑1 (MQ=255)
                                                                  aaaaCACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtata                              >  1:564361/1‑100 (MQ=255)
                                                                  aaaaCACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtata                              <  2:1518605/100‑1 (MQ=255)
                                                                  aaaaCACCAAAGAGCGTTGAGATATATATTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtata                              <  2:20053/100‑1 (MQ=255)
                                                                  aaaaCAACAAATATCGTTGAGATATATATTTCTAA‑‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtata                              <  2:347139/100‑1 (MQ=255)
                                                                   aaaCACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATaa                             >  1:2228962/1‑100 (MQ=255)
                                                                   aaaCACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATaa                             <  1:2310765/100‑1 (MQ=255)
                                                                    aaCACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAg                            <  1:852340/100‑1 (MQ=255)
                                                                    aaCACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAg                            <  2:424773/100‑1 (MQ=255)
                                                                    aaCACCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAg                            <  2:2381076/100‑1 (MQ=255)
                                                                     acacCAAAGAGCGTTGAGATATATAATTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGt                           <  1:863377/100‑1 (MQ=255)
                                                                      cacCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGta                          >  2:1206207/1‑100 (MQ=255)
                                                                      cacCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGta                          <  1:2481213/100‑1 (MQ=255)
                                                                      cacCAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGta                          >  1:52926/1‑100 (MQ=255)
                                                                        ccAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGtata                        >  2:5056/1‑100 (MQ=255)
                                                                        ccAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGtata                        >  1:2302757/1‑100 (MQ=255)
                                                                         cAAAGAGCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATaa                       >  2:334228/1‑100 (MQ=255)
                                                                            agagCGTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATgg                    >  2:1784200/1‑100 (MQ=255)
                                                                                 gTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaat               <  2:1480412/100‑1 (MQ=255)
                                                                                 gTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaat               >  1:2108549/1‑100 (MQ=255)
                                                                                 gTTGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaat               >  2:330423/1‑100 (MQ=255)
                                                                                 gTTGAGATATATCATTATGAAAAAAAATCTTATTTATAAAATAATTAATGACCAAAAAATTAATATAGGGAATGATAATTTTATAAATATAATGGTtaat               >  2:1056612/1‑100 (MQ=255)
                                                                                   tGAGATATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaataa             <  1:328378/100‑1 (MQ=255)
                                                                                     agaTATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAAtt           <  1:2159402/100‑1 (MQ=255)
                                                                                     agaTATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAAtt           <  1:1905657/100‑1 (MQ=255)
                                                                                      gaTATATCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTa          <  1:281085/100‑1 (MQ=255)
                                                                                       atatatCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTa          <  2:1773622/99‑1 (MQ=255)
                                                                                       atatatCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTAc         >  1:258323/1‑100 (MQ=255)
                                                                                        tatatCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACa        >  1:437252/1‑100 (MQ=255)
                                                                                            tCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGAc    <  1:1873812/100‑1 (MQ=255)
                                                                                            tCTTTCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGAc    >  1:1638600/1‑100 (MQ=255)
                                                                                              tttCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGc  <  2:258320/100‑1 (MQ=255)
                                                                                              tttCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGc  >  1:375657/1‑100 (MQ=255)
                                                                                              tttCTGAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGc  >  1:1458145/1‑100 (MQ=255)
                                                                                                     |                                                                                               
TGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCTG‑‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGC  >  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: