Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 780,841:1 +T 7.0% intergenic (+1/‑2) valT → / → lysW tRNA‑Val/tRNA‑Lys

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913780,8411.T 7.0% 752.7 / 51.5 198intergenic (+1/‑2)valT/lysWtRNA‑Val/tRNA‑Lys
Reads supporting (aligned to +/- strand):  ref base . (66/118);  new base T (6/6);  total (74/124)
Fisher's exact test for biased strand distribution p-value = 3.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.63e-01

Only 100 of 199 total aligned reads displayed.

CGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGC  >  NC_000913/780695‑780984
                                                                                                                                                   |                                                                                                                                               
cGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAATGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCgg                                                                                                                                              >  2:225674/1‑150 (MQ=255)
cGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCgg                                                                                                                                              <  1:847759/150‑1 (MQ=255)
  agaCGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCcaccactcggg                                                                                                                                             >  1:1017580/1‑145 (MQ=255)
   gaCGAAGTCGAGTCATCCTGCACGACCCACCACTAACATCGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACCAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTc                                                                                                                                           <  2:1595083/150‑1 (MQ=255)
     cGAAGTCGAGTCATCCTGCACGACCCACCACTACCATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGt                                                                                                                                         <  2:1795249/150‑1 (MQ=255)
     cGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGt                                                                                                                                         >  1:2042024/1‑150 (MQ=255)
       aaGTCGAGTCATCCTGCACGACCCACCACTCCCATAGTTAGTTGTAGTATCCAGCGTTGTATCGGGTGATTAGCTCAGCTGGGAGAGCTCCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATGCCGTCATCACCCACCACCGGGTCGTTa                                                                                                                                       <  2:1697752/150‑1 (MQ=18)
       aaGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTa                                                                                                                                       <  1:3708753/150‑1 (MQ=255)
        aGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAg                                                                                                                                      >  1:3507617/1‑150 (MQ=255)
        aGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAg                                                                                                                                      >  1:3325738/1‑150 (MQ=255)
        aGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAg                                                                                                                                      <  1:2700569/150‑1 (MQ=255)
        aGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAg                                                                                                                                      >  2:3836917/1‑150 (MQ=255)
         gTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGTGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGc                                                                                                                                     >  1:3652111/1‑150 (MQ=255)
            gAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTcc                                                                                                                                  >  2:1622507/1‑149 (MQ=255)
            gAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCa                                                                                                                                  >  2:4011196/1‑150 (MQ=255)
             aGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAg                                                                                                                                 >  2:1663353/1‑150 (MQ=255)
               tCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGtt                                                                                                                               <  1:2296364/150‑1 (MQ=255)
                 aTCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTgg                                                                                                                             <  1:2903375/150‑1 (MQ=255)
                   ccTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGt                                                                                                                            >  2:3029957/1‑150 (MQ=255)
                     tGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTaga                                                                                                                         <  2:3249877/150‑1 (MQ=255)
                     tGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTaga                                                                                                                         <  2:2617444/150‑1 (MQ=255)
                     tGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTaga                                                                                                                         <  1:11785/150‑1 (MQ=255)
                      gCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTagag                                                                                                                        <  2:3495158/150‑1 (MQ=255)
                       cACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGc                                                                                                                       <  1:2934789/150‑1 (MQ=255)
                        aCGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGTGTGTTTAGCTCAGCTGTGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCa                                                                                                                      <  2:1868767/150‑1 (MQ=255)
                          gACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAg                                                                                                                     <  1:1312525/150‑1 (MQ=255)
                          gACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGt                                                                                                                    <  1:884605/150‑1 (MQ=255)
                          gACCCACCACTAACAGAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGt                                                                                                                    <  2:722120/150‑1 (MQ=255)
                             ccaccaCTAACCTAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGt                                                                                                                                         >  1:973980/1‑127 (MQ=255)
                             ccaccaCTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGt                                                                                                                                         <  2:973980/127‑1 (MQ=255)
                             ccaccaCTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGa                                                                                                                 >  2:2928837/1‑150 (MQ=255)
                              caccacTACCATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGAc                                                                                                                <  2:173166/150‑1 (MQ=255)
                               accacTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACt                                                                                                               >  2:2670428/1‑150 (MQ=255)
                                ccacTAACCTAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCGGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACtt                                                                                                              <  2:3634696/150‑1 (MQ=255)
                                ccacTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACACGGGACTTAAGAAAAGTTGGAAGACCAGtttata                                                                                                               >  2:2635821/1‑146 (MQ=14)
                                ccacTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACtt                                                                                                              <  1:2452063/150‑1 (MQ=255)
                                 cacTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACttt                                                                                                             <  1:3943734/150‑1 (MQ=255)
                                 cacTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACttt                                                                                                             <  1:2858067/150‑1 (MQ=255)
                                    tAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTaa                                                                                                          <  2:804874/150‑1 (MQ=255)
                                    tAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATGAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTaa                                                                                                          <  2:1353033/150‑1 (MQ=255)
                                     aaCATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAAt                                                                                                         >  1:4163160/1‑150 (MQ=255)
                                     aaCATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAAt                                                                                                         >  1:2361263/1‑150 (MQ=255)
                                        aTAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCaa                                                                                                      <  1:390751/150‑1 (MQ=255)
                                        aTAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCaa                                                                                                      <  1:3188991/150‑1 (MQ=255)
                                         tagtttgtTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCaa                                                                                                      <  1:3029957/150‑1 (MQ=2)
                                         tagttagtTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCaa                                                                                                      <  2:2479475/150‑1 (MQ=17)
                                         tagttagtTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAAt                                                                                                     <  2:1462292/150‑1 (MQ=255)
                                         tagttagtTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAAt                                                                                                     <  2:1462296/150‑1 (MQ=255)
                                         tagttagtTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAAt                                                                                                     <  2:747306/150‑1 (MQ=255)
                                          agttagttGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAAtt                                                                                                    <  2:489841/150‑1 (MQ=40)
                                              agttGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCTCCTCCCTTACCAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTc                                                                                                <  2:189509/150‑1 (MQ=255)
                                                ttGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGc                                                                                              >  1:1053868/1‑150 (MQ=35)
                                                  gtagtaTCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAg                                                                                            <  2:4126534/150‑1 (MQ=35)
                                                  gtagtaTCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAg                                                                                            <  2:308104/150‑1 (MQ=35)
                                                      taTCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTc                                                                                        >  1:3799001/1‑150 (MQ=35)
                                                       aTCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACGCGGGGCGTACGCGCAGTGGGTCGGGCAGTTGAATTTTAAGCAATTGGTCGCAGGtac                                                                                        >  2:2062218/1‑148 (MQ=1)
                                                        tCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGa                                                                                      >  1:1229307/1‑150 (MQ=35)
                                                        tCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGGTGGTAGAGCAGTTTACATTTAATCAAATGGTAGCAGGTTAGa                                                                                      >  2:1519424/1‑150 (MQ=11)
                                                             cGTAGTATCGGGTGATTAGCTCAGCTGGGAGGGCCCCAACCTTGACAGGAGGGGGTCGGCGGTTCGATACCGTCAACAACAACCACCGCGGCGTTAGATTAGTTGGTAGTGCGGGTGACTTTTGGTGAGGTGGTCGCGGGTTCGATTcgt                                                                                 >  2:653092/1‑148 (MQ=2)
                                                               tagtaTCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGc                                                                               >  1:4093252/1‑150 (MQ=34)
                                                               tagtaTCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGc                                                                               <  2:2743812/150‑1 (MQ=34)
                                                               tagtaTCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGACATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGc                                                                               >  1:1985282/1‑150 (MQ=21)
                                                                agtaTCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACCAGGAGTGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCa                                                                              <  2:247696/150‑1 (MQ=18)
                                                                    tCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGAc                                                                          >  1:1803009/1‑150 (MQ=33)
                                                                             tAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGtt                                                                 <  2:2042024/150‑1 (MQ=33)
                                                                              aGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGttt                                                                >  1:2304800/1‑150 (MQ=33)
                                                                              aGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGttt                                                                <  2:1128621/150‑1 (MQ=33)
                                                                                cTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTa                                                              >  1:446714/1‑150 (MQ=33)
                                                                                 tCAGCTTGGAGATCACCTCCCTTACAAGGAGGGGGTCGTAGGTTCTATCCCGTCATCACCAACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTaa                                                             <  1:4076961/150‑1 (MQ=11)
                                                                                 tCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCgg                                                                                                                                              <  2:2871221/69‑1 (MQ=35)
                                                                                 tCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTaa                                                             <  1:507470/150‑1 (MQ=33)
                                                                                    gCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACAt                                                          <  2:2245681/150‑1 (MQ=33)
                                                                                    gCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACAt                                                          <  2:898047/150‑1 (MQ=33)
                                                                                        ggAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGaa                                                      <  1:2678235/150‑1 (MQ=35)
                                                                                          tgagCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATAACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCTAATCCTGCACGACCCACCAGTTTTAACATCGAAGa                                                    <  1:2008783/149‑1 (MQ=255)
                                                                                              cACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCACAGGTGGTAGAGCAGTTGACTTTTAATCAAATGGTCGCACGTTAGAATCCTGCACCACCCACCAGTTTTAACAAAGAACACAGa                                                >  2:360474/1‑150 (MQ=11)
                                                                                                    ccTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGTCTTTTAAACAAATTTTAGCAGGTTCGAATACGGCACCACCCAACCCTTTTAAAAACAAAGACACATGTTaa                                          >  1:1932059/1‑150 (MQ=11)
                                                                                                          aaGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGTGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCgtgt                                    >  2:1999922/1‑150 (MQ=255)
                                                                                                           aggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTCGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTCTTAACATCGAAGACAGATGTTAAGCGTGTa                                   <  1:1355373/150‑1 (MQ=255)
                                                                                                            ggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAg                                  <  1:3836917/150‑1 (MQ=255)
                                                                                                            ggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAg                                  <  2:3326173/150‑1 (MQ=255)
                                                                                                            ggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAg                                  <  1:3327866/150‑1 (MQ=255)
                                                                                                            ggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAg                                  <  1:3723904/150‑1 (MQ=255)
                                                                                                             gaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAgg                                 <  1:1663353/150‑1 (MQ=255)
                                                                                                                ggggTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGGGCCGGTGGTAGAGCAGTTGACATTTAATAAATTGGTCGCAGGTTCGAATCCAGCACGACCCACCAGTTTTAAAAAAAAAAAAACGTGTTAAGCGTGTAGGATa                              >  1:917084/1‑150 (MQ=11)
                                                                                                                 gggTCGGCGGTTCGATCCCGTCATCACCCACCAGCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATaa                             <  2:3943671/150‑1 (MQ=255)
                                                                                                                 gggTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATaa                             <  1:3326404/150‑1 (MQ=255)
                                                                                                                    tCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGt                          >  2:2977634/1‑150 (MQ=255)
                                                                                                                       gcggTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGc                       >  2:1575103/1‑150 (MQ=255)
                                                                                                                        cggTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCg                      <  1:4091751/150‑1 (MQ=255)
                                                                                                                         ggTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGGTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGGTCGAATCCGGCACGACCCACCAGTTTTAACAAAGAAGACAGAGGGTAAGAGTGTAGGATAAAGGTGCGt                     >  1:2573477/1‑150 (MQ=14)
                                                                                                                           ttCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGTGGAGTTGACTTTTCATCTATTGGCCTCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCa                   <  1:650616/150‑1 (MQ=14)
                                                                                                                               aTCCCGTCATCACCCACCACCGGGTCGTTAGGTCAGTTGGGAGAGCAGGTGTCTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCACATTTAAAAACGAAGAAAGAAGTTAAGCGTGTAGGATAAAGTTGCGTCAGCaa               >  1:3351063/1‑150 (MQ=14)
                                                                                                                               aTCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCaa               >  2:2740632/1‑150 (MQ=255)
                                                                                                                                tCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAAc              <  1:2977634/150‑1 (MQ=255)
                                                                                                                                   cggcatcaCCCACCACCGGGGCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGc           <  2:1361637/147‑1 (MQ=255)
                                                                                                                                    gTCATCACCAACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGcc          <  2:2153561/150‑1 (MQ=255)
                                                                                                                                        tcaCCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTa      <  2:4080836/150‑1 (MQ=255)
                                                                                                                                        tcaCCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTa      <  2:2304800/150‑1 (MQ=255)
                                                                                                                                            ccaccaCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGgc  >  2:3242012/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
CGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGC  >  NC_000913/780695‑780984

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