Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,185,1990TG30.9% 126.0 / 12.6 69T35P (ACC→CCC) thiMhydroxyethylthiazole kinase
Reads supporting (aligned to +/- strand):  ref base T (8/39);  new base G (20/1);  total (29/40)
Fisher's exact test for biased strand distribution p-value = 6.69e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTTCATGACAAAAATGGCGGCGATAATCGAGCGCACCCACCGCTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATG  >  NC_000913/2184956‑2185440
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       
gTTCATGACAAAAATGGCGGCGATAATCGAGCGCACCCACCGCTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGcacc                                                                                                                                                                                                                                                 <  1:116312/251‑1 (MQ=255)
  gcATGACAAAAATGGCGGCGTTAATCGAGCGCACCCGCCGCTACGGGATCAAGCGTCCAGGGTGTTGGAGAGCTTTTTGCTTGCTCAACGGCATCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGcacca                                                                                                                                                                                                                                                <  2:125217/249‑1 (MQ=255)
   cATGACAAAAAGGGCGGCGATAATCGAGCGCACCCACCGCTACTGGATCAAGCGTCCAGGGTGGTTGAGAGCTGTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGcaccac                                                                                                                                                                                                                                               <  2:387436/250‑1 (MQ=255)
        cAAAAATGGCGGCGTTAATCGGGCGCACCCACCGCTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCAtt                                                                                                                                                                                                                                        <  2:164813/251‑1 (MQ=255)
         aaaaaTGGCGGCGATAATCGAGCGCACCCACCGCTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTTACCCACATCAttt                                                                                                                                                                                                                                       >  1:411806/1‑250 (MQ=255)
          aaaaTGGCGGCGATAATCGAGCGCACCCACCGCTACTGGATAAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAAGGTTGGACCAcatacttgg                                                                                                                                                                                                                                       >  2:312438/1‑245 (MQ=255)
           aaaTGGCGGCGATAATCGAGCGCACCCACCGCTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCGTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGt                                                                                                                                                                                                                                     <  2:173395/251‑1 (MQ=255)
            aaTGGCGGGGATAATCGAGCCCCCCCCCCGCTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTc                                                                                                                                                                                                                                    <  2:414104/251‑1 (MQ=255)
                gcggcgATAATCGAGCGCACCCACCGCTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGc                                                                                                                                                                                                                               <  1:210100/251‑1 (MQ=255)
                   gcgATAATCGAGCGCACCCCCCGCTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGt                                                                                                                                                                                                                            <  2:444680/251‑1 (MQ=255)
                     gATAATCGAGCGCACCCACCGCTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGc                                                                                                                                                                                                                          <  1:397293/251‑1 (MQ=255)
                         aTCGAGCGCACCCACCGCTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGGCCCACATCCATTggcgtgccgggccca                                                                                                                                                                                                                         >  1:192588/1‑237 (MQ=255)
                               cgcACCCACCGCTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAgggg                                                                                                                                                                                                               <  1:393754/251‑1 (MQ=255)
                                     caccGCTACTGGATCAAGGTTCCAGGGTGTTTGAGAGCTTTTGGCTTGCGCAACGGCAGCCCGCATCGTCTGAGCGCGTGGCGGCCTCAGTTTGCCAACGTTAATCAACAAGGAACTGGCGATAGCCGCAAACTGACGGGCCTCTTCGTTTTCGATAACCAGCGCTGGCGATGCACCGAGCGCCAGCAGGGTATTGGCGGGAAAGATTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGtt                                                                                                                                                                                                         <  2:386460/251‑1 (MQ=255)
                                     caccGCTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGCTAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGCAAAGGGTTGGCCCACATCNCTGGTccatacgtgcgacacaggggggtggt                                                                                                                                                                                                           >  2:219488/1‑226 (MQ=255)
                                         gCTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCCTTGGGCATGCAGTGCGCAAGAGGGGGATGTTgggg                                                                                                                                                                                                     >  1:60727/1‑249 (MQ=255)
                                            aCTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGaaa                                                                                                                                                                                                 <  1:326288/251‑1 (MQ=255)
                                                  tCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTGAAGGGTTGGACCACAGCCTTGGGCAGTGAGTGCACAAGAGGGGGAGGTTGGGGGAAAAgggg                                                                                                                                                                                           >  1:344732/1‑249 (MQ=255)
                                                   cAAGCGTCCAGGGTGGTTGAAAGCTTTGTGCTTGCTCAACGGCAGCCCGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGAGAGCCGCAAACTGACTGGCCTCTGCGGTTTCGATAACCAGCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGgtgt                                                                                                                                                                                          <  1:138853/251‑1 (MQ=255)
                                                                     gagaGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAAGTTTGGCCCCCATCCTTGGGCATGCAGTGGACAAGGAGGGGATGTTGGGGGAAAAGgtggctcggtttgggtttcgtt                                                                                                                                                                        >  1:89786/1‑232 (MQ=255)
                                                                                gcttgctCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCca                                                                                                                                                             <  2:84163/250‑1 (MQ=255)
                                                                                   tgctCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCcagc                                                                                                                                                           <  1:79463/249‑1 (MQ=255)
                                                                                      ccAACGGCCGCACGCATCGCCTGAGCGGGTGGGTGCGGCAGTGTTCCCACGTGAAGCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTTGCGGGAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAATAGGTGAATGTTGGTGAAAAAGGTGTAACGCGTACGCAGATTGCTCTGAACCCAGCAGGTc                                                                                                                                                      <  2:20104/250‑1 (MQ=255)
                                                                                       cAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTc                                                                                                                                                      <  1:473489/250‑1 (MQ=255)
                                                                                        aaCGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTTCACCCCAACATTGGTCATGCAGTGCCCAAGAGGGGGATGGTGGTGAAAAAGGTGTAACCCGCGCGCAGATTTCGCCGAAACCCGCAgggcgg                                                                                                                                                    >  1:69813/1‑247 (MQ=255)
                                                                                              agcaCGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCATTGGTCATGCAGGGCACAAGAGGGGGATGTTGGGGAAAAAGGGGTAACGCGTGCGCAGATTGGGCTGAACCCAGCAGGGCGGCTTGc                                                                                                                                               >  1:192385/1‑250 (MQ=255)
                                                                                                aaCGCATCGCCTGGGCGCGGGGCGGCGTCAGTGTGCAAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTTTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATa                                                                                                                                            <  2:56738/250‑1 (MQ=255)
                                                                                                 aCGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGGAAAGGTGTGGCACCCCTCCATTGTCATTCCGGGGACAAGAGGGGAATTGTGGGGGAAAAGGTGGAACGCGTGCGCAGATTGTGCTGTACCCAGAGGGTCGagtttcaaat                                                                                                                                           >  2:374608/1‑242 (MQ=255)
                                                                                                   gCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGGAAAGGGTTGCACCACACTATTTGTCATTCCCGGCGCGAGAGGGGGAAGTTGGGGGAACAAGGGGGAAGCGTGCGCGAATTGCGCTGAACCCAGCAGGGCGGCTTGCATAGtt                                                                                                                                         >  1:13900/1‑251 (MQ=255)
                                                                                                            gAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCt                                                                                                                                  <  2:434895/249‑1 (MQ=255)
                                                                                                                          gTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTgaagaa                                                                                                                  <  1:422084/251‑1 (MQ=255)
                                                                                                                          gTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGACAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTTCACCACATCATTGGTCATGCAGGGCACAAGAGGGGGATGGTGGGGGAAAAAGTGGAAAGCGTGCGCCGGTTGCGCCGAACCCCGCCGGGCGGCTTGCATTGTTTTCTCCTTCCATaaacggaaga                                                                                                                   >  1:487274/1‑242 (MQ=255)
                                                                                                                            cAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCATTGGTCATGCCGGGCACAAGAGGGGAATTTTGGGGAAAAAAGTGGAAAGCGTGCGCAGATTGCGCTGAACCCAGCCGGTCGGCTTGCCAAGGTTGCTCCTTCCCTAACGGGAAAAAGc                                                                                                                >  1:68963/1‑251 (MQ=255)
                                                                                                                            cAGTGTGCCAACGTAAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGc                                                                                                                <  2:363300/251‑1 (MQ=255)
                                                                                                                             aGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACCGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCATTGGGCACGCCGAGCACCACGGGGGGATGTTGGGGGAAAAGGTGTTAAGCGGGCGCAGATTGCGCTTGAACCACCAGGGCGGCTTGCATAGGTTGCTCCTTCCACAACggggcgaagg                                                                                                                >  1:278271/1‑241 (MQ=255)
                                                                                                                                 tgCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCCTTGGTGATGCCGTGCACAAGAGGGGAATGGTGGTGGAAAAGGTGTAACGCGGGCGCAGATTGCGCTGAACCCCGCAGGGCGACCTTCATAGTTTGCTCCTGCCATAAAGTGAAGAAGCAATg                                                                                                            >  1:295206/1‑250 (MQ=255)
                                                                                                                                          tAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGGAACACATCCTTGGTCATTCAGTGCACAAGAGGGGGAAGTTGGTGAAAAAAGTGTAACGCGTGCGCCGAATGGGCTGAACCCAGCAGGGCCGCTTGCATAGTTTGCTCCTGCCATAActggcagaaggcacgacccggtggg                                                                                                  >  1:94714/1‑227 (MQ=255)
                                                                                                                                           aaTCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACCTCATTGGTCATGCAGGGCCCAAGAGGGGGATGTTGGGGAAAAAGGGGTAACGCGGGCGCAGAtt                                                                                                                                                                          >  1:402379/1‑179 (MQ=255)
                                                                                                                                           aaTCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGAtt                                                                                                                                                                          <  2:402379/179‑1 (MQ=255)
                                                                                                                                               aacaaGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGGAAAGGGTTGCGCCACATCATTGGTCATGCAGTGCACAAGAGGGGAAAGTTGGTGGAAAAGGGGTAAAGCGTGCGCAGATTGCGCTGGAACCAGCAGGGCGACTTGCATAGTTTGCTCCCGCCCTAAAGTGGAGAAGCCATGACCTGGTGGGCCGTg                                                                                             >  1:29618/1‑251 (MQ=255)
                                                                                                                                                 caaGGCCCTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCGGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGAc                                                                                           <  2:192385/251‑1 (MQ=255)
                                                                                                                                                  agggCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACt                                                                                          <  1:196661/249‑1 (MQ=255)
                                                                                                                                                  aaGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCAACACAACATTGGTCACGCCGTGCACAAGAGGGGGGGGTTGGGGAAAAAGGGGTAACGCGTGCGGAAATTGCGGTGAACCCAGGAGGGGGACTTGGATAGTTTGCTCCTGCCAGAAAGGGAAGAAGCAATGAACTGGGGGGCCGTggtt                                                                                          >  1:467588/1‑248 (MQ=255)
                                                                                                                                                       aCTGGCGATAGCCGCAAACTGACCGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCt                                                                                     <  1:428170/251‑1 (MQ=255)
                                                                                                                                                          ggCGATAGCCGCAAACTGGCTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGTGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATGGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGACCCCAGCAGGTCGACTTGCATAGTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTAc                                                                                   <  1:29852/249‑1 (MQ=255)
                                                                                                                                                          ggCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACg                                                                                  <  2:452878/251‑1 (MQ=255)
                                                                                                                                                            cgagaGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGGCCGTGGCTTCCCTACGCt                                                                                <  2:415631/247‑1 (MQ=255)
                                                                                                                                                            cGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCt                                                                                <  2:194976/251‑1 (MQ=255)
                                                                                                                                                             gATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTgcgc                                                                                                                                                                      <  2:31476/165‑1 (MQ=255)
                                                                                                                                                             gATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAAGTGTAACGCGTGCGCAGATTgcgc                                                                                                                                                                      >  1:31476/1‑165 (MQ=255)
                                                                                                                                                                   cGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCATTGGGTATGCAGTGCACAAGAGGGGGATGTTGGGGAAAAAAGGGTAACGCGTGCGCAGATTTGGCCGAAACCAGCAGGGCGACTTGCATAGTTTGCTCCCTCCCTAACGTGAAGAAGCAATGACCCGGTGGGCCGTGACttaccccccccgggctt                                                                          >  2:319558/1‑235 (MQ=255)
                                                                                                                                                                        actgactgGCCGTTTCGGTGTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAg                                                                    <  2:192588/251‑1 (MQ=255)
                                                                                                                                                                         ctgactgGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTGTTGACCACCTCATTTGGCATGCAGGGCACAAGAGGGGGGTGTTGGGGAAAAAGGGGTAACGCGTGCGCAGATTTCGCTGAACCCACCAGGTCGGCTTGCCTAGTTTTCTCCTGCCATAAAGTGCAGAAGGAATGACCAGGTGGGCCGTGACTTCCCCACGCTGGCATTAttccg                                                                    >  2:445528/1‑246 (MQ=255)
                                                                                                                                                                              tgGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGGATGTTGGGGAAAAAGGTGTAACGCGGGCGCAGATTGCGCTGAACCCAGCAGGTCGGCTTTCATAGGTTGCTCCTGCCCTAACGTGGAGAAGGAATGACCTGGGGGGCCGTGACTTCCCCACGCCGGCATTATCCAGATCAgg                                                              >  1:130565/1‑251 (MQ=255)
                                                                                                                                                                               ggCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAg                                                               <  2:411806/249‑1 (MQ=255)
                                                                                                                                                                                      tCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATAc                                                        <  1:243876/249‑1 (MQ=255)
                                                                                                                                                                                      tCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACg                                                       <  1:480421/250‑1 (MQ=255)
                                                                                                                                                                                         gTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCAACCCATCATTGGTGATGCAGTGGACAAGGGGGGGATGTTGGTGAAAAAAGTGTAACGCGTGCGCAGATTGCGCTGAAACCAGCAGGGCGGCTTGCATTGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTTCGCTGGCATTATCCAGATCAGGAGGTACggggt                                                   >  1:108111/1‑249 (MQ=255)
                                                                                                                                                                                             cGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGAAAAGGGTTGGACCACATCCTTGGGCATGCAGGGGACAAGAGGGGGATGTTGGGGGAAAAGGTGTAACGCGTGCGCCGATTGCGCTGAACCCAGCAGGGCGACGTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGGCCGTGGCTTCCCTCCGACGNCATTATCCAGATGAGGTGGTACGGGTAtttt                                               >  2:232850/1‑250 (MQ=255)
                                                                                                                                                                                               aTAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTtctc                                             <  2:411825/251‑1 (MQ=255)
                                                                                                                                                                                                  aCCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAAGGTTTCCCCACATCATTGGTCATGCAGTGGACCAGAGGGGGATGTTGGTGAAAAAAGGGTAACGCGGGGGCCGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAAAGTGAAAAAGCAATGACCTGGGGGTCCGTGAATTTCCTACGCTGGCATTATCCAGATCAGGGGATACGGGGATTTCTCAGc                                          >  1:315234/1‑251 (MQ=255)
                                                                                                                                                                                                          tGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGGTTTCACCCCATCATTTGGCATGCAAGGAACACGAGGGGGATGGTGGGGGAAAAGGTGTAACGCGGGCGCAGAATTCGGTGAACCCCGCCGGGCGGATTGCCTTGTTTGCTCCCTCCATAACCGGAAGAAACAATGACCTGGTGGGCCGTGGCTTCCCTACCCCGGCCTTATCCAGATCAGGGGGTACGGGGTTTTTTCAGCCtttccgc                                  >  1:214766/1‑246 (MQ=255)
                                                                                                                                                                                                                       ggcgcCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCAcccc                     <  2:425394/250‑1 (MQ=255)
                                                                                                                                                                                                                        gcgcCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGt                                                                                                                                                                                                      <  2:49966/75‑1 (MQ=255)
                                                                                                                                                                                                                           ccAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGt                 <  1:146654/251‑1 (MQ=255)
                                                                                                                                                                                                                            cAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCa                                                                                                                                                                                                                              <  1:200954/48‑1 (MQ=255)
                                                                                                                                                                                                                               cAAGGTATTGGCGGTAAAGGGTTGCACCACATCATTGGTCATGCAGGGCACAAAGGGGGAATGTTGGGGGAAAAAGGGTAACGCGTGCGCAAAATGCGCTGAACCCAGCAGGTCGACCTGTATAGTTTTCTCCTGCCCTAACGTGAAAAAGCAATGACCTGGGGGTCCGGGACTTCCCTACGCTGGCATTATCCAGATCAGGGGATACGGGGATTTTTCAGCCTTTACGCAGAAGGGCACCCCGAGTCGtt             >  1:461296/1‑251 (MQ=255)
                                                                                                                                                                                                                                        ggCGGGAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCTAGTCGTTTGGTTGCg     <  2:315234/250‑1 (MQ=255)
                                                                                                                                                                                                                                           ggTAAAGGATTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGCCTTCCCTACGCTGGCATTATCCAGCTCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAACGGCACCCCGCGTCGCTTGGTTGCgatg  >  2:265432/1‑250 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       
GTTCATGACAAAAATGGCGGCGATAATCGAGCGCACCCACCGCTACTGGATCAAGCGTCCAGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCATAGTTTGCTCCTGCCATAACGTGAAGAAGCAATGACCTGGTGGTCCGTGACTTCCCTACGCTGGCATTATCCAGATCAGGTGATACGGGTATTTCTCAGCCTTCACGCAGAAGGGCACCCCGAGTCGTTTGGTTGCGATG  >  NC_000913/2184956‑2185440

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: