Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,538,556 T→G 6.4% intergenic (+73/+116) yfeS → / ← cysM WGR domain protein/cysteine synthase B (O‑acetylserine sulfhydrolase B)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,538,5560TG 6.4% 359.7 / 9.4 119intergenic (+73/+116)yfeS/cysMWGR domain protein/cysteine synthase B (O‑acetylserine sulfhydrolase B)
Reads supporting (aligned to +/- strand):  ref base T (49/62);  new base G (4/4);  total (53/66)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 120 total aligned reads displayed.

ATGGGGTGACGGTCAGGACTCACCAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTT  >  NC_000913/2538408‑2538703
                                                                                                                                                    |                                                                                                                                                   
  ggggTGACGGTCAGGACTCACCAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTAt                                                                                                                                                  >  1:45046/1‑150 (MQ=255)
      tGACGGTCAGGACTCACCAATTCTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCTCGAATATTAAATAGTGCATTAGTCGGCAGAAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTTACGCCGCATCCgg                                                                                                                                              <  2:1059594/150‑1 (MQ=38)
        aCGGTCAGGACTCACCAATACTTCAAAAAATGCTTGATGACCTTACGGCGTTTCCTCACGAAACTTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                                                                                            <  1:813791/150‑1 (MQ=255)
        aCGGTCAGGACTCACCAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCg                                                                                                                                               >  1:265967/1‑147 (MQ=255)
        aCGGTCAGGACTCACCAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCg                                                                                                                                               <  2:265967/147‑1 (MQ=255)
        aCGGTCAGGACTCACCAATACATCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                                                                                            <  2:1193291/150‑1 (MQ=255)
          ggTCAGGACTCACCAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa                                                                                                                                          >  2:1099988/1‑150 (MQ=255)
          ggTCAGGACTCACCAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa                                                                                                                                          >  2:156760/1‑150 (MQ=255)
          ggTCAGGACTCACCAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGAGGCGTAAACGCCTTATCCGGCCTa                                                                                                                                          >  2:1328208/1‑150 (MQ=255)
               ggACTCACCAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaa                                                                                                                                     >  2:1502171/1‑150 (MQ=255)
                gACTCACCAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAt                                                                                                                                    <  1:572643/150‑1 (MQ=255)
                   tCACCAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGt                                                                                                                                 <  1:1473097/150‑1 (MQ=255)
                    cACCAATACTTCAAAAAATGATTGCTGACCTTACGGCGTTTCCTCACGAAAATTAATTACAGCATTTGTCGGCAGCAACACCTGTTAAAAATGGGCCCTTTGTATATTCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTg                                                                                                                                >  1:78097/1‑150 (MQ=37)
                     aCCAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATTCAGCATTTGTCGGCAGCATCAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGc                                                                                                                               <  2:182535/150‑1 (MQ=255)
                     aCCAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTACGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGc                                                                                                                               <  1:788208/150‑1 (MQ=255)
                      ccAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCt                                                                                                                              >  1:356954/1‑150 (MQ=255)
                      ccAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCt                                                                                                                              <  2:603778/150‑1 (MQ=255)
                       cAATACTTCAAAAACTGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTa                                                                                                                             <  1:282179/150‑1 (MQ=255)
                           aCTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGCTGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTc                                                                                                                         <  1:283099/150‑1 (MQ=255)
                           aCTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTc                                                                                                                         <  1:279264/150‑1 (MQ=255)
                             ttCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGAAGCGGCGTAATCGCCTTATCCGGCCTACATAATCGTGCTAATTCaa                                                                                                                       >  1:88189/1‑150 (MQ=255)
                              tCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAat                                                                                                                      <  2:511757/150‑1 (MQ=255)
                              tCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAat                                                                                                                      <  2:70521/150‑1 (MQ=255)
                               cAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAata                                                                                                                     <  1:1621258/150‑1 (MQ=255)
                                    aaTGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTg                                                                                                                <  2:1529203/150‑1 (MQ=255)
                                     aTGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGc                                                                                                               <  2:616815/150‑1 (MQ=255)
                                     aTGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGc                                                                                                               >  1:393165/1‑150 (MQ=255)
                                           gatgaCCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAAc                                                                                                         >  1:156061/1‑150 (MQ=255)
                                           gatgaCCTTACGGCGTGTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGAAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAAc                                                                                                         <  2:1220820/150‑1 (MQ=255)
                                            atgaCCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAAcc                                                                                                        <  1:801435/150‑1 (MQ=255)
                                             tgaCCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCt                                                                                                       >  1:397647/1‑150 (MQ=255)
                                                ccTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGt                                                                                                    >  1:720115/1‑150 (MQ=255)
                                                ccTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGt                                                                                                    <  2:1306575/150‑1 (MQ=255)
                                                 cTTACGGCGTTTCCTCACGAAAATTAAATACAGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTa                                                                                                   <  2:1441027/150‑1 (MQ=255)
                                                       gCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGCAACCTTGTAGGCCTg                                                                                             >  2:887400/1‑150 (MQ=255)
                                                          tttCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATa                                                                                          <  2:1495378/150‑1 (MQ=255)
                                                              cTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCg                                                                                      >  2:1304416/1‑150 (MQ=255)
                                                                 aCGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAg                                                                                   >  1:1046364/1‑150 (MQ=255)
                                                                 aCGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCCGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAg                                                                                   >  1:556451/1‑150 (MQ=255)
                                                                  cGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAgc                                                                                  >  1:101590/1‑150 (MQ=255)
                                                                  cGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAgc                                                                                  <  1:1063317/150‑1 (MQ=255)
                                                                   gAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGACGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAgcg                                                                                 <  1:1016235/150‑1 (MQ=255)
                                                                              aCTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGt                                                                      >  1:571611/1‑150 (MQ=255)
                                                                               cTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGtt                                                                     >  1:630857/1‑150 (MQ=255)
                                                                               cTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGtt                                                                     >  2:1113969/1‑150 (MQ=255)
                                                                                tGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGttt                                                                    <  1:727356/150‑1 (MQ=255)
                                                                                 gCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGtttt                                                                   >  1:1494265/1‑150 (MQ=255)
                                                                                 gCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGtttt                                                                   >  2:367393/1‑150 (MQ=255)
                                                                                  cATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTg                                                                  >  2:67515/1‑150 (MQ=255)
                                                                                   aTTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGc                                                                 <  2:45046/150‑1 (MQ=255)
                                                                                    tttGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCg                                                                <  1:1620941/150‑1 (MQ=255)
                                                                                    tttGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCg                                                                <  1:1614530/150‑1 (MQ=255)
                                                                                     ttGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGt                                                               <  1:1366203/150‑1 (MQ=255)
                                                                                       gTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGttt                                                             >  2:1143659/1‑150 (MQ=255)
                                                                                       gTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGttt                                                             >  1:672908/1‑150 (MQ=255)
                                                                                        gcGGCAGCAACAACTGTTAAAAACGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTg                                                            <  2:1307579/149‑1 (MQ=255)
                                                                                         gggCATCAACTACTGTTAAAAAAGTGCGCTTTGTTTATGCGGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTCTTGCGTTTGt                                                           <  1:237225/149‑1 (MQ=255)
                                                                                          ggCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGtc                                                          >  2:1105299/1‑150 (MQ=255)
                                                                                          ggCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGAGGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGtc                                                          <  1:1296811/150‑1 (MQ=255)
                                                                                           gcagcaACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGtca                                                         <  1:1246686/150‑1 (MQ=255)
                                                                                            cagcaACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGtcat                                                        <  2:1576132/150‑1 (MQ=255)
                                                                                               caacaaCTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAg                                                     <  2:397647/150‑1 (MQ=255)
                                                                                                 acaacTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGtc                                                   <  2:11805/150‑1 (MQ=255)
                                                                                                 acaacTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGta                                                   >  2:91804/1‑149 (MQ=255)
                                                                                                    acTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAACCGCCTTATCCGGCCTACAAAATCGTGCTAATTCAAAATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCTTTTGTCATCAGTCTcc                                                >  1:382670/1‑150 (MQ=255)
                                                                                                      tGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGa                                              <  1:1105299/150‑1 (MQ=255)
                                                                                                      tGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCCGa                                              >  2:122955/1‑150 (MQ=255)
                                                                                                       tttAAAAAAGTGCGCTTTGTTTATGCCGGTTGCGGCGTAAACGCCTTATCCGGCCTACAAATTCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGAt                                             <  1:1004556/149‑1 (MQ=255)
                                                                                                       gTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGAt                                             >  1:616536/1‑150 (MQ=255)
                                                                                                        ttAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATg                                            >  2:1728270/1‑150 (MQ=255)
                                                                                                           aaaaaGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTa                                         >  2:405290/1‑150 (MQ=255)
                                                                                                           aaaaaGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTa                                         <  2:1610911/150‑1 (MQ=255)
                                                                                                           aaaaaGTGAGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCTATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGcga                                         <  1:592013/150‑3 (MQ=255)
                                                                                                            aaaaGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTAt                                        <  1:1572738/150‑1 (MQ=255)
                                                                                                             aaaGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTAtt                                       <  1:1396378/150‑1 (MQ=255)
                                                                                                                gTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAAt                                    <  1:1255615/150‑1 (MQ=255)
                                                                                                                 tGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATc                                   >  1:65307/1‑150 (MQ=255)
                                                                                                                 tGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATc                                   <  1:1336367/150‑1 (MQ=255)
                                                                                                                     cTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATACAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTa                               <  1:1337921/150‑1 (MQ=255)
                                                                                                                     cTTTGTTTATGCCAGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGAGTTTGTCATCAGTCTCCGATGCTATTAATCCTTa                               >  2:547887/1‑150 (MQ=255)
                                                                                                                        tGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAAt                            <  2:1698342/150‑1 (MQ=255)
                                                                                                                             atcccGGATGCGGCGTAAACGCCGTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCg                       <  2:350617/147‑1 (MQ=255)
                                                                                                                             aTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCg                       <  2:393165/150‑1 (MQ=255)
                                                                                                                              tGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGc                      <  2:630857/150‑1 (MQ=255)
                                                                                                                              ggCCGGATAAGGCGTTTACGCCGCATCCGGCATACAAAATCGTGCTAATTCAATATATTGCATAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGc                      >  2:338049/2‑150 (MQ=21)
                                                                                                                                     tGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTgg               <  2:1547678/150‑1 (MQ=255)
                                                                                                                                      gcggcgTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGc              <  2:1540468/150‑1 (MQ=255)
                                                                                                                                        ggcgTAAACTCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTGATCCTTAATTCCCCGCCCCCTGGCTa            >  2:73390/1‑150 (MQ=255)
                                                                                                                                         gcgTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTaa           <  1:1499350/150‑1 (MQ=255)
                                                                                                                                         gcgTAAACGCCTTATCCGGCCTACAAAATCGGGCTAATTCAATATATTGCAGACACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCAGCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCtca           >  1:136812/1‑148 (MQ=255)
                                                                                                                                          cgTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTaaa          >  2:15515/1‑150 (MQ=255)
                                                                                                                                          cgTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTaaa          <  2:1223109/150‑1 (MQ=255)
                                                                                                                                            tttaCGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGc                      >  1:1550120/4‑136 (MQ=255)
                                                                                                                                              aaCGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGc      <  1:1581985/150‑1 (MQ=255)
                                                                                                                                              aaCGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGc      <  2:195942/150‑1 (MQ=255)
                                                                                                                                              aaCGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGc      >  1:1276771/1‑150 (MQ=255)
                                                                                                                                                cGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGctc    >  2:302628/1‑150 (MQ=255)
                                                                                                                                                  ccTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGTCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCACCTGGCTAAAATGATCtt  >  2:180441/1‑150 (MQ=255)
                                                                                                                                                  ccTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCtt  >  2:1677098/1‑150 (MQ=255)
                                                                                                                                                  ccTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGACCCCTGGCTAAAATGCTCtt  >  2:294184/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ATGGGGTGACGGTCAGGACTCACCAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTT  >  NC_000913/2538408‑2538703

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