Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,070,469 A→T 7.2% intergenic (‑68/+46) yihR ← / ← yihS putative sulphoquinovose mutarotase/sulphoquinovose isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,070,4690AT 7.2% 1604.5 / 13.5 525intergenic (‑68/+46)yihR/yihSputative sulphoquinovose mutarotase/sulphoquinovose isomerase
Reads supporting (aligned to +/- strand):  ref base A (391/67);  new base T (2/47);  total (402/123)
Fisher's exact test for biased strand distribution p-value = 3.98e-32
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 571 total aligned reads displayed.

AATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATAT  >  NC_000913/4070320‑4070618
                                                                                                                                                         |                                                                                                                                                     
aaTCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGg‑‑a                                                                                                                                                       <  1:3439004/150‑1 (MQ=255)
aaTCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGg‑‑a                                                                                                                                                       >  1:2827505/1‑150 (MQ=255)
 aTCCCCCGCGGCAAGTGAAACGGTCTGTCCCCTGCAATGCATATTGGTGAGTTGCAGACAAGGCACCATAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAACACTGCCCCGGACGGTCCCCTCGCCCC‑GGGGGGg‑‑ag                                                                                                                                                      >  2:6719523/1‑150 (MQ=255)
 aTCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGg‑‑ag                                                                                                                                                      <  2:10005380/150‑1 (MQ=255)
 aTCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGg‑‑ag                                                                                                                                                      <  1:4738000/150‑1 (MQ=255)
  tCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGg‑‑aga                                                                                                                                                     <  2:6553082/150‑1 (MQ=255)
   cccccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑agag                                                                                                                                                    <  2:1899069/150‑1 (MQ=255)
   cccccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑agag                                                                                                                                                    <  1:7496343/150‑1 (MQ=255)
     cccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCTCCCC‑AAAGGGG‑‑CGAggg                                                                                                                                                  <  1:9422170/150‑1 (MQ=255)
      ccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGt                                                                                                                                                 >  2:2837996/1‑150 (MQ=255)
      ccGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGt                                                                                                                                                 >  2:3038130/1‑150 (MQ=255)
        gcgGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAACACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTa                                                                                                                                               <  2:961277/150‑1 (MQ=255)
         cgGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTACCCTCTCCCC‑AAAGGGG‑‑CGAggggacc                                                                                                                                              <  1:6333476/150‑5 (MQ=255)
           gCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAggg                                                                                                                                            >  1:3923284/1‑150 (MQ=255)
            cAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGCCATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGga                                                                                                                                           <  2:9126204/150‑1 (MQ=255)
             aaGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TGTGGGG‑‑AGAGGGTTAGGgag                                                                                                                                          >  1:2005691/1‑150 (MQ=255)
              aGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGgaga                                                                                                                                         >  2:7056358/1‑150 (MQ=255)
               gTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGagag                                                                                                                                        <  2:3497421/150‑1 (MQ=255)
               gTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGagag                                                                                                                                        >  1:8058342/1‑150 (MQ=255)
                tGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAATCCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGTgg                                                                                                                                       >  2:9539336/1‑150 (MQ=255)
                tGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAgg                                                                                                                                       <  1:2696898/150‑1 (MQ=255)
                tGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TGTGGGG‑‑AGAGGGTTAGGGAGAgg                                                                                                                                       >  2:1029578/1‑150 (MQ=255)
                 gAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAggg                                                                                                                                      >  1:290269/1‑150 (MQ=255)
                  aaaCGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAgggg                                                                                                                                     >  2:986161/1‑150 (MQ=255)
                  aaaCGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAgggg                                                                                                                                     >  2:5656919/1‑150 (MQ=255)
                  aaaCGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAgggg                                                                                                                                     >  1:10810989/1‑150 (MQ=255)
                   aagggTCTGTCCACTGAAATGAATTTTGGTGTTTTGCATAAACGGCACATTATTTAACGACATACAGCCTCCCTTGTGGCTG‑ATGTAAAACCAAAAAACTGCACGGGACGGTCCCCTCGCCCC‑TTTGGGT‑‑TGAGGGTTCGGTTGAGGGGa                                                                                                                                    <  2:3205217/147‑1 (MQ=255)
                    ccGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTACTTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑GGTGGGG‑‑AGAGGGTTAGGGAGAGGGGaa                                                                                                                                   >  2:445854/2‑150 (MQ=255)
                       gtgtGTCCAGTGAAATGTATTTGGTTTTTTTGCATACATGGCACCTTAATTAACGACATACCGCCCCCCTTGAGGCTG‑ATGTAAAACCCAAACACTGCCCCGGAGGGTCCCCTCGCCCC‑TTTGGGT‑‑TGTGGGTTAGGTTGTGGGGAAAAc                                                                                                                                <  2:404140/147‑1 (MQ=255)
                          tGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑GGTGGGG‑‑AGAGGGGTAGGGAGAGGGGAAAACCCg                                                                                                                             >  2:4971300/1‑150 (MQ=255)
                            tCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCa                                                                                                                           >  2:5292004/1‑150 (MQ=255)
                            tCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑GTTGGGG‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCa                                                                                                                           >  2:8203500/1‑150 (MQ=255)
                              cACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCAcc                                                                                                                         >  1:5706467/1‑150 (MQ=255)
                                cTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCt                                                                                                                       >  1:10524786/1‑150 (MQ=255)
                                     aTGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAtttt                                                                                                                  >  2:5497940/1‑150 (MQ=255)
                                        cATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAg                                                                                                               >  2:9498108/1‑150 (MQ=255)
                                         attttGGTTATTTGCATAAAAGGCACATTAATAAACGACATACAGCCCCCCTTGTGGATG‑ATGTAAAACCCAAAAACTCCCCCGGACGGGCCCCTCGCCCC‑TTTGGGG‑‑AGTGGTTTCGGGAGTGGGGAAAACCCGCACCCTATTTTCAGa                                                                                                              <  2:7378760/147‑1 (MQ=255)
                                          tatTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGT‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGAt                                                                                                             <  1:11090886/150‑1 (MQ=255)
                                             tGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAGGGGACAACCCGCACCCTATTTTCAGATCAg                                                                                                          >  1:8848833/1‑150 (MQ=255)
                                              ggTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGt                                                                                                         >  1:1095600/1‑150 (MQ=255)
                                              ggTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGt                                                                                                         >  1:3612520/1‑150 (MQ=255)
                                               gTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGtt                                                                                                        >  1:10069678/1‑150 (MQ=255)
                                                tGATTTGCATACAAGGCAGCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑TTGTCAAACCCAAAAACTGCACCGGTCGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAGGGGTAAACCCGCACCCTATTTTCAGTTCAGTTa                                                                                                       >  2:5430807/1‑150 (MQ=39)
                                                 gATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTAt                                                                                                      >  1:2400023/1‑150 (MQ=255)
                                                     tcaaTACACGCCACATTAATTAACCACATACAGCCCCCCGTGCGGCTG‑ATTAAAAACCCAAAAACGGCACGGGACGGGCCCCTCCCCCC‑TTTGGGG‑‑AGTGGGTTTGGGTGCGGGCAAAACCCGCACCCTATTTTCAGATCAGTTATTTCg                                                                                                  <  2:10928990/147‑1 (MQ=255)
                                                       cATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCa                                                                                                >  2:835516/1‑150 (MQ=255)
                                                        aTACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCAt                                                                                               >  1:5442196/1‑150 (MQ=255)
                                                          aCAAGCCACTTAAATAAACCACAAACACCCTCCCGTGTGGCTT‑TTTAAAAACCAAAACACGCCACCGGAGGGCCCCCCCCCCCC‑TTTGGGG‑‑AGCGGTTTCGGGTGTGGGTAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTa                                                                                             <  1:9112405/150‑1 (MQ=255)
                                                                cACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGCCGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGAGCAGTTATTTCGCATTAATATcc                                                                                       >  1:1539393/1‑150 (MQ=255)
                                                                     taattaaCGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCcagcag                                                                                  >  2:1021409/1‑150 (MQ=255)
                                                                      aattaaCGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGa                                                                                 >  1:1373496/1‑150 (MQ=255)
                                                                       attaaCGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGAc                                                                                >  2:6353609/1‑150 (MQ=255)
                                                                        ttaaCGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCCATTTCCAGATCAGTTATTTCGCATTAATATCCAGCAGAcc                                                                               >  2:8702456/1‑150 (MQ=255)
                                                                         taaCGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGAccc                                                                              >  1:4283649/1‑150 (MQ=255)
                                                                            cGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCt                                                                           >  2:1643269/1‑150 (MQ=255)
                                                                             gACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTg                                                                          >  2:5847260/1‑150 (MQ=255)
                                                                               acatcAGCCTCACGGGAGGCTGTATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGc                                                                         >  1:9247976/5‑150 (MQ=39)
                                                                                aTACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCaa                                                                       >  2:6919159/1‑150 (MQ=255)
                                                                                 tACCGCCTCCCGTGTGGCTT‑ATGACAAACCAAAAAACGCCACGGGAGGGCCCCCTCGCCCC‑TTGGGGG‑‑TGGGGTTTTGGGAGTGGGCAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAAc                                                                      <  1:8188047/150‑1 (MQ=255)
                                                                                 tACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAAc                                                                      >  2:7736532/1‑150 (MQ=255)
                                                                                 tACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAAc                                                                      >  1:8697176/1‑150 (MQ=255)
                                                                                     gCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCt                                                                  >  1:4415684/1‑150 (MQ=255)
                                                                                      ccTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTGATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTg                                                                 >  2:2172398/1‑150 (MQ=255)
                                                                                          ccGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGc                                                             >  2:5090246/1‑150 (MQ=255)
                                                                                          ccGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGc                                                             >  2:6243919/1‑150 (MQ=255)
                                                                                           cGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTCGGGAGAGGGGAAAACCCGCCCCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGcc                                                            >  1:590538/1‑150 (MQ=255)
                                                                                             tgcggCTT‑ATGACAAACCAAAAACCCCCCCCGGCGGGCCCCCGCGCCCC‑TTGGGGG‑‑TGTGGTTTTGGGAGCGGGCAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCAt                                                          <  1:6879578/147‑1 (MQ=255)
                                                                                                         tCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTa                                               >  1:3241946/1‑150 (MQ=255)
                                                                                                         tCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTa                                               >  1:9697319/1‑150 (MQ=255)
                                                                                                           aaaCCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAAc                                             >  2:1461578/1‑150 (MQ=255)
                                                                                                               ccAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGa                                         >  1:2199859/1‑150 (MQ=255)
                                                                                                                cAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGAt                                        >  1:2990699/1‑150 (MQ=255)
                                                                                                                cAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TATGGGG‑‑AGAGGGTAAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGAACAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGAt                                        >  2:6603948/1‑150 (MQ=255)
                                                                                                                 aaaaaCTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAGCGGGATa                                       >  1:7049198/1‑150 (MQ=255)
                                                                                                                 aaaaaCTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGAtc                                       >  2:4636585/1‑149 (MQ=255)
                                                                                                                 aaaaaCTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATa                                       >  1:4940174/1‑150 (MQ=255)
                                                                                                                 aaaaaCTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATa                                       >  1:1039058/1‑150 (MQ=255)
                                                                                                                 aaaaaCTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATa                                       >  1:6077031/1‑150 (MQ=255)
                                                                                                                 aaaaaCTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATa                                       >  1:6666483/1‑150 (MQ=255)
                                                                                                                  aaaaCTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATAc                                      >  1:10254358/1‑150 (MQ=255)
                                                                                                                       tGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCggg                                 >  2:1536863/1‑150 (MQ=255)
                                                                                                                       ccccccGGGCGGGCCCCCCCCCCCGTTGGGGG‑‑TGTGGTTGTGGGTGGGGGCAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCgg                                  <  1:9847756/146‑1 (MQ=255)
                                                                                                                            cGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCAc                            >  2:5988716/1‑150 (MQ=255)
                                                                                                                             ggccGGTCCCCTCGCCCT‑TTTGGGT‑‑TGTGGGTTTGGTTGTGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCAcc                           <  2:2575005/147‑1 (MQ=255)
                                                                                                                              cggcggccccccGCCCC‑CTTGGGGGTTGTGGGTTGGGGCGGGGCCAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCAc                            <  1:7833990/144‑1 (MQ=255)
                                                                                                                                  taacccTCTCCCC‑AAAGGGG‑‑CGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCa                      >  2:7556103/4‑150 (MQ=255)
                                                                                                                                   tCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCaa                     >  1:737646/1‑150 (MQ=255)
                                                                                                                                     cccTCTCCCC‑AAAGGGG‑‑CGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATg                   >  2:10507725/1‑150 (MQ=255)
                                                                                                                                     cccTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATg                   >  1:8466541/1‑150 (MQ=255)
                                                                                                                                       cccGCCCC‑TGGGGGT‑‑TGTGGGTTTGGGTGCGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATgca                 <  2:3751885/148‑1 (MQ=255)
                                                                                                                                       cTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATgca                 >  2:5215951/1‑150 (MQ=255)
                                                                                                                                       cTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATgca                 >  2:4257894/1‑150 (MQ=255)
                                                                                                                                         cGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATgcagc               >  2:4897889/1‑150 (MQ=255)
                                                                                                                                          gCCCC‑TTGGGGG‑‑AGTGGGTTTGGGAGCGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATgcagca              <  2:10755894/150‑1 (MQ=255)
                                                                                                                                                tttGGGG‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGa         >  2:11080581/1‑150 (MQ=255)
                                                                                                                                                ttggggg‑‑TGTGGTTTGGGGTGGGGGCAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGa         <  1:6444829/147‑1 (MQ=255)
                                                                                                                                                     gg‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAat    >  2:8087169/1‑150 (MQ=255)
                                                                                                                                                     gg‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAat    >  2:1148855/1‑150 (MQ=255)
                                                                                                                                                      g‑‑agaGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAata   >  2:7904477/1‑150 (MQ=255)
                                                                                                                                                         agagGGTTAGGGAGAGGGGAAAACCCGCTCCCTATTTTGAGATCAGTTATTTCGCATTAATATCCAGCAGTCGCGCTGCTACCGGTGGAGCCATGGCGGGGGCTATCGGGATACGCGGGATCTCCAGGCAATGGAGCTGGTGATAAttat  >  2:2037248/1‑150 (MQ=255)
                                                                                                                                                         |                                                                                                                                                     
AATCCCCCGCGGCAAGTGAAACGGTCTGTCCACTGCAATGCATATTGGTGATTTGCATACAAGGCACCTTAATTAACGACATACAGCCTCCCGTGAGGCTG‑ATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCC‑TTTGGGG‑‑AGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTTCAGATCAGTTATTTCGCATTAATATCCAGCAGACCCGCTGCAACCGCTGGAGCCATGCCGGGCGCTAACGGGATACGCGGGATCACCAGGCAATGCAGCAGGTGATAAATAT  >  NC_000913/4070320‑4070618

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑