Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 78,520 C→A F300L (TTC→TTA)  setA → sugar exporter SetA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091378,5200CA98.5% 191.6 / NA 67F300L (TTC→TTA) setAsugar exporter SetA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (27/38);  minor base G (1/0);  total (29/38)
Fisher's exact test for biased strand distribution p-value = 4.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

ATTTGTCGCCTCGACGTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGACT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTCTTTAATAGCCGTATGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTC  >  NC_000913/78280‑78756
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                             
aTTTGTCGCCTCGACGTTAATGTGGTCCTTCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGcc                                                                                                                                                                                                                                     <  1:724043/251‑1 (MQ=255)
aTTTGTCGCCTCGACGTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGcc                                                                                                                                                                                                                                     <  1:826818/251‑1 (MQ=255)
     gcGCCTCGACGTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTTGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGtctg                                                                                                                                                                                                                                <  1:26311/250‑4 (MQ=255)
      cGCCTCGACGTTAATGTGGACCTGCAAAACCATGTACATTATTGATCTGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGGTGGGGACCGACGCTGGGCT‑GGAAATACCAGCAATGGTTCTGGGTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAAAGATGGTCATAGCAGTGGCGGCAGGGG‑GACTGTTTTAAA‑CCGG‑GTTGTTTTTATTTTATAGCCGtcttg                                                                                                                                                                                                                               >  1:365023/1‑247 (MQ=255)
        ccgcggcGTTAATGTGGACCTTCAACACCATGTACATTATTGATATGCCGTTGGGGATCAGGAGCGAGGTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGGCGGCAGCTGGTCT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCg                                                                                                                                                                                                                             <  1:115374/245‑1 (MQ=255)
           cGACGTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑tt                                                                                                                                                                                                                          <  2:895513/250‑1 (MQ=255)
                     gtgtgcctgcactcccATGTACATTATTGATTTGCCGTTGTGTATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGGCT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGc                                                                                                                                                                                                               <  2:805934/237‑1 (MQ=255)
                        ggcctgcaaccccattgtcaatattGAAATGGCGGTGTGGGTCAGTAGCGAGGTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAAc                                                                                                                                                                                                            <  1:855559/230‑1 (MQ=255)
                             gCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCTGGTTTCCTGATGGGGACGGCAGCTGGTCT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACtttt                                                                                                                                                                                                        <  1:895058/250‑1 (MQ=255)
                             gCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGAGGGGGACGGCAGCTGGTCTCGGAAA‑ACCAGCAATGATTCTCGCCAGCTACTATGTCAAAAGATTTGGTACGCGGCGAATGGGGGGCATAGACGTGGCGGGAGGAGGTA‑TATTTTACAACCGG‑‑TTGATTTTATTTCCTAGCCCTCTGGCG‑Tttatgtctctgtgacttttt                                                                                                                                                                                                       >  1:947285/1‑236 (MQ=255)
                                acacCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAAc                                                                                                                                                                                                    <  1:346017/251‑1 (MQ=255)
                                     aTGTACATTATTGATATGCCGTGTTGGATCAGTAGCGAGTGAGGATTGCCAGACAAACTGGCTGGTTTCCTGATGGGGACGGCAGCTTGTCT‑GGGAATACCAGCAATGATTCTGTCTGGCTACTAGGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGGAATTTCGGCGGAAG‑TACTGTTTTGCA‑CCGG‑ATTGATGTTAGTGCATAGCCGGCTGACG‑TTGATGACGCTGCAATTTTTTAACGCTGt                                                                                                                                                                                               <  2:503253/251‑1 (MQ=255)
                                             tagtGATATGCCGTTGTGGATCAGTAGCGAGTTAGGAGTTCCAGACAAATTGTCGGGTTTCCTGATGGGGAGGGCAGCTGGTCT‑GGAAATACCAGTAATGATTCTGGCGTGCTACGATGTCATACGTTATGGTAAGCGGCGAATGATTGTCAGAGTAGTGGCTGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGAGTGTATTTCAGAGCCGGCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTTTATTTATc                                                                                                                                                                                        <  2:222350/247‑1 (MQ=255)
                                                 gATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCAt                                                                                                                                                                                   <  1:691682/251‑1 (MQ=255)
                                                    atGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCAttg                                                                                                                                                                                 <  1:445877/250‑1 (MQ=255)
                                                           tgtgGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑GACTGTTTTTCA‑CCGG‑ATTGATTTTATTTTATAGCCGCCTGGCG‑TTGATGACGGTGCCACTTTTTAACGCTGTATTTTTCGGGATTTTTGCgggg                                                                                                                                                                         >  1:617770/1‑250 (MQ=255)
                                                            gtgGATCATTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTGCCTGATGGGGACGGCAGCTTGTCT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATTGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTa                                                                                                                                                                        <  2:852372/251‑1 (MQ=255)
                                                              ggATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATATCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCACTGTTACGGGTAt                                                                                                                                                                       <  2:427912/250‑1 (MQ=255)
                                                                     tAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGGG‑TACTGTTTTTCA‑CCGG‑GTTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTTACGCTGTATTTATCGGCATTGTTGCGGGGATTGGGATgg                                                                                                                                                               >  1:898372/1‑250 (MQ=255)
                                                                          aGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTAAA‑CCGGGTTTGATTTT‑TTTCATAGCCGTCTGGCG‑TTGATGGCGCTTCAACTTTTTTACGCTGTTTTTTTCGGCATTGTTGGGGGGAtttgggagctttgg                                                                                                                                                          >  1:185990/1‑239 (MQ=255)
                                                                               ggATTGCCAGACAAACTGGCGGGTTTCCTTATGGGGACGGCAGCTGGTCT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGGAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑GACTTTTTTACA‑CCGG‑ATTGATTTTATTTTATAGCCGTCTGGCGGTTG‑TGGCGCTGCAACTTTTTAACGCTGTATTTTTCGGCATTGTTGCGGGGATTTGGATGCTATGGTTTc                                                                                                                                                      >  1:898703/1‑250 (MQ=255)
                                                                                gATTGCCAGACAAACTGGCGGGTTCCCTGATGGGGACGGCAGCTGGTCT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAg                                                                                                                                                    <  1:206560/251‑1 (MQ=255)
                                                                                            aaCTGGCGGGTTTCCTGATGGGGACGGCAGCTGGTCT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCt                                                                                                                                        <  2:898372/251‑1 (MQ=255)
                                                                                                  cGGGTTTCCTGATGGGGACGGCAGCTGGTCT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAaga                                                                                                                                  <  1:567028/251‑1 (MQ=255)
                                                                                                       ttCCTGATGGGGACGGCAGCTGGTCT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCTTTGTTGCGGGTATTGGGGTGCTATGGTTTCAGGATTTAATGCCTGGAAGGGCggg                                                                                                                             >  2:411814/1‑251 (MQ=255)
                                                                                                           tGATGGGGAAGGCAGCTGGTCT‑GAAAATACCAGCAATGATTATGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGTG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTAATAACTGTCTGGAG‑TTGGTGACGGTGCAAATTTTTAACGCTGTTTTTATCGGCATTGTTGAGGGGATTGGGATGCTATGGTTTCAGGAGTTCATGCCTGGAAGGGGGGGGGc                                                                                                                          >  2:913607/1‑250 (MQ=255)
                                                                                                            gATGGGGACGGCAGCTGGTCT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAAAGTTGTATTTATCGGATTTGTTGCGGGGATTGGGATGGTATGGTTTCAGGATTTAATGCCTGGAAGAGCggggga                                                                                                                          >  1:698149/1‑248 (MQ=255)
                                                                                                               ggggatggccggtgaact‑gtaAATACCTGAAATGATTGTGGCTGGCGTCTAGGTCAAACGTTATGGTAAGTGGCGAATGATGGTTATAGCAGTGGCGGCAGGAG‑TACCGTTTTACA‑CCGG‑ATTGAGTGTATTTCATAGCCGTCTGGCG‑TTGATGACGCTTCAACTTTTTAACGCTGTATGTATCGGCATTGTTGCGGGGATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCtt                                                                                                                     <  2:575355/236‑2 (MQ=255)
                                                                                                                gggACGGCAGCTGGTCT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGACTGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGATGCGGGTATTGGGTGGCTTTGGTTTCAGGTGTTAATGCCTggtaggaggggggaagctac                                                                                                                    >  2:474392/1‑233 (MQ=255)
                                                                                                                gggACGGCAGCTGGGCT‑GGAAATACCAGCAATGATGCTGGCTGGCTACTATGTCAAACGTTATTGTAAGCGGCGAAGGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTac                                                                                                                    <  2:117932/251‑1 (MQ=255)
                                                                                                                     ggCAGCTGGTAT‑GGAAATACCAGCAATGATTCTGGCTGGCTAATATGTCAACCGTTATGGTAAGCGGCGAATGATGGTAATAGCCGTGGCGGCAGGGG‑TGCTGTTTTCCA‑CAGG‑ATTGCTTTTGTTTCATAGCAGTCTGGCG‑TTGGTGACGCAGCACCTTTTTAATGCTTTATTTATTGGCATTGGTGAGGGTATTGGGCTGATATGGTTTCATGGTTTAATGACTGGAAGAGCGGGagcagcaaaaaact                                                                                                               >  2:874838/1‑243 (MQ=255)
                                                                                                                        aGCTGGTCT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGGCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATTGTTTCAGGATTTAATGCCTGGAAGGGCGGGGGCAGCTACCCCCTTAt                                                                                                            >  2:304552/1‑251 (MQ=255)
                                                                                                                          ctggtct‑ggaAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATGGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCgg                                                                                                                              <  1:316024/231‑1 (MQ=255)
                                                                                                                          ctggtct‑ggaAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGAAGGAG‑TACTATTTTACC‑CCGG‑CTTGATTTTATTTCATAGCCGTATGGCG‑TTGATTACGCTGCTAATTTTTAACGCTGTATTTATCGGCATTGTTGGGGGTATTGGGATGATTTGGATTTAGGTATTAATGCATGGAAGAGAgggagtgagaagaggggagtgt                                                                                                          >  2:316024/1‑232 (MQ=255)
                                                                                                                           tggtct‑ggaAATACCAGCAATGATTCTAGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑AACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTAGCG‑TTGATGACGCTGCAACTTTTTTACGCAGTATTTATAGGCATTGTTGCGGGTGTTGGGATGCTATGGTTTCAGGATTTAATACCTGGAAGAGCGGGGGCAGcataccccttctttt                                                                                                         >  2:80177/5‑237 (MQ=255)
                                                                                                                              tct‑ggaAATACCAGCAATGATTCTGGCTGGCTACTATGTCAACCGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACAGTTTTACA‑CCGG‑GTTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTTTTTATCGGCCATGTTGCGGGGATTGGGAAGCTATGGGTTCAGGGTTTAAATCCTGGAAGGGCGGGGGGAGCTacaccctttattact                                                                                                       >  2:911051/2‑237 (MQ=255)
                                                                                                                                  ggaAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAAc                                                                                                    >  1:265185/1‑250 (MQ=255)
                                                                                                                                      aTACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAAAGCTGTATTTATCGGCATTGTTGCGGGTATTGGGTTGTTATTGTTTCAGGATTTATTGCCAGGAAGAGCGGGGGCAGGTACCACCTTATTTACTAACAGttt                                                                                               >  2:572793/1‑249 (MQ=255)
                                                                                                                                          cAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGGTGCTATGGTTTCAGGATTTAATGCCTGGCAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCt                                                                                           >  1:210756/1‑251 (MQ=255)
                                                                                                                                          cAGCAATGATTCGGGCTGGCTACTATGTCAAACGGTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCt                                                                                           <  1:298967/251‑1 (MQ=255)
                                                                                                                                           aGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTa                                                                                          <  1:758944/251‑1 (MQ=255)
                                                                                                                                             cAATGAGTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTAcc                                                                                        <  2:898703/251‑1 (MQ=255)
                                                                                                                                              aaTGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCg                                                                                       <  2:400660/251‑1 (MQ=255)
                                                                                                                                                tGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCACCTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGCTGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAACTACCACCTTATTTACTAACAGTATTTCTACCggg                                                                                     >  1:347714/1‑251 (MQ=255)
                                                                                                                                                  aTTCTGGCTGGCTACTCTGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑GACTGGTTTACA‑CCGG‑GTTGATTTTATTTTATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTTTTTATCGGCATTGTTGCGGGTATTGGCGTGCTTTGGTTTCAGGATTTAATTCCTGGAAGAGCGGGGGCAGGTACCACCTTATTTACTAAAAGTATTTATACCggggg                                                                                   >  2:188578/1‑250 (MQ=255)
                                                                                                                                                  aTTCTGGCTGGCTACTATGTCAAACGTTATGGAAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCCTTGTTGCGGGAATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTAACACCTTTTTTACTAACAGTATTTCTACCGGGGt                                                                                   >  1:949142/1‑251 (MQ=255)
                                                                                                                                                       ggctggctACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTaa                                                                                 <  1:109688/248‑1 (MQ=255)
                                                                                                                                                       ggctggctACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTc                                                                              >  1:787733/1‑251 (MQ=255)
                                                                                                                                                                         aCGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCTGGGGAATTTTGGCTGGCGTTATTCAgg                                                            >  2:659957/1‑251 (MQ=255)
                                                                                                                                                                              aTGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCa                                                       <  2:117086/251‑1 (MQ=255)
                                                                                                                                                                               tGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCCTTGTTGCGGGTATTGGGATGCTATGGTTTCAGGTTTTCATGCCTGGAAGAGCGGGGGCCGCTACCTCCTTATTTTCTAACAGTATTCCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCa                                                       >  2:348915/1‑250 (MQ=255)
                                                                                                                                                                                     gcggcgAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTTTTTATCGGCCTTTTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAAAAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTccgg                                                            >  2:891862/1‑236 (MQ=255)
                                                                                                                                                                                     gcggcgAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAgg                                                            <  1:891862/239‑1 (MQ=255)
                                                                                                                                                                                      cggcgAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGcac                                                >  1:656708/1‑250 (MQ=255)
                                                                                                                                                                                       ggcgAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTTCA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACaa                                              <  1:71345/251‑1 (MQ=255)
                                                                                                                                                                                            atgatgGTCATAGCAGGTGCGGCATGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTGTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGTGCAGCGTCCACCTTATTTACTAACAGTATTTCTACCTGGGTAATTCTGGCTGGCGGTATTCAGGGAGCAATGGCACAAATTTg                                         <  2:626798/251‑1 (MQ=255)
                                                                                                                                                                                                   tCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCAc                                  >  1:964614/1‑251 (MQ=255)
                                                                                                                                                                                                       aGCAGTGGCGGTAGGAG‑TACTGTGTTCCA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTAGCGGCAGTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACttt                               <  2:617770/250‑1 (MQ=255)
                                                                                                                                                                                                       aGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACttt                               <  1:523889/250‑1 (MQ=255)
                                                                                                                                                                                                           gTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGt                          <  1:587185/251‑1 (MQ=255)
                                                                                                                                                                                                             ggcggcAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTc                         <  1:411814/250‑1 (MQ=255)
                                                                                                                                                                                                                             gTTTTACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAAGTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCg          <  1:401280/250‑1 (MQ=255)
                                                                                                                                                                                                                                ttACA‑CCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGtt       <  1:304552/250‑1 (MQ=255)
                                                                                                                                                                                                                                  aca‑cCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTAt     <  1:230120/250‑1 (MQ=255)
                                                                                                                                                                                                                                  aca‑cCGG‑AATGATTTTATTTCATAACCGGCTGGCG‑TTGATGACGCTTCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACAGGGGTAATTCTGGCTGGCGTTATTCATGGAGCAATTGCACACAGTTGGGGGCACTTTGCTGTATACTGGGTAATTGCGGTTAtt    >  2:716131/1‑251 (MQ=255)
                                                                                                                                                                                                                                   ca‑cCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCATTTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTAttt   >  2:370498/1‑251 (MQ=255)
                                                                                                                                                                                                                                    a‑cCGG‑ATTGATTTTATTTCATAGCCGTCTGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTGTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTTCCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGTGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGTGTAATTGCGGTTATTTc  <  2:538369/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                             
ATTTGTCGCCTCGACGTTAATGTGGACCTGCAACACCATGTACATTATTGATATGCCGTTGTGGATCAGTAGCGAGTTAGGATTGCCAGACAAACTGGCGGGTTTCCTGATGGGGACGGCAGCTGGACT‑GGAAATACCAGCAATGATTCTGGCTGGCTACTATGTCAAACGTTATGGTAAGCGGCGAATGATGGTCATAGCAGTGGCGGCAGGAG‑TACTGTTTTACA‑CCGG‑ATTGATTTTCTTTAATAGCCGTATGGCG‑TTGATGACGCTGCAACTTTTTAACGCTGTATTTATCGGCATTGTTGCGGGTATTGGGATGCTATGGTTTCAGGATTTAATGCCTGGAAGAGCGGGGGCAGCTACCACCTTATTTACTAACAGTATTTCTACCGGGGTAATTCTGGCTGGCGTTATTCAGGGAGCAATTGCACAAAGTTGGGGGCACTTTGCTGTCTACTGGGTAATTGCGGTTATTTC  >  NC_000913/78280‑78756

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