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% 724.7 / 53.3 199intergenic (+1/‑2)valT/lysWtRNA‑Val/tRNA‑Lys
Reads supporting (aligned to +/- strand):  ref base . (80/105);  new base T (8/4);  total (90/109)
Fisher's exact test for biased strand distribution p-value = 1.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 207 total aligned reads displayed.

GGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCG  >  NC_000913/780693‑780985
                                                                                                                                                     |                                                                                                                                                
ggCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCaccacc                                                                                                                                                 <  1:1881640/150‑1 (MQ=255)
 gCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCg                                                                                                                                                >  1:568269/1‑150 (MQ=255)
 gCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCGGGGAGAGCACCTCCCCTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCg                                                                                                                                                >  1:4042918/1‑150 (MQ=255)
 gCGAGACGAAGTCGAGACATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCTGTAGATAAGCTCAGCTGGGAGAGCACCTCACTTACAAGGAGGGGGTCGGCGGTTCGATCCCTTAATCACCcacaacag                                                                                                                                                 >  2:3944042/1‑148 (MQ=21)
  cGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCgg                                                                                                                                               <  1:437964/150‑1 (MQ=255)
   gagaCGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCggg                                                                                                                                              <  1:2221904/150‑1 (MQ=255)
    agaCGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGt                                                                                                                                             >  2:3082783/1‑150 (MQ=255)
        gAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGt                                                                                                                                          >  2:278539/1‑150 (MQ=255)
         aaGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTa                                                                                                                                        <  1:2068328/150‑1 (MQ=255)
              gAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTc                                                                                                                                    >  1:4029268/1‑150 (MQ=255)
              gAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCa                                                                                                                                   >  2:1107615/1‑150 (MQ=255)
              gAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCa                                                                                                                                   <  1:2038735/150‑1 (MQ=255)
              gAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCa                                                                                                                                   <  1:2371854/150‑1 (MQ=255)
              gAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCa                                                                                                                                   <  1:2269882/150‑1 (MQ=255)
               aGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACGGGGTCGTTAGCTCAg                                                                                                                                  >  1:356896/1‑150 (MQ=255)
               aGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGt                                                                                                                                          >  1:3288212/1‑142 (MQ=255)
               aGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAg                                                                                                                                  >  1:356807/1‑150 (MQ=255)
               aGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAg                                                                                                                                  <  2:3013610/150‑1 (MQ=255)
                gTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAg                                                                                                                                  >  2:2356397/1‑150 (MQ=255)
                gTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGt                                                                                                                                 >  2:2712319/1‑150 (MQ=255)
                 tCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGtt                                                                                                                                <  1:3305496/150‑1 (MQ=255)
                 tCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCGCAGtt                                                                                                                                >  1:3509250/1‑150 (MQ=255)
                  cATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTCGCACAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCGCAGTTg                                                                                                                               >  1:3783055/1‑150 (MQ=255)
                   aTCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTgg                                                                                                                              >  2:2647575/1‑150 (MQ=255)
                   aTCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTgg                                                                                                                              <  1:1573263/150‑1 (MQ=255)
                   aTCCTGCACGACCCACCACTAACAGAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTgg                                                                                                                              <  1:128598/150‑1 (MQ=255)
                    tCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGt                                                                                                                             <  1:130864/150‑1 (MQ=255)
                     ccTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTa                                                                                                                            <  1:1971913/150‑1 (MQ=255)
                        gCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTagag                                                                                                                         <  2:1146281/150‑1 (MQ=255)
                        gCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTagag                                                                                                                         <  1:62173/150‑1 (MQ=255)
                        gCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTagag                                                                                                                         >  1:2924516/1‑150 (MQ=255)
                          aCGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCa                                                                                                                       <  2:3667521/150‑1 (MQ=255)
                            gACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGt                                                                                                                     <  2:2798979/150‑1 (MQ=255)
                             aCCCACCACTACCATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGtt                                                                                                                    <  2:2834253/150‑1 (MQ=255)
                             aCCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGtt                                                                                                                    <  2:2625138/150‑1 (MQ=255)
                               ccaccaCTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGa                                                                                                                  >  1:2721373/1‑150 (MQ=255)
                               ccaccaCTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGa                                                                                                                  <  2:3043533/150‑1 (MQ=255)
                               ccaccaCTAACAGAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACCAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGa                                                                                                                  <  2:2125188/150‑1 (MQ=255)
                                   cacTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACttt                                                                                                              <  1:666403/150‑1 (MQ=255)
                                    acTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGGAGAGCAGTTGACttt                                                                                                              >  2:497040/1‑150 (MQ=255)
                                    acTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACtttt                                                                                                             >  2:3303000/1‑150 (MQ=255)
                                    acTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACtttt                                                                                                             <  1:2633813/150‑1 (MQ=255)
                                    acTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGGCGTTAGGTCAGTTGGTAGTGCAGTTTACtttt                                                                                                             >  2:2544553/1‑150 (MQ=255)
                                     cTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTa                                                                                                            <  1:2499627/150‑1 (MQ=255)
                                     cTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTa                                                                                                            <  2:2028603/150‑1 (MQ=255)
                                      taacccagttagttGTAGTACCCAGCCTATTATCGGGTGCTTATCTCAGAATTGAGAGAAACTCCCTTACCTGGACGGGGTATCCGGTTCGATCCGGTCAACACCCACAACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTaa                                                                                                           <  2:709226/144‑1 (MQ=9)
                                      tAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTaa                                                                                                           <  1:3499925/150‑1 (MQ=255)
                                          ctagttagtTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACTCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGAc                                                                                                                 <  1:1080422/140‑1 (MQ=17)
                                           tagttagtTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAAt                                                                                                      >  1:1707438/1‑150 (MQ=255)
                                            agttagttGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAAtt                                                                                                     <  1:2806127/150‑1 (MQ=40)
                                               tagttGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGt                                                                                                  <  1:1843856/150‑1 (MQ=38)
                                                 gttGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGTGAGAGCACCTCCCTTACATGGAGGGGGTCGGCGGTACGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCg                                                                                                <  2:1356257/150‑1 (MQ=255)
                                                  ttGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACACGGGACGTAAGCTCAGTTGGTAGACAGGTTGAATTTTAAGAATTGGtgcg                                                                                               >  2:1117918/1‑147 (MQ=1)
                                                       gtaTCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGtt                                                                                          <  1:1388867/150‑1 (MQ=35)
                                                        taTCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTc                                                                                         <  1:2986852/150‑1 (MQ=35)
                                                        taTCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTc                                                                                         >  1:2551546/1‑150 (MQ=35)
                                                        taTCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTc                                                                                         <  2:2151465/150‑1 (MQ=35)
                                                          tCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGa                                                                                       <  2:1343313/150‑1 (MQ=35)
                                                            cAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAAt                                                                                     >  2:1689595/1‑150 (MQ=35)
                                                             aGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATc                                                                                    >  1:2179991/1‑150 (MQ=35)
                                                              gCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATcc                                                                                   >  1:2183482/1‑150 (MQ=35)
                                                               attattaTCGGGAGATTATCTCTTCCGGGAGAGCACCTCCCTTCGAAGTAGTGGGTATGCGTTTCTATCCCCTCTTCACCCAACACAGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCt                                                                                  <  2:3172943/145‑1 (MQ=11)
                                                                       cGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGAcc                                                                          <  2:818114/150‑1 (MQ=33)
                                                                        gggTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCGCAGTTGGTAGAGCCGTTGGCTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGAccc                                                                         >  1:2109466/1‑150 (MQ=16)
                                                                          gTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACccac                                                                       <  2:568269/150‑1 (MQ=33)
                                                                            gATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACccacca                                                                     <  1:4045202/150‑1 (MQ=33)
                                                                               tAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGtt                                                                  >  2:2993198/1‑150 (MQ=33)
                                                                               tAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGCTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGtt                                                                  <  2:3118006/150‑1 (MQ=18)
                                                                                     aGCTGGGAGAGCACCTCCCTTACAAGGTGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACa                                                            <  1:1220676/150‑1 (MQ=18)
                                                                                       cTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATc                                                          <  1:3396729/150‑1 (MQ=33)
                                                                                          ggAGAGCACCTCCCTTACAAGGTGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGaa                                                       <  1:2371565/150‑1 (MQ=25)
                                                                                          ggAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCaaa                                                       <  2:2183482/150‑1 (MQ=14)
                                                                                            agagCACCTCCCTTACTAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGa                                                     <  2:437163/150‑1 (MQ=38)
                                                                                              agCTCCTCCCTTACTAGGAGTGGGTCGGCGTTTCGATCCCGTCTTCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACa                                                   <  2:2576807/150‑1 (MQ=14)
                                                                                              agCACCTCCCTTACAAGGAGGGGGTCGGCTGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACa                                                   <  1:762249/150‑1 (MQ=255)
                                                                                              agCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACa                                                   <  2:3975475/150‑1 (MQ=39)
                                                                                                  ccTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATg                                               >  2:1021882/1‑150 (MQ=40)
                                                                                                      ccTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTaa                                           >  1:676855/1‑150 (MQ=40)
                                                                                                       cTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAg                                          <  1:3374045/150‑1 (MQ=40)
                                                                                                        ttACACGGAGGGGGTCGGCGTTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGTGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGc                                         <  2:3824087/150‑1 (MQ=255)
                                                                                                             aggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTa                                    <  1:1994703/150‑1 (MQ=255)
                                                                                                              ggaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTCTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAg                                   <  2:1838328/150‑1 (MQ=255)
                                                                                                               gaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAgg                                  >  2:3879717/1‑150 (MQ=255)
                                                                                                               gaggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCAGCACGACCCACCCGTTTTAACATCCAACACAGATGTTAAGCGTGAAgg                                  >  2:3537942/1‑150 (MQ=11)
                                                                                                                aggGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGa                                 <  1:3082783/150‑1 (MQ=255)
                                                                                                                  ggggTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAATCGTGTATGATa                               >  2:2412455/1‑150 (MQ=255)
                                                                                                                       cggcggTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGtt                          <  2:2924516/150‑1 (MQ=255)
                                                                                                                        ggcggTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTg                         >  2:3450853/1‑150 (MQ=255)
                                                                                                                         gcggTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCAGCACGACCCACCCGTTTTAACATCGAAGACAGATGTTTAGCGTGTAGGATAACGTTGc                        >  1:1813468/1‑150 (MQ=255)
                                                                                                                             ttCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCa                    <  2:3746293/150‑1 (MQ=255)
                                                                                                                             ttCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCa                    <  2:4036605/150‑1 (MQ=255)
                                                                                                                              tCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAg                   <  2:2721373/150‑1 (MQ=255)
                                                                                                                               cGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTTGTTGAGCAGTTGACTTTTAATAAATTTGTCGCAGGTTGGAATACTGCACCACCCCCCAGTTTTTAAATCGAAGAAACATGTTAAGAGGGTAGGATAAAGGTGCgcccgc                  >  1:3797307/1‑145 (MQ=14)
                                                                                                                                 aTCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGAGGTTAAGCGTGTAGGATAACGTTGCGTCAGCaa                >  2:2605593/1‑150 (MQ=255)
                                                                                                                                 aTCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGGTGGTAGAGCAGTTGACTTTTAATCAATTTGTCGCAGGGTCGAATCCTGCACGACCCACCACTTTTAAAATCGAATACAGGTGTTAAGCGTGTAGGATAACGTTGCGTCAGCaa                >  1:3756456/1‑150 (MQ=14)
                                                                                                                                    ccGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGGTAAGCGTGTAGGATAACGTTGCGTCAGCAACgg             >  2:2834907/1‑150 (MQ=255)
                                                                                                                                    ccGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGGTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACgg             >  2:767987/1‑150 (MQ=255)
                                                                                                                                         atcaCCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGt        >  2:3273960/1‑150 (MQ=255)
                                                                                                                                           cacccaccACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAg      >  2:161792/1‑150 (MQ=255)
                                                                                                                                           aacgaaccagcGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAg      <  2:2315338/145‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
GGCGAGACGAAGTCGAGTCATCCTGCACGACCCACCACTAACATAGTTAGTTGTAGTATCCAGCGTAGTATCGGGTGATTAGCTCAGCTGGGAGAGCACCTCCCTTACAAGGAGGGGGTCGGCGGTTCGATCCCGTCATCACCCACCACCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAGTTTTAACATCGAAGACAGATGTTAAGCGTGTAGGATAACGTTGCGTCAGCAACGGCCCGTAGGGCG  >  NC_000913/780693‑780985

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: