Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 191,730 A→C 5.0% intergenic (+22/‑125) tsf → / → pyrH translation elongation factor EF‑Ts/uridylate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913191,7300AC 5.0% 354.9 / 8.6 100intergenic (+22/‑125)tsf/pyrHtranslation elongation factor EF‑Ts/uridylate kinase
Reads supporting (aligned to +/- strand):  ref base A (52/43);  new base C (3/2);  total (55/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

Only 100 of 101 total aligned reads displayed.

CTGAAAGAGCATAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTC  >  NC_000913/191592‑191878
                                                                                                                                           |                                                                                                                                                    
ctgAAAGAGCATAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGcttc                                                                                                                                            >  2:161738/1‑149 (MQ=255)
    aaGAGCATAACGCTGAAGTGACTGGCTTCATCCTCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTAAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCtttt                                                                                                                                        <  4:583460/149‑1 (MQ=255)
    aaGAGCATAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCCTCAGGCGGCTCCtttt                                                                                                                                        <  1:61602/149‑1 (MQ=255)
     agagCATAACGCTGAAGTGACTGGCTTCATCCGCTTCGACGTGTGTGAAGGCATCGAGAAAGCTGAGACTGACTTTACGGCAGAAGTTGCTGCGATGTCCAAGCAGTATTAATTATCAAAAAGGAGCCG‑CCTGATGGCGGCTTCttttt                                                                                                                                       >  3:544899/1‑149 (MQ=255)
     agagCATAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCttttt                                                                                                                                       <  1:428399/149‑1 (MQ=255)
     agagCATAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGAAGCCGCCCTCA‑GGCGGCTCCttttt                                                                                                                                       <  2:454520/149‑1 (MQ=255)
     agagCATAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGAAGCCG‑CCCTCAGGCGGCTCCttttt                                                                                                                                       <  1:160011/149‑1 (MQ=255)
          aTAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGcc                                                                                                                                  >  4:414458/1‑149 (MQ=255)
           tAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGccc                                                                                                                                 >  2:216003/1‑149 (MQ=255)
            aaCGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCa                                                                                                                                >  1:772135/1‑149 (MQ=255)
                 tGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTg                                                                                                                           >  2:398599/1‑149 (MQ=255)
                   aaGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTa                                                                                                                         >  4:536953/1‑149 (MQ=255)
                         cTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCa                                                                                                                   >  3:195677/1‑149 (MQ=255)
                            gCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCt                                                                                                                >  4:476967/1‑149 (MQ=255)
                             cTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTa                                                                                                               >  1:156203/1‑149 (MQ=255)
                              ttCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTaa                                                                                                              >  3:58235/1‑149 (MQ=255)
                              ttCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTaa                                                                                                              <  1:170603/149‑1 (MQ=255)
                              ttCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTaa                                                                                                              >  4:754780/1‑149 (MQ=255)
                              ttCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTaa                                                                                                              >  3:872301/1‑149 (MQ=255)
                                           aaGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCt                                                                                                 >  1:295816/1‑149 (MQ=255)
                                           aaGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCt                                                                                                 >  3:593238/1‑149 (MQ=255)
                                               gggTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCtttt                                                                                                                                        >  1:210515/1‑106 (MQ=255)
                                               gggTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCtttt                                                                                                                                        <  2:210515/106‑1 (MQ=255)
                                               gggTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTgcgc                                                                                             >  1:794530/1‑149 (MQ=255)
                                                 gTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTg                                                                                           >  2:252756/1‑149 (MQ=255)
                                                 gTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTg                                                                                           >  2:692922/1‑149 (MQ=255)
                                                   gAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGaa                                                                                         >  3:817197/1‑149 (MQ=255)
                                                     aGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGaaaa                                                                                       <  4:410170/149‑1 (MQ=255)
                                                     aGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGaaaa                                                                                       <  2:658915/149‑1 (MQ=255)
                                                      ggCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAg                                                                                      <  2:101258/149‑1 (MQ=255)
                                                       gCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCttttt                                                                                                                                       <  4:268222/99‑1 (MQ=255)
                                                       gCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCttttt                                                                                                                                       >  3:268222/1‑99 (MQ=255)
                                                       gCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGc                                                                                     <  1:625304/149‑1 (MQ=255)
                                                       gCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGc                                                                                     <  2:437228/149‑1 (MQ=255)
                                                         aTCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGa                                                                                   >  2:506313/1‑149 (MQ=255)
                                                         aTCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGa                                                                                   <  3:754780/149‑1 (MQ=255)
                                                         aTCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGa                                                                                   >  1:126489/1‑149 (MQ=255)
                                                         aTCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGa                                                                                   >  2:310101/1‑149 (MQ=255)
                                                          tCGAGAAAGTTGAGACTGACTTTGCAGCAGTAGTAGCTGCGATGTCCAAGCAGTCTCAATTATCAAAAAGGAGCCG‑CCTCAGGGCGGCTTCTTTTTGTGACCATCTTATAAATTCAGCTAACCCTTGTGGGACTGCGCTGAAAAGCGAc                                                                                  <  4:544899/149‑1 (MQ=255)
                                                          tCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGAc                                                                                  <  1:692922/149‑1 (MQ=255)
                                                            gagaAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGt                                                                                >  3:591619/1‑149 (MQ=255)
                                                                aaGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACaa                                                                            <  3:376058/149‑1 (MQ=255)
                                                                 aGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAAt                                                                           <  1:256128/149‑1 (MQ=255)
                                                                   ttGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGt                                                                         >  4:301725/1‑149 (MQ=255)
                                                                      agaCTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGc                                                                      <  2:206666/149‑1 (MQ=255)
                                                                       gactgactTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCt                                                                     >  3:93886/1‑149 (MQ=255)
                                                                       gactgactTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCt                                                                     >  2:441543/1‑149 (MQ=255)
                                                                        actgactTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTa                                                                    >  2:586602/1‑149 (MQ=255)
                                                                        actgactTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTa                                                                    >  4:282807/1‑149 (MQ=255)
                                                                          tgactTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGt                                                                  <  1:441543/149‑1 (MQ=255)
                                                                          tgactTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGt                                                                  >  3:235803/1‑149 (MQ=255)
                                                                               ttGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTaa                                                             <  1:323047/149‑1 (MQ=255)
                                                                                      agAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCAttt                                                      <  2:167568/149‑1 (MQ=255)
                                                                                         aGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCaa                                                   <  1:216003/149‑1 (MQ=255)
                                                                                          gTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAAt                                                  <  3:627535/149‑1 (MQ=255)
                                                                                              ctgcGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGtt                                              <  1:310101/149‑1 (MQ=255)
                                                                                              ctgcGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGtt                                              <  3:536953/149‑1 (MQ=255)
                                                                                              ctgcGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGtt                                              <  3:476967/149‑1 (MQ=255)
                                                                                               tgcGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAAt                                                                           >  4:722614/1‑119 (MQ=255)
                                                                                               tgcGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAAt                                                                           <  3:722614/119‑1 (MQ=255)
                                                                                               tgcGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTg                                             >  4:588166/1‑149 (MQ=255)
                                                                                               tgcGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTg                                             >  1:668278/1‑149 (MQ=255)
                                                                                               tgcGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGCGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTg                                             >  2:722476/1‑149 (MQ=255)
                                                                                                gcGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGa                                            <  4:73126/149‑1 (MQ=255)
                                                                                                 cGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGAc                                           >  3:149928/1‑149 (MQ=255)
                                                                                                  gATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACa                                          <  4:45521/149‑1 (MQ=255)
                                                                                                      tCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGtct                                      >  3:743415/1‑149 (MQ=255)
                                                                                                       ccAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGtctc                                     <  2:460171/149‑1 (MQ=255)
                                                                                                        cAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCa                                    <  3:701459/149‑1 (MQ=255)
                                                                                                         aaGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAg                                   >  3:513809/1‑149 (MQ=255)
                                                                                                          aGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAgg                                  >  3:652353/1‑149 (MQ=255)
                                                                                                             aGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGgaaa                               <  3:106746/149‑1 (MQ=255)
                                                                                                               tCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGgaaaga                             >  1:247387/1‑149 (MQ=255)
                                                                                                               tCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGgaaaga                             <  1:386682/149‑1 (MQ=255)
                                                                                                                 ttAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGgaaagaaa                           <  3:679794/149‑1 (MQ=255)
                                                                                                                 ttAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGgaaagaaa                           >  1:271845/1‑149 (MQ=255)
                                                                                                                 ttAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGgaaagaaa                           >  3:10670/1‑149 (MQ=255)
                                                                                                                 ttAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGgaaagaaa                           >  1:570908/1‑149 (MQ=255)
                                                                                                                 ttAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGgaaagaaa                           <  4:711863/149‑1 (MQ=255)
                                                                                                                 ttAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGgaaagaaa                           <  3:790990/149‑1 (MQ=255)
                                                                                                                 ttAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGgaaagaaa                           <  3:426325/149‑1 (MQ=255)
                                                                                                                   aaTTATCAAAAAGGAGTCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATACATTACAATCGTTGACATTCTCAGGCAAGAcaca                         >  2:435526/1‑149 (MQ=255)
                                                                                                                     ttATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATg                       <  3:744069/149‑1 (MQ=255)
                                                                                                                     ttATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATg                       >  1:94469/1‑149 (MQ=255)
                                                                                                                        tCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCt                    <  4:284961/149‑1 (MQ=255)
                                                                                                                          aaaaaGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTAc                  >  1:12818/1‑149 (MQ=255)
                                                                                                                          aaaaaGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTAc                  >  2:521352/1‑149 (MQ=255)
                                                                                                                          aaaaaGAAGCCGCCCTCA‑GGCGGCTCCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTAc                  >  2:409073/1‑149 (MQ=255)
                                                                                                                          aaaaaGAAGCCG‑CCCTCAGGCGGCTCCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTAc                  >  2:251380/1‑149 (MQ=255)
                                                                                                                           aaaaGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATc                                                                                                                              >  4:24521/1‑40 (MQ=255)
                                                                                                                           aaaaGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATc                                                                                                                              <  3:24521/40‑1 (MQ=255)
                                                                                                                           aaaaGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTAACAGTCTCa                                    <  2:739486/130‑1 (MQ=255)
                                                                                                                           aaaaGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTAACAGTCTCa                                    >  1:739486/1‑130 (MQ=255)
                                                                                                                           aaaaGAAGCCG‑CCCTCAGGCGGCTCCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTAcc                 >  4:465842/1‑149 (MQ=255)
                                                                                                                            aaaGAAGCCG‑CCCTCAGGCGGCTCCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCa                >  2:377450/1‑149 (MQ=255)
                                                                                                                              aGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAAt              <  2:94469/149‑1 (MQ=255)
                                                                                                                                  gccg‑ccTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCaa          >  2:745614/1‑149 (MQ=255)
                                                                                                                                     g‑ccTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAAc       >  2:396354/1‑149 (MQ=255)
                                                                                                                                          gAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGt   <  2:570908/149‑1 (MQ=255)
                                                                                                                                           aGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTc  <  2:794530/149‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                    
CTGAAAGAGCATAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTC  >  NC_000913/191592‑191878

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