Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913634,3150AG30.0% 95.2 / 13.1 50K244E (AAA→GAA) ybdLmethionine transaminase
Reads supporting (aligned to +/- strand):  ref base A (10/25);  new base G (15/0);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 2.90e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GAGCGCACCAGACTGGTGATCCTCAACACTCCGCATAACCCCAGTGCAACTGTCTGGCAGCAGGCTGATTTCGCCGCTTTGTGGCAGGCGATCGCCGGGCACGAGATTTTTGTCATTAGCGATGAAGTCTACGAGCACATCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAATACCCCGGCACAGCTGGCGCTTGCTGATATGCTACGTGCAGAACCTGAGCATTATCTTGCGTTACCGGACTTTTATCGCCAGAAGCGCGATATTCTGGTGAATGCTTTAAATGAAAGCCGGCTGGAGATTTTACCGTGTGAAGGTACATACTTTTTGCTGGTGGATTAC  >  NC_000913/634069‑634563
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                           
gagcgCACCAGACTGGTGATCCTCAACACTCCGCATAACCCCAGTGCAACTGTCTGGCAGCAGGCTGATTTCGCCGCTTTGTGGCAGGCGATCGCCGGGCACGAGATTTTTGTCATTAGCGATGAAGTCTACGAGCACATCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGt                                                                                                                                                                                                                                                         <  2:210931/251‑1 (MQ=255)
     cACCAGACTGGTGATCCTCAACACTCCGCATAACCCCAGTGCAACTGTCTGGCAGCAGGCTGATTTCGCCGCTTTGTGGCAGGCGATCGCCGGGCACGAGATTTTTGTCATTAGCGATGAAGTCTACGAGCACATCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGtt                                                                                                                                                                                                                                                    <  2:65459/251‑1 (MQ=255)
     cACCAGACTGGTGATCCTCAACACTCCGCATAACCCCAGTGCAACTGTCTGGCAGCAGGCTGATTTCGCCGCTTTGTGGCAGGCGATCGCCGGGCACGAGATTTTTGTCATTAGCGATGAAGTCTACGAGCACATCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGtt                                                                                                                                                                                                                                                    <  1:275102/251‑1 (MQ=255)
          gACTGGTGATCCTCAACACTCCGCATAACCCCAGTGCAACTGTCTGGCAGCAGGCTGATTTCGCCGCTTTGTGGCAGGCGATCGCCGGGCACGAGATTTTTGTCATTAGCGATGAAGTCTACGAGCACATCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATtgt                                                                                                                                                                                                                                               <  2:441658/251‑1 (MQ=255)
             tGGTGATCCTCAACACTCCGCATAACCCCAGTGCAACTGTCTGGCAGCAGGCTGATTTCGCCGCTTTGTGGCAGGCGATCGCCGGGCACGAGATTTTTGTCATTAGCGATGAAGTCTACGAGCACATCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGGTTTTTTCTTTTGGAAGACCTAACATTTGAACGGCTGGGAAGAGGGtttgatcggg                                                                                                                                                                                                                                            >  2:511351/1‑243 (MQ=255)
             tGGTGATCCTCAACACTCCGCATAACCCCAGTGCAACTGTCTGGCAGCAGGCTGATTTCGCCGCTTTGTGGCAGGCGATCGCCGGGCACGAGATTTTTGTCATTAGCGATGAAGTCTACGAGCACATCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGtt                                                                                                                                                                                                                                                    <  1:511351/243‑1 (MQ=255)
              ggTGATCCTCAACACTCCGCATAACCCCAGTGCAACTGTCTGGCAGCAGGCTGATTTCGCCGCTTTGTGGCAGGCGATCGCCGGGCACGAGATTTTTGTCATTAGCGATGAAGTCTACGAGCACATCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGtt                                                                                                                                                                                                                                            <  2:514726/250‑1 (MQ=255)
                  aTCCTCAACACTCCGCATAACCCCAGTGCAACTGTCTGGCAGCAGGCTGATTTCGCCGCTTTGTGGCAGGCGATCGCCGGGCACGAGATTTTTGTCATTAGCGATGAAGTCTACGAGCACATCATCTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGcc                                                                                                                                                                                                                                       <  1:355497/251‑1 (MQ=255)
                  aTCCTCAACACTCCGCATAACCCCAGTGCAACTGTCTGGCAGCAGGCTGATTTCGCCGCTTTGTGGCAGGCGATCGCCGGGCACGAGATTTTTGTCATTAGCGATGAAGTCTACGAGCACATCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGcc                                                                                                                                                                                                                                       <  2:50558/251‑1 (MQ=255)
                         acacTCCGCATAACCCCAGTGCAACTGTCTGGCAGCAGGCTGATTTCGCCGCTTTGTGGCAGGCGATCGCCGGGCACGAGATTTTTGTCATTAGCGATGAAGTCTACGAGCACATCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTTTTTATTTGGGAAGAACTCTTATATGGACGGGTGGGAAGTGGGGTCATGGGTTGCGGCAGCGccc                                                                                                                                                                                                                                >  1:504376/1‑251 (MQ=255)
                                  aTAACCCCAGTGCAACTGTCTGGCAGCAGGCTGATTTCGCCGCTTTGTGGCAGGCGATCGCCGGGCACGAGATTTTTGTCATTAGCGATGAAGTCTACGAGCACATCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGGTTTTTCTTTTTGCCAGACCCCTCATATTGCCGGCTGGGAAGGGGGGTTTTTTGTTGTGCCCGCGcccctctccgcc                                                                                                                                                                                                                       >  2:368686/1‑242 (MQ=255)
                                     aCCCCAGTGCAACTGTCTGGCAGCAGGCTGATTTCGCCGCTTTGTGGCAGGCGATCGCCGGGCACGAGATTTTTGTCATTAGCGATGAAGTCTACGAGCACATCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTTTTCATTTTGGAAGACCCCTCATATGACCCGCTGGAAAGTGGGGTTTTGTGttttggccgggcccctcaacgcccg                                                                                                                                                                                                                     >  1:128520/1‑227 (MQ=255)
                                                       ggCAGCAGGCTGATTTCGCCGCTTTGTGGCAGGCGATCGCCGGGCACGAGATTTTTGTCATTAGCGATGAAGTCTACGAGCACATCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCAAAATTCGCAAGGTACATCAg                                                                                                                                                                                                 <  2:378833/251‑1 (MQ=255)
                                                                   aTTTCGCCGCTTTGTGGCAGGCGATCGCCGGGCACGAGATTTTTGTCATTAGCGATGAAGTCTACGAGCACATCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCttt                                                                                                                                                                                   <  2:129166/251‑1 (MQ=255)
                                                                       cgccgcTTTGTGGCAGGCGATCGCCGGGCACGAGATTTTTGTCATTAGCGATGAAGTCTACGAGCACATCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTTTTGCGCCAGCGCCCATCAgc                                                                                                                                                                                                                          <  1:330138/211‑1 (MQ=255)
                                                                       cgccgcTTTGTGGCAGGCGATCGCCGGGCACGAGATTTTTGTCATTAGCGATGAAGTCTACGAGCACATCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCgg                                                                                                                                                                               <  1:494217/251‑1 (MQ=255)
                                                                         ccgcTTTGTGGCAGGCGATCGCCGGGCACGAGATTTTTGTCATTAGCGATGAAGTCTACGAGCACATCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGCTTCTTCATTTGGCAAGACCTCTCATATGAACGGCTGGAAAGGGGGGTATTGGGTTGGGCCAGCGCCCATCAGCGCCGAAAATCGCAAGGGAAATCAGTATCTGACCCTTTCgggg                                                                                                                                                                             >  1:471772/1‑249 (MQ=255)
                                                                           gcTTTGTGGCAGGCGATCGCCGGGCACGAGATTTTTGTCATTAGCGATGAAGTCTACGAGCACATCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGaa                                                                                                                                                                           <  2:476037/251‑1 (MQ=255)
                                                                                  ggCAGGCGATCGCCGGGCACGAGATTTTTGTCATTAGCGATGAAGTCTACGAGCACATCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAATAcccc                                                                                                                                                                     <  1:455867/250‑1 (MQ=255)
                                                                                            cGCCGGGCACGAGATTTTTGTCATTAGCGATGAAGTCTACGAGCACATCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAATACCCCGGCACAGCTg                                                                                                                                                           <  2:202312/250‑1 (MQ=255)
                                                                                             gCCGGGCACGAGATTTTTGTCATTAGCGATGAAGTCTACGAGCACATCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTTTTCATTTTGCAAGACCTATCATATGAACGGCTGGGAAGGGGGGTATTTTTTTGCGCCAGCGCCCCTCAGCGCCCAAAATCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAATACCCCGGCACAggctgg                                                                                                                                                         >  2:139283/1‑246 (MQ=255)
                                                                                               cGGGCACGAGATTTTTGTCATTAGCGATGAAGTCTACGAGCACATCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCCAGACCTATCATATGACCGGCTGGGAAGTGGGTTTTTGTGTTGCGCCAGCGcccctc                                                                                                                                                                                                                             >  1:86498/1‑181 (MQ=255)
                                                                                               cGGGCACGAGATTTTTGTCATTAGCGATGAAGTCTACGAGCACATCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATc                                                                                                                                                                                                                             <  2:86498/184‑1 (MQ=255)
                                                                                                  gCACGAGATTTTTGTCATTAGCGATGAAGTCTACGAGCACATCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTTTTCATTTGGCAAGACCTATCAAATGACCGGCTGGGAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAAATTGCTAGGTACATCAGGATACGACCTTTTCGGTGCATATACCGGCACAgactgcgcttg                                                                                                                                                    >  1:104940/1‑241 (MQ=255)
                                                                                                   cACGAGATTTTTGTCATTAGCGATGAAGTCTACGAGCACATCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTTATTTTGCAAGACCTATCATATGACCGGGTGGGAAGAGGGGTATTGTGTTTCGCCAGCGCCCATCCGCGCCGAAATTCGCAAGGTACATCAGGATCTCACCCTTTTCGTGAAAACCCCGGCACACCTGGCGCTTGc                                                                                                                                                   >  2:68713/1‑251 (MQ=255)
                                                                                                    aCGAGATTTTTGTCATTAGCGATGAAGTCTACGAGCACATCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTTGCCAGACCTATCATATGGCCGGGTGGAAAGTGGGGTATTGTGTTGCGCCAGCGCCCATCCGCGCCGGAAATTGCCAGGTAAATCAGGATCTTGCCTTTTTGGTGGATAACCCCGCACAGCCGGCGCTTGcc                                                                                                                                                  >  2:165441/1‑250 (MQ=255)
                                                                                                         aTTTTTGTCATTAGCGATGAAGTCTACGAGCACATCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTTGCCAGACCTATCATATGACCGGCTGGGAAGGGGGGTATTGTGTTGCTCCCACGCCCATCAACGCCGAAATTCGCAAGGGACATCAGTATATGACCTTTTCGGGGAATAACCCCGCCCAGCTGGCGCTTGCTGatat                                                                                                                                             >  1:202186/1‑251 (MQ=255)
                                                                                                             ttGTCATTAGCGATGAAGTCTACGAGCACATCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGATTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGGAAGGGGGGTATTGTGTTGCGCCAACGCCCATCAGCGCCGGAATTTGCAAGGTACATCAGTATCTGACCTTTTCGGTGGAAACCCCGGCACAACTGGCGCTTTCTGATATGCtc                                                                                                                                         >  1:56072/1‑250 (MQ=255)
                                                                                                              tGTCATTAGCGATGAAGTCTACGAGCACATCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTTGCAAGACCTATCATATGACCGGCTGGGAAGGGGGTTATTGTGTTGCGCCCGCGCCCCTCAGCGCCGAAATTCGCAAGGGACATCAGTATCTGACCTTTTCGGTGAATACCCCCGCACAACTGGCGCTTGCTGATATTCTAc                                                                                                                                        >  1:512145/1‑251 (MQ=255)
                                                                                                                           gAAGTCTACGAGCACATCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGGAAGACCTATCATATGACCGGCTGGAAAGTGGGGTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGGACATCAGTATCTTACCTTTTCGGTGAATACCCCGGCACAGCTGGCGCTTGCTGATATTCTACGTGCAGAACCTg                                                                                                                            >  1:113611/1‑250 (MQ=255)
                                                                                                                                        acaTCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGGTTTTTCATTTGGCAAAACCTCTCATATGACCGGCTGGGAAAGGGGGTATTTTGGTGCGCCCGCCCCCACCCGCGCCGAAATTTGCAAGGTACATCCGTAACTGGCCCTTTCCGTGAAAAACCCGGCACAGCTGGCGCCTGCTGATATTCTACTGGCCGAACCTGGACAATATCTTGCg                                                                                                              >  2:377482/1‑251 (MQ=255)
                                                                                                                                            cAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCCTTTTGCCAAACCTATCATATGACCCGCTGGGAAGGGGGGTGTTGTGTTGCGCCCGCGCCCCTCCGCGCCCCAATTCCCCAGGGACATCCGTATTCTACCTTTTTGGTGTATACCCCCGCACAGCCGGCGCTTGCTTGTATTCCACATGCAGAACCCGAGCACTATCTTTCGtttc                                                                                                          >  2:463862/1‑249 (MQ=255)
                                                                                                                                            cAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAATACCCCGGCACAGCTGGCGCTTGCTGATATGCTACGTGCAGAACCTGAGCATTATCTTGCGTTa                                                                                                           <  1:165441/250‑1 (MQ=255)
                                                                                                                                            cAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAATACCCCGGCACAGCTGGCGCTTGCTGATATGCTACGTGCAGAACCTGAGCATTATCTTGCGTTAc                                                                                                          <  1:157082/251‑1 (MQ=255)
                                                                                                                                                       aacaGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGGAAGGGGGGTATTGTGTTGCGCCAGCGCCCATCAGCGCCGCAATTCGCAAGGGACATCAGTATCTGACCCTTTCGGGGAATACCCCGGCACAACTGGGGCTTGCTGATTTGCTACGTGCAGAACCTGAGCATTATCTTGCGTTACCGGACTTTTAt                                                                                               >  1:201370/1‑251 (MQ=255)
                                                                                                                                                          aGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCCTTTTGCCAGACCTATCCTATGACCGGCGGGGAAGGGGGGTATTGTGTTGCGCCCGCGCCCATCAACGGCGAAAATTGGAAGGTACATCAGTATATGACCTTTTTGGTGAATACCCCGGCACAGCAGGGGCTTGCTGTTATGCTACGTGCCGAACCTTAGCATTATTTTGTGTTTCCCGACTTTTTTCGc                                                                                            >  2:118350/1‑251 (MQ=255)
                                                                                                                                                           gggCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGTTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAATACCCCGGCACAGCTGGCGCTTGCTGATATGCTACGTGCAGAACCTGAGCATTATCTTGCGTTACCGGACTTTTATCGcc                                                                                           >  1:55048/1‑251 (MQ=255)
                                                                                                                                                                  gCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGTTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAATACCCCGGCACAGCTGGCGCTTGCTGATATGCTACGTGCAGAACCTGAGCATTATCTTGCGTTACCGGACTTTTATCGCCAGAAgcg                                                                                    <  2:504376/251‑1 (MQ=255)
                                                                                                                                                                        gtgCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGGAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAATACCCCGGCACAGCTGGCGCTTGCTGATATGCTACGTGCAg                                                                                                                                  >  1:86474/1‑199 (MQ=255)
                                                                                                                                                                        gtgCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAATACCCCGGCACAGCTGGCGCTTGCTGATATGCTACGTGCAg                                                                                                                                  <  2:86474/199‑1 (MQ=255)
                                                                                                                                                                          gCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCCTTTTGCAAGACCTATCATATGACCGGGTGGGAAGGGGGGTATTGTGTTGCGCCCGCGCCCATCAACGCCGAAAATCGCAAGGGACATCAGGATCTGACCCTTTTGGTGAATACCCCCGCACAGCCGGGGCCTGCTGATATGCTAAGTGCAGAAACTGAGCATTTTCTTGCGTTTCCTGGCTTTTATCGcccagagcgcgggattt                                                                            >  2:378976/1‑236 (MQ=255)
                                                                                                                                                                                  aTCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAATACCCCGGCACAGCTGGCGCTTGCTGATATGCTACGTGCAGAACCTGAGCATTATCTTGCGTTACCGGACTTTTATCGCCAGAAGCGCGATATTCTGGTGAAt                                                                    <  2:362362/251‑1 (MQ=255)
                                                                                                                                                                                         gCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGTCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAATACCCCGGCACAGCTGGCGCTTGCTGATATGCTACGTGCAGAACCTGAGCATTATCTTGCGTTACCGGACTTTTATCGCCAGAAGCGCGATATTCTGGTGAATGCTTTaa                                                             <  2:43252/251‑1 (MQ=255)
                                                                                                                                                                                                               gTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAATACCCCGGCACAGCTGGCGCTTGCTGATATGCTACGTGCAGAACCTGAGCATTATCTTGCGTTACCGGACTTTTATCGCCAGAAGCGCGATATTCTGGTGAATGCTTTAAATGAAAGCCGGCTGGAGAtttt                                       <  2:543477/251‑1 (MQ=255)
                                                                                                                                                                                                                     ttcATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAATACCCCGGCACAGCTGGCGCTTGCTGATATGCTACGTGCAGAACCTGAGCATTATCTTGCGTTACCGGACTTTTATCGCCAGAAGCGCGATATTCTGGTGAATGCTTTAAATGAAAGCCGGCTGGAGATTTTACCgt                                  >  1:322093/1‑251 (MQ=255)
                                                                                                                                                                                                                           tGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAATACCCCGGCACAGCTGGCGCTTGCTGATATGCTACGTGCAGAACCTGAGCATTATCTTGCGTTACCGGACTTTTATCGCCAGAAGCGCGATATTCTGGTGAATGCTTTAAATGAAAGCCGGCTGGAGATTTTACCGTGTGa                              >  1:470525/1‑249 (MQ=255)
                                                                                                                                                                                                                                      aTCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAATACCCCGGCACAGCTGGCGCTTGCTGATATGCTACGTGCAGAACCTGAGCATTATCTTGCGTTACCGGACTTTTATCGCCAGAAGCGCGATATTCTGGTGAATGCTTTAAATGAAAGCCGGCTGGAGATTTTACCGTGTGAAGGTACATACtt                  <  1:225493/250‑1 (MQ=255)
                                                                                                                                                                                                                                        cATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAATACCCCGGCACAGCTGGCGCTTGCTGATATGCTACGTGCAGAACCTGAGCATTATCTTGCGTTACCGGACTTTTATCGCCAGAAGCGCGATATTCTGGTGAATGCTTTAAATGAAAGCCGGCTGGAGATTTTACCGTGTGAAGGTACATACttttt               <  2:201370/251‑1 (MQ=255)
                                                                                                                                                                                                                                          tatGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAATACCCCGGCACAGCTGGCGCTTGCTGATATGCTACGTGCAGAACCTGAGCATTATCTTGCGTTACCGGACTTTTATCGCCAGAAGCGCGATATTCTGGTGAATGCTTTAAATGAAAGCCGGCTGGAGATTTTACCGTGTGAAGGTACATACTTTTTGc             >  1:276252/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      gAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAATACCCCGGCACAGCTGGCGCTTGCTGATATGCTACGTGCAGAACCTGAGCATTATCTTGCGTTACCGGACTTTTATCGCCAGAAGCGCGATATTCTGGTGAATGCTTTAAATGAAAGCCGGCTGGCGATTTTACCGTGTGAAGGTACATACTTTTTGCTGGTGGATTAc  >  2:545947/1‑250 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                           
GAGCGCACCAGACTGGTGATCCTCAACACTCCGCATAACCCCAGTGCAACTGTCTGGCAGCAGGCTGATTTCGCCGCTTTGTGGCAGGCGATCGCCGGGCACGAGATTTTTGTCATTAGCGATGAAGTCTACGAGCACATCAACTTTTCACAACAGGGCCATGCCAGTGTGCTGGCGCATCCGCAGCTGCGTGAGCGGGCAGTGGCGGTTTCTTCATTTGGCAAGACCTATCATATGACCGGCTGGAAAGTGGGTTATTGTGTTGCGCCAGCGCCCATCAGCGCCGAAATTCGCAAGGTACATCAGTATCTGACCTTTTCGGTGAATACCCCGGCACAGCTGGCGCTTGCTGATATGCTACGTGCAGAACCTGAGCATTATCTTGCGTTACCGGACTTTTATCGCCAGAAGCGCGATATTCTGGTGAATGCTTTAAATGAAAGCCGGCTGGAGATTTTACCGTGTGAAGGTACATACTTTTTGCTGGTGGATTAC  >  NC_000913/634069‑634563

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