Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,388,202 A→T F309I (TTT→ATT)  aaeA ← aromatic carboxylic acid efflux pump membrane fusion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,388,2020AT96.2% 158.5 / NA 54F309I (TTT→ATT) aaeAaromatic carboxylic acid efflux pump membrane fusion protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  major base T (24/27);  minor base G (0/2);  total (25/29)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGG‑CCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATAC‑GAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGCGTTGGTGACCCCTGCGGCAACACTATC  >  NC_000913/3387955‑3388448
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        
aCAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCCGTTTTACCGCACAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAtctc                                                                                                                                                                                                                                                       >  1:358656/1‑251 (MQ=255)
             gtgCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTGAGCAATGGAGAAAATTCCCATCGTGATTAACCAATCTCACGCAGGCGATgg                                                                                                                                                                                                                                          <  1:214430/251‑1 (MQ=255)
                 cGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAATCTCACGCAGGCGATGGG‑CCa                                                                                                                                                                                                                                     <  2:360658/251‑1 (MQ=255)
                     gatgCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAATCTCACGCAGGCGATGGG‑CCAtttt                                                                                                                                                                                                                                 <  1:220236/251‑1 (MQ=255)
                           cAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAATCTCACGCAGGCGATGGG‑CCATTTTACGgaa                                                                                                                                                                                                                           <  2:358656/251‑1 (MQ=255)
                           cAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAATCTCACGCAGGCGATGGG‑CCATTTTACGga                                                                                                                                                                                                                            >  1:301746/1‑250 (MQ=255)
                              aaaGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAATCTCACGCAGGCGATGGG‑CCATTTTACGgaagaa                                                                                                                                                                                                                        <  2:323416/251‑1 (MQ=255)
                                         aGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCCGTTTTTCCGCAAAGCGAATATGTTGGGTAGCAAAGGAGAAAATAACCATAGTGATTAACCAATCTCACACAGGCGGTGGG‑CCATTTTACGGAAGAAAGAATCCTGGc                                                                                                                                                                                                             >  2:297952/1‑251 (MQ=255)
                                            gAATAGCGCCAGAATACGGTTCACCTCCCGCAGCCAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTAACGCAAAGCGGATATGTTGGTTAGCAATGGGGAAAATACCCAACGGGATTTACCAATCTAAAGCCAGAGAAGGG‑CCATTTTACGGAAGAACGAATCCTGGCttt                                                                                                                                                                                                          >  1:296413/1‑251 (MQ=255)
                                            gAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGGTGGTTAGCAATGGAGAAAATACCCATCGGGATTAACCAATCTCACGCAGGCGATGGG‑CCATTTTACGGAAGAACGAATCCTGGCttt                                                                                                                                                                                                          <  1:37284/251‑1 (MQ=255)
                                             aaTAGCCCCAGAATTCGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAGCTCACGCAGGCGATGGG‑CCATTTTACTGAAGAACGAAGCCGGGCTTtc                                                                                                                                                                                                         <  2:42978/251‑1 (MQ=255)
                                                aGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAATCTCACGCAGGCGATGGG‑CCATTTTACGGAAGAACGAATCCTGGCTTtcgtc                                                                                                                                                                                                      <  2:222823/251‑1 (MQ=255)
                                                  cgcCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAATCTCACGCAGGCGATGGG‑CCATTTTACGGAAGAACGAATCCTGGCTTTCGTcgc                                                                                                                                                                                                    <  2:117116/251‑1 (MQ=255)
                                                  cgcCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTACAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGGTTTTCCGCAAAGCGAAAATGTTGGGTTGCAATGGTGAAAATACCCATCGTGATTAAACAATCTCACCCAGGCGGTGGG‑CCCTTTTTCGGAAGCACGAATCCTGGCTTTTGTcgc                                                                                                                                                                                                    >  1:265273/1‑251 (MQ=255)
                                                                   ccTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGGCAGGTTTACCCCAAAGCGAAAATGTTGGTTTGCCATGGAGAAAATACCCATCGTGATTAAACAATCTCACGCAGGGGATGGGGCCTTTTTCGGGAAGA‑CGAATCCCGGCTTTCGTCGCGCTATTGTTTGTCAGTg                                                                                                                                                                                   >  2:141204/1‑251 (MQ=255)
                                                                     tCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTTGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAATCTCACGCAGGCGAGGGG‑CCATTTTACGGAAAAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGGGac                                                                                                                                                                                 >  1:294995/1‑251 (MQ=255)
                                                                            gcaAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAATCTCACGCAGGCGATGGG‑CCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACt                                                                                                                                                                          <  1:113259/251‑1 (MQ=255)
                                                                              aaaGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCCGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAATCTCACGCAGGCGATGGG‑CCATTTTACGGAAGAAAGAATCCTGGCTTTCGTGTCGATCTTGTTTGcccgggaacacccccgt                                                                                                                                                                        >  1:413992/1‑236 (MQ=255)
                                                                                                      cgccgcTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAATCTCACGCAGGCGATGGG‑CCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGAt                                                                                                                                                <  1:193305/251‑1 (MQ=255)
                                                                                                          gcTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAATCTCACGCAGGCGATGGG‑CCATTTTACGGAAGAAAGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCCGGCCAGATGttt                                                                                                                                            >  1:374765/1‑250 (MQ=255)
                                                                                                               cAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAATCTCACGCAGGCGATGGG‑CCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCtt                                                                                                                                        <  2:391764/250‑1 (MQ=255)
                                                                                                               cAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAATCTCACGCAGGCGATGGG‑CCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCtt                                                                                                                                        <  1:337834/250‑1 (MQ=255)
                                                                                                               cAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAATCTCACGCAGGCGATGGG‑CCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTg                                                                                                                                       <  2:76055/251‑1 (MQ=255)
                                                                                                                               ggCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAATCTCACGCAGGCGATGGG‑CCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACg                                                                                                                        <  1:52265/250‑1 (MQ=255)
                                                                                                                                gCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAATCTCACGCAGGCGATGGG‑CCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGaa                                                                                                                      <  1:174419/251‑1 (MQ=255)
                                                                                                                                 cGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAATCTCACGCAGGCGATTGG‑CCATTTTACGGAAGAACGAATCCTGGCTGTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGTTTGTCGAGACGAAt                                                                                                                     <  2:395747/251‑1 (MQ=255)
                                                                                                                                       cAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAATCTCACGCAGGCGATGGG‑CCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGGGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATAC‑gcac                                                                                                              >  2:115139/1‑248 (MQ=255)
                                                                                                                                        aGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAATCTCACGCAGGCGATGGG‑CCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATAC‑GAAcc                                                                                                             <  2:89396/251‑1 (MQ=255)
                                                                                                                                             gAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAATCTCACGCAGGCGATGGG‑CCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATAC‑GAACCGGAac                                                                                                        <  2:301746/251‑1 (MQ=255)
                                                                                                                                                gTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAATTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAATCTCACGCAGGCGATGGG‑CCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATAC‑GAACCTGAACACg                                                                                                     <  1:83681/251‑1 (MQ=255)
                                                                                                                                                   aaaGCCAACAAACAGCGCCAGTACAATGGCGGTCACCCGTTTTACAGCAAAGAGAATCTGTTGGGTAGCCATGGAGAAAATACCCATCGTGATTAACCAATCTCACGCAGGCGATGGG‑CCATTTTACGGAAGAACGAATCCTTGCTTTCGTCGAGCTATTGTTTGCCAGTGACCACCCATGTCGCAGTGGTGCCCGCAGGACAGATTTTCTCTTGCTGGTTGTCGAGACGACTACCGAACCgaaccccttt                                                                                                   >  2:269074/1‑242 (MQ=255)
                                                                                                                                                                 gcgcCAGTACAATGTCGTTCGCCAGTTTTACCGCAAAGCGAATATTTTGTTTAGCAATGGAGAAAATACCCATCGTGATTAACCAATCTCACGCAGGCGATGGC‑CCATTTTACGGAAGAACGATTCCTGGCTTTCGTCGCGATCTTTTTTGCCAGTGACCACCACTGTATCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATAC‑GAACCGGAACACGTTGCGCAAGACGCAccc                                                                                    <  1:18215/251‑1 (MQ=255)
                                                                                                                                                                    ccAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAATCTCACGCAGGCGATGGG‑CCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATAC‑GAACCGGAACACGTTGCGCAAGACGCACCCAtt                                                                                 <  1:313851/251‑1 (MQ=255)
                                                                                                                                                                                  gTCGCCAGTTTTACCGCAACGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAATCTCACGCAGGCGATGGG‑CCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATAC‑GAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCt                                                                   >  2:427174/1‑251 (MQ=255)
                                                                                                                                                                                    cGCCAGTTTTAACGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCAATCGTGATTAACCAATCTCACGCCGGCGATGGG‑CCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGACAGTGACCACCACTGTAGCAGTGGTGCACGCAGGCCCGATGTTCTCTTGCTGGTTGTCGAGACGCATAC‑GAACCGGCACACGTTGCGCAAGACGCCCCCATTCAAGGTTAGCGTCtat                                                                 >  2:238780/1‑251 (MQ=255)
                                                                                                                                                                                       cAGTTTTACCGCAAAGCGAATATGTTGGTTCGCAATGGAGACAATACCCCTCGCGATTCACCAATCTCACGCAGGCGCTGGG‑CCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTTACCACCACTGTAGCAGTGGCGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGCATAC‑GAACCGGAACACGTTGCGCAAGACGCACCCCTTCAAGGTTAGAGTCTATAGt                                                              >  2:192154/1‑251 (MQ=255)
                                                                                                                                                                                       cAGTTTTACCGAAAAGCGAATATGTTGGTTAGCAATGAAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGCTGGG‑CCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATAC‑GAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTCTAGt                                                              >  2:392043/1‑251 (MQ=255)
                                                                                                                                                                                             tACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAATCTCACGCAGGCGATGGG‑CCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATAC‑GAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCAt                                                        >  1:199750/1‑251 (MQ=255)
                                                                                                                                                                                                cGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCACCGTGATTAACCAATCTCACGCAGGCGATGGG‑CCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATAC‑GAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATccc                                                     <  1:313812/251‑1 (MQ=255)
                                                                                                                                                                                                 gCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAATCTCACGCAGGCGATGGG‑CCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATAC‑GAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCt                                                    >  1:352310/1‑251 (MQ=255)
                                                                                                                                                                                                      gCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAATCTCACGCAGGCGATGGG‑CCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATAC‑GAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTgtc                                               >  1:344751/1‑251 (MQ=255)
                                                                                                                                                                                                          atatGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAATCTCACGCAGGCGATGTG‑CCATTTTACGGAAGAACGAAGCCTGGCTTTCGTCGCGATCGTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTGCTCTTGCTGGTTGTCGAGACGAATAC‑GAACCGGAACACGTTGCGCAAGACGCTCCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCa                                           <  1:423034/251‑1 (MQ=255)
                                                                                                                                                                                                            atGTTGGGTAGCAATGGATAAAATCACCACCGTGATTAACCAATCTCACGCAGGCGATGGG‑CCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGTTGGTTGTCGAGACGAATAC‑GAACCGGAACACGTTGCGCAAGACGCCCCCATTCAAGGTGAGAGTCTATAGTCGCCATCCCTTTGTCGTCAc                                          <  2:66004/250‑1 (MQ=255)
                                                                                                                                                                                                             tGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAATCTCACGCAGGCGATGGG‑CCATTTTACGGAAGAACGAATCATGGCTTTCGTCGCGCTCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATAC‑GAACCGGAACACGATGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCAcgc                                        >  1:276560/1‑251 (MQ=255)
                                                                                                                                                                                                               ttGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAATCTCACGCAGGCGATGGG‑CCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATAC‑GAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGt                                      >  1:188629/1‑251 (MQ=255)
                                                                                                                                                                                                                  gTTATCAATGGAGAAAATACCCATCGTGATTAACCAATCTCACGCAGGCGATGGG‑CCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATAC‑GAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGtgct                                   <  2:11064/251‑1 (MQ=255)
                                                                                                                                                                                                                   ttAGCCATGGAGACAATACCCATCGTGATTAACCAATCTCACGCAGGCGCTGGG‑ACATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCAGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATAC‑GAACCGGCACACGTTGCGCAAGACGCACCCATTCAAGGATAGAGACTATAACCGCCATCCCTTTGTCGTCCCGCGtgctg                                  >  2:278651/1‑251 (MQ=255)
                                                                                                                                                                                                                    tAGCAATGGAGAAAATACCCATCGTGATTAACCAATCTCACGCAGGCGATGGG‑CCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGCCGAATAC‑GAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGtgctgc                                 >  1:15778/1‑251 (MQ=255)
                                                                                                                                                                                                                          tGGAGAAAATACCCATCGTGATTAACCAATCTCACGCAGGCGATGGG‑CCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATAC‑GAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGCGt                           >  1:397039/1‑251 (MQ=255)
                                                                                                                                                                                                                           ggAGAAAATACCCATCGTGATTAACCAATCTCACGCAGGCGATGGG‑CCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATAC‑GAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGCGtt                          >  2:359787/1‑251 (MQ=255)
                                                                                                                                                                                                                               aaaaTACCCAGCGTGATTAACCAAGCTCACGCAGGCGATGGG‑TCATTTTACGGAAGAACGAATCCGGGCTTTCGTCGCGATCTGGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGGTCTCTTGCTGGTTGTCGATACGAATAC‑GAACCGGAACACGTTGCTCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGCGTTGGTg                      <  2:323576/251‑1 (MQ=255)
                                                                                                                                                                                                                                     cccATCGTGATTAACCAATCTCACGCAGGCGATGGG‑CCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATAC‑GAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGCGTTGGTGACCCCt                >  1:333912/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 aaCCAATCTCACGCAGGCGATGGG‑CCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCGGGCCAGATGTTCTCTGGCTGGTTGTCGAGACGAATAC‑GAACCGGAACACGTTGCGCAAGACGCGCCCATTCAAGGGTAGAGTCTATAGTCGCCAGCCCATTGTCGGCACGCGTGCTGCTGGCGTTGGTGACCCCTGCGGCAACACt     <  2:265273/250‑1 (MQ=255)
                                                                                                                                                                                                                                                 aaCCAATCTCACGCAGGCGATGGG‑CCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATAC‑GAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGCGTTGGTGACCCCTGCGGCAACACt     >  1:125759/1‑250 (MQ=255)
                                                                                                                                                                                                                                                   ccAATCTCACGCAGGCGATGGG‑CCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTTACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATAC‑GAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGCGTTGGTGACCCCTGCGGCAACACTATc  <  1:67174/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        
ACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGG‑CCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCACTGTAGCAGTGGTGCCCGCAGGCCAGATGTTCTCTTGCTGGTTGTCGAGACGAATAC‑GAACCGGAACACGTTGCGCAAGACGCACCCATTCAAGGTTAGAGTCTATAGTCGCCATCCCTTTGTCGTCACGCGTGCTGCTGGCGTTGGTGACCCCTGCGGCAACACTATC  >  NC_000913/3387955‑3388448

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