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.

TGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑G‑‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGC  >  NC_000913/1529970‑1530162
                                                                                                     |                                                                                               
tGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTGTTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑g                                                                                                  <  2:1418266/100‑1 (MQ=255)
tGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑g                                                                                                  >  2:2429154/1‑100 (MQ=255)
tGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑g                                                                                                  <  1:606764/100‑1 (MQ=255)
 gTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑G‑‑a                                                                                               <  2:111113/100‑1 (MQ=255)
 gTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑G‑‑a                                                                                               <  1:1260427/100‑1 (MQ=255)
  tttATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑G‑‑aa                                                                                              >  1:1518616/1‑100 (MQ=255)
    tATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑G‑‑aaaa                                                                                            >  1:1078774/1‑100 (MQ=255)
    tATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑G‑‑aaaa                                                                                            >  1:1078846/1‑100 (MQ=255)
        tgttgtTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑G‑‑aaaa‑aaaa                                                                                       >  1:1061786/1‑99 (MQ=255)
         gttgttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑G‑‑aaaa‑aaaat                                                                                      >  2:2097848/1‑98 (MQ=255)
            ntttACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑G‑‑aaaa‑aaaatatt                                                                                   >  2:2159388/2‑95 (MQ=255)
            gttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑G‑‑aaaa‑aaaatatt                                                                                   >  2:2554929/1‑95 (MQ=255)
             tttACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAAtatta                                                                                   <  2:1295218/100‑1 (MQ=255)
             tttACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAAtatta                                                                                   <  1:1629781/100‑1 (MQ=255)
               tACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATattatt                                                                                 >  1:682806/1‑100 (MQ=255)
               tACCTTGGGGTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGAGTTTAGATATATCTTTCT‑GA‑AAAA‑AAATattatt                                                                                 >  2:1113000/1‑100 (MQ=255)
                aCCTTGGGTTTGATAACAATGCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATAttattt                                                                                <  2:545139/100‑1 (MQ=255)
                 ccTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATAttattta                                                                               <  1:282832/100‑1 (MQ=255)
                 ccTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATAttattta                                                                               >  1:1907638/1‑100 (MQ=255)
                   ttGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGGGATATATCTTTCT‑GA‑AAAA‑AAATAatatttat                                                                              >  2:629685/1‑99 (MQ=255)
                    tGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATaa                                                                            <  1:2256350/100‑1 (MQ=255)
                     gggTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATaaa                                                                           >  1:1524075/1‑100 (MQ=255)
                        tttGATAACAATTCAAATTCTGACATCGTATATTATGGACATAAAACACCAAAGAGCGTTGATATATATCTTTCT‑G‑‑AAAA‑AAATATTATTTATaaaa                                                                          <  2:1433118/97‑1 (MQ=255)
                         ttaaTAACAAATCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAAAAAAAA‑GA‑AAAA‑AAATATTATTTATAAaata                                                                        <  2:2434482/96‑1 (MQ=255)
                         ttGATAACAATTGAAATTATGACATTGTGTTTTATGGACATAAAAAACAAAACAGATTTGAAATAAATATTTCA‑GA‑AAAA‑AAATATTATTTATaaaa                                                                          <  1:1231658/97‑1 (MQ=255)
                         ttGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAaataa                                                                       >  1:2188018/1‑100 (MQ=255)
                          tGATAACAATTCAAATTCTGACATCGAATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAaataat                                                                      >  1:1561980/1‑100 (MQ=255)
                            aTAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAAtaatta                                                                    >  1:717630/1‑100 (MQ=255)
                            aTAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAAtaatta                                                                    >  2:458917/1‑100 (MQ=255)
                            aTAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAAtaatta                                                                    >  2:1637844/1‑100 (MQ=255)
                             tAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAAtaattaa                                                                   >  1:1657000/1‑100 (MQ=255)
                             tAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAAtaattaa                                                                   >  1:2060668/1‑100 (MQ=255)
                              aaCAAGTCAAATTCAGATATCGTATTTTATGAACATAAAACATCAAAGAACTTTGAGATATATCTTTCT‑G‑‑AAAA‑AAATATTATTTATAAAATAATTAATg                                                                 <  2:2157379/100‑1 (MQ=255)
                                 aaTTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGAc                                                               <  1:1986815/100‑1 (MQ=255)
                                  aTTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGAcc                                                              >  1:1135890/1‑100 (MQ=255)
                                   ttCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCa                                                             >  1:1610509/1‑100 (MQ=255)
                                   ttCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCa                                                             >  2:1193114/1‑100 (MQ=255)
                                     cAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCaaa                                                           >  1:181584/1‑100 (MQ=255)
                                     cAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCaaa                                                           >  1:1334260/1‑100 (MQ=255)
                                     cAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCCTTAT‑G‑‑AAAACAAATATTATTTATAAAATAATTAATGACCaaa                                                           >  2:1798696/1‑100 (MQ=255)
                                        aTTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCaaaaaa                                                        >  2:663323/1‑100 (MQ=255)
                                         ttCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAAt                                                       <  2:2060654/100‑1 (MQ=255)
                                          tCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAAtt                                                      <  1:147796/100‑1 (MQ=255)
                                            tGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATtag                                                    <  1:180219/100‑1 (MQ=255)
                                            tGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATtag                                                    <  1:2239155/100‑1 (MQ=255)
                                            tGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTCATAAAATAATTAATGACCAAAAAATtag                                                    <  1:1414324/100‑1 (MQ=255)
                                             gACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATtagt                                                   <  1:2389043/100‑1 (MQ=255)
                                             gACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATtagt                                                   >  1:2109207/1‑100 (MQ=255)
                                               cATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATtagta                                                  <  2:505266/99‑1 (MQ=255)
                                               cATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑G‑‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATtagta                                                  >  2:151813/1‑98 (MQ=255)
                                                 tCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTagag                                               <  1:135255/100‑1 (MQ=255)
                                                 tCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAAAATTATTTATAAAATAATTAATGAACAAAAAATTAGTagag                                               >  2:2315547/1‑100 (MQ=255)
                                                  cGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAgg                                              >  2:151168/1‑100 (MQ=255)
                                                   gTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTagagag                                             >  1:1322554/1‑98 (MQ=255)
                                                   gTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAggg                                             >  1:407594/1‑100 (MQ=255)
                                                    tATTTTATGGAGATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGa                                            >  2:1967054/1‑100 (MQ=255)
                                                       tttATGGACATAAAACACAAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATg                                         <  2:2109193/100‑1 (MQ=255)
                                                         tATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGt                                       <  1:1112992/100‑1 (MQ=255)
                                                         tATGGACATAAAACACCAAAGAGCGTTGAGATATAAATTTCT‑GAAAAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATgg                                        <  2:1078766/100‑1 (MQ=255)
                                                         tATGGACATAAAACACCAAAGAGCATTGAGATATATCTTTCT‑GAAAAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATgg                                        <  2:1078838/100‑1 (MQ=255)
                                                             gACATAAAACACCAAAGAGCTTTGAGATATATCTATAT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtt                                   <  2:546498/100‑1 (MQ=255)
                                                              aCATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAttt                                  <  1:782889/100‑1 (MQ=255)
                                                              aCATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAttt                                  <  1:2390335/100‑1 (MQ=255)
                                                                aTAAAACAGCAAAGCGAGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtttt                                 <  2:2225775/99‑1 (MQ=255)
                                                                aTAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTta                                <  1:1165322/100‑1 (MQ=255)
                                                                  aaaaCACCAAAGAGCGTTGATATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtata                              <  2:346781/100‑1 (MQ=255)
                                                                  aaaaCACCAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATTTTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtata                              <  1:1116755/100‑1 (MQ=255)
                                                                  aaaaCACCAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtata                              >  1:564361/1‑100 (MQ=255)
                                                                  aaaaCACCAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtata                              <  2:1518605/100‑1 (MQ=255)
                                                                  aaaaCACCAAAGAGCGTTGAGATATATATTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtata                              <  2:20053/100‑1 (MQ=255)
                                                                  aaaaCAACAAATATCGTTGAGATATATATTTCTAA‑‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtata                              <  2:347139/100‑1 (MQ=255)
                                                                   aaaCACCAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATaa                             >  1:2228962/1‑100 (MQ=255)
                                                                   aaaCACCAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATaa                             <  1:2310765/100‑1 (MQ=255)
                                                                    aaCACCAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAg                            <  1:852340/100‑1 (MQ=255)
                                                                    aaCACCAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAg                            <  2:424773/100‑1 (MQ=255)
                                                                    aaCACCAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAg                            <  2:2381076/100‑1 (MQ=255)
                                                                     acacCAAAGAGCGTTGAGATATATAATTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGt                           <  1:863377/100‑1 (MQ=255)
                                                                      cacCAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGta                          >  2:1206207/1‑100 (MQ=255)
                                                                      cacCAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGta                          <  1:2481213/100‑1 (MQ=255)
                                                                      cacCAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGta                          >  1:52926/1‑100 (MQ=255)
                                                                        ccAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGtata                        >  2:5056/1‑100 (MQ=255)
                                                                        ccAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGtata                        >  1:2302757/1‑100 (MQ=255)
                                                                         cAAAGAGCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATaa                       >  2:334228/1‑100 (MQ=255)
                                                                            agagCGTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATgg                    >  2:1784200/1‑100 (MQ=255)
                                                                                 gTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaat               <  2:1480412/100‑1 (MQ=255)
                                                                                 gTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaat               >  1:2108549/1‑100 (MQ=255)
                                                                                 gTTGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaat               >  2:330423/1‑100 (MQ=255)
                                                                                 gTTGAGATATATCATTAT‑GA‑AAAA‑AAATCTTATTTATAAAATAATTAATGACCAAAAAATTAATATAGGGAATGATAATTTTATAAATATAATGGTtaat               >  2:1056612/1‑100 (MQ=255)
                                                                                   tGAGATATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaataa             <  1:328378/100‑1 (MQ=255)
                                                                                     agaTATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAAtt           <  1:2159402/100‑1 (MQ=255)
                                                                                     agaTATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAAtt           <  1:1905657/100‑1 (MQ=255)
                                                                                      gaTATATCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTa          <  1:281085/100‑1 (MQ=255)
                                                                                       atatatCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTa          <  2:1773622/99‑1 (MQ=255)
                                                                                       atatatCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTAc         >  1:258323/1‑100 (MQ=255)
                                                                                        tatatCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACa        >  1:437252/1‑100 (MQ=255)
                                                                                            tCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGAc    <  1:1873812/100‑1 (MQ=255)
                                                                                            tCTTTCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGAc    >  1:1638600/1‑100 (MQ=255)
                                                                                              tttCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGc  <  2:258320/100‑1 (MQ=255)
                                                                                              tttCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGc  >  1:375657/1‑100 (MQ=255)
                                                                                              tttCT‑GA‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGc  >  1:1458145/1‑100 (MQ=255)
                                                                                                     |                                                                                               
TGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATATCTTTCT‑G‑‑AAAA‑AAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGC  >  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: ‑