Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913836,7560TG30.8% 143.2 / 11.4 65F136V (TTT→GTT) hcxBhydroxycarboxylate dehydrogenase B
Reads supporting (aligned to +/- strand):  ref base T (7/38);  new base G (20/0);  total (27/38)
Fisher's exact test for biased strand distribution p-value = 3.13e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTACGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGTCCCCGTGCCGCCAGGTTGCCTGATTGACGTTAACGGCGTGCCGACGACCAATCCG  >  NC_000913/836517‑837001
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                        
aTGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTAt                                                                                                                                                                                                                                                <  1:497636/250‑1 (MQ=255)
               gTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTg                                                                                                                                                                                                                                <  1:569811/251‑1 (MQ=255)
                        cccTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTAtt                                                                                                                                                                                                                       <  1:315872/250‑1 (MQ=255)
                                   aaCCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTAGCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTcg                                                                                                                                                                                                             <  1:481537/250‑1 (MQ=255)
                                        gcatGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGt                                                                                                                                                                                                     <  2:602635/250‑1 (MQ=255)
                                                    cGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTcgcg                                                                                                                                                                                         <  1:329390/251‑1 (MQ=255)
                                                         aaGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGc                                                                                                                                                                                    <  1:698300/251‑1 (MQ=255)
                                                          agagGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGcc                                                                                                                                                                                   <  1:359567/251‑1 (MQ=255)
                                                            agGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGGTTGTCTCTATCCACTTTGTTAGCGGGGTCGGGATTCCAATGGTCCCGCCGGTCCCCGGTCGCGGCAACCGc                                                                                                                                                                                 >  1:593650/1‑251 (MQ=255)
                                                                 ggggCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGGTTGTCTCTATCCACTTTGTTTGCGTGGTCGGGATTTCAATTGTCGCGCCGCTCCACGGGCGCGACAGCCCCTTTgg                                                                                                                                                                           >  2:733197/1‑251 (MQ=255)
                                                                  gggCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGGTTGTCTCTTTCCACTTTGTTAGCGTGGTCGGGGTTCCAACGGTCGCGCCGTTCCCCGGGCCCGACAGCCGCGTTTgg                                                                                                                                                                           >  2:96052/1‑251 (MQ=255)
                                                                    gcggcgGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCAc                                                                                                                                                                        <  1:358506/251‑1 (MQ=255)
                                                                     gggcggtcccgtgcGATGGCGGTCGCGCATTTGATATGGTGGCGGCACATGAAGCGAGGGCGGGGGGGATGGAGAAGGCGCATCAGCACGGTTGTGCCGCCGTGGCTCTACATAACTCGCTTCATATCTGCCGTATCGGTTACTGGGCGGAGCAGTGTGCGGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCTCGCCGTTCCACGTTCGCGACAGCCGCTTTGGCAcc                                                                                                                                                                       <  2:348641/238‑1 (MQ=255)
                                                                     gggcggTCACGCTCGATGGCGTGCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGAGGGCGCTGGGGATTGAGAAAGCGAATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCAcc                                                                                                                                                                       <  2:741514/250‑1 (MQ=255)
                                                                     cggcggTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGGTTTTCTCTATCCCCTTTGTTAGCGTGGGCGGGGTTTCAATGGGCGCGCCGTTCCCCGGGCGCGACAGCCGCttttggccc                                                                                                                                                                       >  1:608490/1‑245 (MQ=255)
                                                                          gTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGGTTGTCTTTTTTCACTTTGTTTGAGGGGTCGGGCTTCCAATGGTCGCGCCGTTCCCCGGTCGCGACCGCCGCTTTgggcccactcc                                                                                                                                                                  >  2:40568/1‑242 (MQ=255)
                                                                          gTCACCCTCGATGGCGATCGGGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATACCTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATcc                                                                                                                                                                  <  2:101926/251‑1 (MQ=255)
                                                                            cACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGt                                                                                                                                                                <  1:28160/251‑1 (MQ=255)
                                                                              cGCTCGATGGCGGTCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGGTTGTCTCTATCCCCTTTGTTAGCGTGGGCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCCCCAATCCGTTc                                                                                                                                                              >  1:522546/1‑251 (MQ=255)
                                                                                   gATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGGTTGTCTCTATACACTTTGTTAACGGGGGCGGTATTCCAATAGTCGCGCCGTTCCACGGGCGCGACAGCCGCTTTTGCACCAAACCATTCtgtgt                                                                                                                                                         >  1:318678/1‑251 (MQ=255)
                                                                                        cGATCGCGCATTTGTTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCCCGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCTCCAATCCGTTCTGTGTGGttt                                                                                                                                                    <  1:582207/251‑1 (MQ=255)
                                                                                          aTCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGGTTGTCTCTATCCCCTTTGTTCGGGGGGGCGGGATTTCAATGGTCGGGCCGCTCCCCGGTCCGGGCAACCGCCTTGGCACCCATCCGCTCTGGGTGGGTTTc                                                                                                                                                  >  2:50987/1‑251 (MQ=255)
                                                                                              cgcATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTc                                                                                                                                              <  1:575846/251‑1 (MQ=255)
                                                                                              cgcATTTGGTCAGGGCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTc                                                                                                                                              <  2:21161/251‑1 (MQ=255)
                                                                                                    tGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTaaa                                                                                                                                         <  1:172889/250‑1 (MQ=255)
                                                                                                    tGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGGTTGTCTCTATCCCCTTTTTTAGCGGGGGCGGGATTCCAAAGGGGGCGCCCTTCCACGGTCCCGACAGCCCCCTTGGCACCAATCCGGTCTGTGGGGGTTTCCCCCGTaaa                                                                                                                                         >  1:693244/1‑250 (MQ=255)
                                                                                                      gTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGAt                                                                                                                                      <  2:608490/251‑1 (MQ=255)
                                                                                                                 gCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCt                                                                                                                           <  1:684497/251‑1 (MQ=255)
                                                                                                                    caTGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGtt                                                                                                                        <  1:733197/251‑1 (MQ=255)
                                                                                                                    caTGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGtt                                                                                                                        <  2:167877/251‑1 (MQ=255)
                                                                                                                        aaGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCtt                                                                                                                    <  1:51008/251‑1 (MQ=255)
                                                                                                                        aaGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCtt                                                                                                                    <  2:29170/251‑1 (MQ=255)
                                                                                                                        aaGCGAAGGAGCTGGAGATTGAGAAAGCTCAACAGCACGGTATTGCCGACGTGGCGCTACATAACTCGCATCATATCGGCTGTATCGGTTACTGGCCGGACCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCGGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGGGTGGTTTTCCCTCGGAAAGATAATTTCCCGCTGTTGCtt                                                                                                                    <  2:368563/251‑1 (MQ=255)
                                                                                                                          gCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGGTTGTTTCTTTCCACCTTGTTAGCGTGGGCGGTATTCCCATTGTCGCGCCCGTTCACGGTTGCGGCCGCCGCCTTGGCACCAAACCCTTCTGTGTGGTTTTTCCCCGTTAAGATAATTTCCCGCCGTTGCTTg                                                                                                                   >  2:726067/1‑250 (MQ=255)
                                                                                                                             aTGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGGTTGTCTCTATCCACTTTTTTAGCGGGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACCGCCGCTTTTGCACCAATCCGTTCTGTGTGGGTTTCCCTCGTAAAAATAATTTCCCGCTGTTGCTTGAtt                                                                                                                >  1:437930/1‑250 (MQ=255)
                                                                                                                                  gcTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGGTTGTCTCTATCCCCTTTGTTTGCGTGGTCGGGTTTCCAATAGTCGCGCCGTTCCACGGGCGCGGCAGCCGCTTTGGCCCCAATCCGTTCTGTGTGGTTTTTCCCCCTAAAGATAATTTCCCCCTGTTGCTTGATTaatccc                                                                                                          >  2:66825/1‑246 (MQ=255)
                                                                                                                                                aaGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTACGCCACCAGCGCCATTGCa                                                                                            <  1:22236/251‑1 (MQ=255)
                                                                                                                                                aaGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGGCGGGATTCCCATGGTCGCGCCGTTCCACGGGCGCGGCAGCCGCTTTGGCCCCAATCCGTTCTGTGTGGTTTTTCCCCGTAAAGAAAATTTCCCGCTGTTGCTTTATTACGCCACCAGCGCCATTGcc                                                                                            >  1:240105/1‑250 (MQ=255)
                                                                                                                                                      aTCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGGTTGTTTCTATCCCCTTTTTTAGCGTGGTCGGGATTTCAAAGGGCGCGCCGGTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGGTTTTCCTCGTTAAGATAATTTTCCCCTGTTGCTTGGTTACGCCCCCAGCGCcatttcattttgc                                                                                      >  2:39699/1‑248 (MQ=255)
                                                                                                                                                              ggTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGGTTGTCTCTTACCACTTTGTTAGCGTGGTCGGGATTTCAAAGGGTGGGCCGGTCTACGGTCGGGCCAGCCGCCTTGGGACCAAATCGTTCCGTGTGGTGTTTCCTCGGAAAGAAAATTTCCCGCTGTTGCTTGATTTGGCCAACACGGCCATTTGATTTTGCAAAAcccc                                                                              >  2:717776/1‑250 (MQ=255)
                                                                                                                                                                  ttGCCGCCGTGGGGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCTCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTACGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGt                                                                           <  2:318678/250‑1 (MQ=255)
                                                                                                                                                                    gccgccGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGGATTCCAATTGTCGCGCCGGTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTTTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTTCTTGATTACGCCCCCAGCGCCATTGCCTTTGGCAAAACCCGCGTCGc                                                                        >  1:684493/1‑251 (MQ=255)
                                                                                                                                                                      cgctgTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTGCTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCATTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTACGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCt                                                                      <  1:172390/247‑1 (MQ=255)
                                                                                                                                                                       gccgTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGGCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTACGCCACCAGCGCCATTGCCTTTGGCAAAACCCGCGGCGCCtt                                                                     >  1:623837/1‑250 (MQ=255)
                                                                                                                                                                         cgTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGAAGCGGCGGGGGTTTTCTATATACACTTTGTTAGAGTGGGAGGTATTCCAATATAGGCGCCGGTCCCACGCGGCGGCACTCCCGTTAGCACAAATCCCTTCTGATTGGTTTTTCCTCGTTAAAAAAATTTTCCGCTGTTGCTTGATTTAGCCCCCAGCGCCAATGCAATTGGCAAAAACCGCGGTGCCTggg                                                                   >  1:396623/1‑250 (MQ=255)
                                                                                                                                                                           tGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGGCGGGATTTCAATGGTCGCGCCGTTCCCCGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTACGCCACCAGCGCCATTGCATTTTGCAAAACCCGCGTCGCCTGGCAt                                                                 >  1:736724/1‑251 (MQ=255)
                                                                                                                                                                               gcTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGGTTGTCTCTATCCACCTTGTTTGCGTGGGCGGTATTCCAATGGTCGCGCCGTTCCACGGTCCCGACAGCCGCCTTGGCACCAATCCGTTCTGTGTGGGTTTCCCTCGTAAAGATAATTTTCCGCTGTTGCTTGATTACGCCACCAGCGCCCTTGCATTTTGCAAAACCCGCGTCGCCTGGCCTaaa                                                              >  1:333537/1‑250 (MQ=255)
                                                                                                                                                                                   cATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTACGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGt                                                         <  1:39699/251‑1 (MQ=255)
                                                                                                                                                                                    aTAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCCCTTTTGTAGCGTGGGCGGAATTTCAATGGTCGCGCCGTTCCCCGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCCCGTAAAGATAATTTCCCGCTGTTGCTTGATTACGCCACCAGCGCCCTAGCATTTGGCAAAAACCGCGTCGCCTGGCATAAAGGccgc                                                        >  1:40128/1‑248 (MQ=255)
                                                                                                                                                                                    aTAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTACGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGtn                                                        <  2:343631/251‑2 (MQ=255)
                                                                                                                                                                                    aTAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTACGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGTc                                                        <  1:377883/251‑1 (MQ=255)
                                                                                                                                                                                      accTCGCATCATAGCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTACGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGTcc                                                       <  2:592887/248‑1 (MQ=255)
                                                                                                                                                                                          cccatcatATCGGCCGTTTCGGTTACTTGGCGGATCATTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCCCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTACGCCACCAGCGCCATTTCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGTCCCCg                                                    <  2:232389/247‑1 (MQ=255)
                                                                                                                                                                                           gCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTACGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGTCCCCGTgc                                                 <  2:736724/251‑1 (MQ=255)
                                                                                                                                                                                           gCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGGTTGTCTCTATCCACTTTGTTAGCGTGGTCGGGATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTTCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTACGCCACCAGCGCCCTTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGTCCCCGTgc                                                 >  1:309449/1‑251 (MQ=255)
                                                                                                                                                                                             atcatATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGTTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTACGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGTCCCCGTgccg                                               <  1:630918/251‑1 (MQ=255)
                                                                                                                                                                                                atatCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGGTTGTCTCTATCCCCTTTGTTAGCGGGGTCGGGATTCCAATGGGCGCGCCGGTCCACGGGCGCGACAGCCGCTTTGGCACCAATCCCTTCTGTGGGGTTTTCCCCCGTAAAGATAATTTTCCGCTGTTGCTTTGTTACGCCCCCAGCGCCCTTGCCTTTGGCAAAAACCCCGTCCCCCGGCCTAAAAGCGTCCCCggtccgccc                                            >  2:349477/1‑243 (MQ=255)
                                                                                                                                                                                                 tatCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTCTCCACTTTGTTAGCGTGGTCGGTATTCCAAAGGTCGCGCCGTTTCCCGGTCGCGACAGCCGCTTTTGCACCAATCCGTTCTGTGTTGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTTCTTGCTTACTCCACCAGCGCCATTTCATTTTGCAAAAACCGCGTCGCCTGGCATAAAGGCGTCCCCGTgccgcc                                             >  1:710014/1‑249 (MQ=255)
                                                                                                                                                                                                        cGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTACGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGTCCCCGTGCCGCCAGGTTGCCt                                    <  1:66825/251‑1 (MQ=255)
                                                                                                                                                                                                             cGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTACGCCACCAGCGCTATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGTCCCCGTGCCGCCAGGTTGCCTGAtt                                <  2:65340/250‑1 (MQ=255)
                                                                                                                                                                                                               gTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGGTTGTCTCTATCCACTTTGTTAGCGTGGGCGGGATTTCAATTGTCGCGCCGTTCCACGGTCGCGACCGCCGCTTTTGCACCAATCCGTTCTGTGTGGGTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTACGCCACCAGCGCCCTTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGTCCCCGTGGCGCCAGGTTGCCTGAtttac                             >  1:175887/1‑248 (MQ=255)
                                                                                                                                                                                                                  aCTGGGCGGAGCAGTGTGCAGCGGCGGGGGTTGTCTCTATCCCCTTTGTTAGCGTGGGCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGGTTACGCCAACAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGTCCCCGTGCCGCCCGGGTGCCTGATTGACGtt                          >  1:238341/1‑251 (MQ=255)
                                                                                                                                                                                                                                       cggcggGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTACGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGTCCCCGTGCCGCCAGGTTGCCTGATTGACGTTAACGGCGTGCCGACGACCAAt     <  1:342800/251‑1 (MQ=255)
                                                                                                                                                                                                                                         gcggGGTTTGTCTCTATCCACTTTGTTAGCGGGGTCTGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTACGCCGCCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGTCCCCGTGCCGCCAGGTTGCCTGATTGACGTTAACGGCGTGCCGACGACCAATcc   <  1:636388/251‑1 (MQ=255)
                                                                                                                                                                                                                                          cggGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTACGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGTCCCCGTGCCGCCAGGTTGCCTGATTGACGTTAACGGCGTGCCGACGACCAATCcg  >  2:625088/1‑251 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                        
ATGTACGCTCCTGGAGTCAGGGGCACCTGCAAATTAACCATCATGCCAAAACCGTTAAAGAGGCGGGGGCGGCGGTCACGCTCGATGGCGATCGCGCATTTGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTACGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGTCCCCGTGCCGCCAGGTTGCCTGATTGACGTTAACGGCGTGCCGACGACCAATCCG  >  NC_000913/836517‑837001

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