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.A94.3% 287.9 / 10.3 88coding (148/483 nt)ydcDputative immunity protein for RhsE
Reads supporting (aligned to +/- strand):  ref base . (1/4);  new base A (46/37);  total (47/41)
Fisher's exact test for biased strand distribution p-value = 1.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 105 total aligned reads displayed.

TGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTG‑‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTG  >  NC_000913/1529970‑1530168
                                                                                                     |                                                                                                    
tGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTg                                                                                                       >  1:2898673/1‑100 (MQ=255)
tGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTg                                                                                                       <  2:3538763/100‑1 (MQ=255)
 gTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTG‑‑a                                                                                                    >  1:2557395/1‑100 (MQ=255)
  tttATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTG‑‑aa                                                                                                   >  2:2574791/1‑100 (MQ=255)
   ttATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTG‑‑aaa                                                                                                  <  2:2239310/100‑1 (MQ=255)
   ttATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTG‑‑aaa                                                                                                  <  1:3269195/100‑1 (MQ=255)
       ttgttgTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTG‑‑aaaaaaa                                                                                              >  1:1563004/1‑100 (MQ=255)
       ttgttgTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTG‑‑aaaaaaa                                                                                              <  1:1760446/100‑1 (MQ=255)
       ttgttgTTTACCTTGAGCTTGATCACAATTCAAATTCTGACATCGTATTTTATGGACATAACACACCAAAGAGCGCTGAGTTATAT‑CTGTCTG‑‑a                                                                                                    >  2:3886559/1‑94 (MQ=255)
         gttgttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTG‑‑aaaaaaaat                                                                                            <  2:3181885/100‑3 (MQ=255)
         gttgttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTG‑‑aaaaaaaat                                                                                            <  1:1724732/100‑3 (MQ=255)
            gttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTG‑‑aaaaaaaatatt                                                                                         >  1:1926016/1‑95 (MQ=255)
            gttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTG‑‑aaaaaaaatatt                                                                                         >  2:1417862/1‑95 (MQ=255)
            gttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTG‑‑aaaaaaaatatt                                                                                         >  1:4067813/1‑95 (MQ=255)
            gttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTG‑‑aaaaaaaatatt                                                                                         >  2:3804941/1‑95 (MQ=255)
            gttTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTG‑‑aaaaaaaatatt                                                                                         >  1:1597877/1‑95 (MQ=255)
             tttACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAAtatta                                                                                         >  1:3319005/1‑100 (MQ=255)
               tACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATattatt                                                                                       >  2:1621579/1‑100 (MQ=255)
               tACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTG‑‑AAAAAAATAttattt                                                                                      <  1:2115682/100‑1 (MQ=255)
                aCCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATAttattt                                                                                      <  1:234553/100‑1 (MQ=255)
                aCCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATAttattt                                                                                      <  1:234604/100‑1 (MQ=255)
                aCCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATActattt                                                                                      >  2:865115/1‑100 (MQ=255)
                 ccTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATAttattta                                                                                     >  1:3365725/1‑100 (MQ=255)
                  cTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATAttatttat                                                                                    >  1:741496/1‑100 (MQ=255)
                     gggTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATaaa                                                                                 >  2:3999944/1‑100 (MQ=255)
                      ggTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATaaaa                                                                                <  1:789611/100‑1 (MQ=255)
                      ggTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATaaaa                                                                                >  1:2210521/1‑100 (MQ=255)
                       gTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAaat                                                                               <  1:3244597/100‑1 (MQ=255)
                       gTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAaat                                                                               >  1:2349301/1‑100 (MQ=255)
                       gTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAaat                                                                               >  2:3736296/1‑100 (MQ=255)
                          tGATAACAATTCAAATTCTGACATCGTATTTTGTGGACAGAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAaataat                                                                            <  2:2905722/100‑1 (MQ=255)
                          tGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAaataat                                                                            >  1:2831115/1‑100 (MQ=255)
                          tGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAaataat                                                                            >  1:1940178/1‑100 (MQ=255)
                            aTAACAAATCAAAATCTGACATCGTATTTTATGGAAATAAAACACCAAAGAGCGTTAAAATATAT‑CATTATGA‑AAAAAAATATTATTTATAAAAtaatt                                                                           <  2:2129090/99‑1 (MQ=255)
                             tAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAAAtaattaa                                                                         >  1:3068359/1‑100 (MQ=255)
                              aaCAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAAAtaattaa                                                                         >  2:206316/1‑99 (MQ=255)
                               aCAATTCAAATTCTGACATCGTATTTTATGGACATAAAACAACAAAGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAAATAATTAATg                                                                       >  2:3652034/1‑100 (MQ=255)
                               aCAATTCAAATTCTGACATCGTATTTTATCGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAAATAATTAATg                                                                       <  2:1252245/100‑1 (MQ=255)
                                 aaatCAAATTGTGACATCGCATTTTAAAGAAATAAAACAATAAAGATCGTAGAGATAAAA‑AAAAAAAA‑AAAAAAATATTATTTATAAAATAATTAATGAc                                                                     <  2:1097752/97‑1 (MQ=255)
                                 aaTTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAAATAATTAATGAc                                                                     <  2:20592/100‑1 (MQ=255)
                                  aTTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAAATAATTAATGAcc                                                                    <  2:2912235/100‑1 (MQ=255)
                                    tCAAATTCTGACATCGTATTTTATGGACATAAAACACCACAGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAAATAATTAATGACCaa                                                                  >  1:1231929/1‑100 (MQ=255)
                                     cAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAAATAATTAATGACCaaa                                                                 >  2:1515828/1‑100 (MQ=255)
                                     cAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTG‑‑AAAAAAATATTATTTATAAAATAATTAATGACCaaaa                                                                >  1:659714/1‑100 (MQ=255)
                                      aaaTTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAAATAATTAATGACCaaaa                                                                >  2:4024626/1‑100 (MQ=255)
                                       aaTTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAAATAATTAATGACCaaaaa                                                               <  1:1521793/100‑1 (MQ=255)
                                       aaTTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAAATAATTAATGACCaaaaa                                                               >  2:2073155/1‑100 (MQ=255)
                                           cTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAG‑CTTTCTGAAAAAAAAATATTATTTATAAAATAATTAATGACCAAAAAAtt                                                            >  1:2873170/1‑100 (MQ=255)
                                            tGACATCGTATTTTATGGACATAAAAAACCAAAGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtag                                                          <  2:4103716/100‑1 (MQ=255)
                                             gACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtag                                                          <  2:1748710/99‑1 (MQ=255)
                                             gACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATta                                                           >  2:2348869/1‑98 (MQ=255)
                                              aaaaaatATTTTATGGACATAAAACACAAAAGAGCGTTGAAATATAA‑ATTTCTAA‑AAAAAAATATTATTTATAAAATAAATAATAAAAAAAAAATtagta                                                        <  2:4054611/94‑1 (MQ=255)
                                              aCATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATtagta                                                        >  1:1845538/1‑100 (MQ=255)
                                                aTCGTATTTTATGGAAATAAAACACCAAAGAGCGTTTAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTaga                                                      <  2:1120614/100‑1 (MQ=255)
                                                 tCGTATTTTATGGACATAAAACACCAAAGAGCGTTTATATATAT‑ATTTCTGA‑AAAAAAATATTATTTATAAAATAATTATTAAAAAAAAAATtagt                                                         >  2:2856557/1‑96 (MQ=255)
                                                   gTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAggg                                                   >  1:2722870/1‑100 (MQ=255)
                                                    tATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGa                                                  >  1:3368125/1‑100 (MQ=255)
                                                      ttttATGGACATAAAACACCAAAGAGCGTTGAGATATAT‑AATTCTG‑‑AAAAAAATATTATTTATAAAATAATTAATGAACAAAAAATTAGTAGAGGGAATg                                               <  2:659708/100‑1 (MQ=255)
                                                        ttATGGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATgg                                              >  1:2636188/1‑100 (MQ=255)
                                                        ttATGGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATgg                                              <  2:3289484/100‑1 (MQ=255)
                                                        ttATGGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATgg                                              >  1:2636799/1‑100 (MQ=255)
                                                          aTGGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAAt                                                >  2:2513441/1‑96 (MQ=255)
                                                            ggACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAt                                          <  1:258900/100‑1 (MQ=255)
                                                            ggACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAt                                          >  1:2836215/1‑100 (MQ=255)
                                                             gACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtt                                         >  1:448506/1‑100 (MQ=255)
                                                               cATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtttt                                       <  1:3569252/100‑1 (MQ=255)
                                                               cATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCAtttt                                       >  1:22281/1‑100 (MQ=255)
                                                                  aaaaCACCAAAGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTtata                                    >  1:3965585/1‑100 (MQ=255)
                                                                   aaaCACCAAAAATCTTTGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAACATTAGTAGAGGGAATGGTCATTTTATaa                                   <  2:1141203/100‑1 (MQ=255)
                                                                    atcacCAATGAGCGTTGAGATATATAAATAAAA‑‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAg                                  <  2:1315840/98‑1 (MQ=255)
                                                                    aaCACCAAAGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAg                                  >  1:3938334/1‑100 (MQ=255)
                                                                     acacGAAAGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAAATAATTGATTACCAAAAAATTAGTAGAGGGAAAGGTCATTTTATAGGt                                 >  1:493048/1‑100 (MQ=255)
                                                                      cacCAAAGAGCGTTGAGATATAT‑ATTTATGA‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGtn                                <  2:2369452/100‑2 (MQ=255)
                                                                       acCAAAGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGtat                               >  2:2094327/1‑100 (MQ=255)
                                                                        ccAAAGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGtata                              >  2:3259161/1‑100 (MQ=255)
                                                                         cAAAGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATaa                             >  2:157925/1‑100 (MQ=255)
                                                                          aaaGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAAt                            >  1:2612850/1‑100 (MQ=255)
                                                                           aaGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATg                           >  2:1629442/1‑100 (MQ=255)
                                                                           aaGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATg                           <  1:1942762/100‑1 (MQ=255)
                                                                           aaGAGCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATATTACTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATg                           <  2:2241991/100‑1 (MQ=255)
                                                                           aaGAGCGTTGAGATATAT‑ATTTCTGA‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATg                           <  2:1087328/100‑1 (MQ=255)
                                                                            agagCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATgg                          <  2:2621595/100‑1 (MQ=255)
                                                                            agagCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATgg                          <  1:1341248/100‑1 (MQ=255)
                                                                               gcggtaagaTATAA‑CTTTCTGA‑AAAAAAATATTATGTATAAAATAATTAATGACCAAAAAATGAGTAGAGGGAATGGTCATTTTATAAGTATAATGGtt                        <  2:3910029/93‑1 (MQ=255)
                                                                               gCGTTGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTta                       >  2:681091/1‑100 (MQ=255)
                                                                                  ttGAGATATAT‑CTTTCTTA‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGTGAATGGTCATTTTATAAGTATAATGGtt                        <  1:2223698/96‑1 (MQ=255)
                                                                                  ttGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaata                    >  2:3652003/1‑100 (MQ=255)
                                                                                  ttGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaata                    >  1:1361756/1‑100 (MQ=255)
                                                                                  ttGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaat                     <  2:438336/99‑1 (MQ=255)
                                                                                  ttGAGATATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaa                      <  2:1097509/98‑1 (MQ=255)
                                                                                  ttGAGATATAT‑CTTTCTAA‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaata                    <  2:271511/100‑1 (MQ=255)
                                                                                   tGAGATATAT‑CTTTCTG‑‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaataa                   <  1:1372446/99‑1 (MQ=255)
                                                                                   tGAGATATAT‑CTTTAAAA‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTtaataa                   <  2:1212829/100‑1 (MQ=255)
                                                                                     agaTATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGATCATTTTATAAGTATAATGGTTAATAAtt                 >  2:408746/1‑100 (MQ=255)
                                                                                      gaTATAT‑CTTTCTGA‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTa                >  2:2425778/1‑100 (MQ=255)
                                                                                        tatat‑CTTTCTGA‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTAc               <  2:217550/99‑1 (MQ=255)
                                                                                         atat‑CTTTCTGA‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGCAATGGTCATTTTATAAGTATAATGGTTAATAATTACAg             <  2:2557380/100‑1 (MQ=255)
                                                                                            t‑CTTTCTGA‑AAAAAAATATTATTTATAAAATAATTACTGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATCATTACAGGAc          >  1:3946937/1‑100 (MQ=255)
                                                                                               tttCTGA‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTATTAGAGTGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGAc          >  2:3931452/1‑98 (MQ=255)
                                                                                                    g‑‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACtgtg  >  2:2872416/1‑100 (MQ=255)
                                                                                                     |                                                                                                    
TGTTTATTTGTTGTTTACCTTGGGTTTGATAACAATTCAAATTCTGACATCGTATTTTATGGACATAAAACACCAAAGAGCGTTGAGATATAT‑CTTTCTG‑‑AAAAAAATATTATTTATAAAATAATTAATGACCAAAAAATTAGTAGAGGGAATGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTG  >  NC_000913/1529970‑1530168

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑