Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,538,556 T→G 7.0% 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 7.0% 286.0 / 8.0 95intergenic (+73/+116)yfeS/cysMWGR domain protein/cysteine synthase B (O‑acetylserine sulfhydrolase B)
Reads supporting (aligned to +/- strand):  ref base T (43/45);  new base G (4/3);  total (47/48)
Fisher's exact test for biased strand distribution p-value = 7.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

AATGGGGTGACGGTCAGGACTCACCAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCC  >  NC_000913/2538407‑2538705
                                                                                                                                                     |                                                                                                                                                     
aaTGGGGTGACGGTCAGGACTCACCAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCt                                                                                                                                                       <  2:643669/150‑1 (MQ=255)
   ggggTGACGGTCAGGACTCACCAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCCCGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTAt                                                                                                                                                    >  2:271658/1‑150 (MQ=255)
   ggggTGACGGTCAGGACTCACCAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTAt                                                                                                                                                    >  2:446313/1‑150 (MQ=255)
         aCGGTCAGGACTCACCAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCACCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                                                                                              <  1:389632/150‑1 (MQ=255)
             tCAGGACTCACCAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACa                                                                                                                                          <  2:521494/150‑1 (MQ=255)
              cAGGACTCACCAATACTTCAAAAAATGATTGAGGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTTACGCCGCATCCGGCATaaac                                                                                                                                         <  1:615200/150‑4 (MQ=38)
                 gACTCACCAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAt                                                                                                                                      <  2:1824113/150‑1 (MQ=255)
                   cTCACCAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCAGTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCg                                                                                                                                    <  1:317023/150‑1 (MQ=255)
                     cACCAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTg                                                                                                                                  >  1:1855386/1‑150 (MQ=255)
                     cACCAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTg                                                                                                                                  <  1:814765/150‑1 (MQ=255)
                      aCCAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGGTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGc                                                                                                                                 <  2:1315893/150‑1 (MQ=255)
                      aCCAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGc                                                                                                                                 <  1:304224/150‑1 (MQ=255)
                      aCCAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGc                                                                                                                                 <  2:1404601/150‑1 (MQ=255)
                      aCCAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGc                                                                                                                                 <  1:1136273/150‑1 (MQ=255)
                       ccAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAGAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCt                                                                                                                                <  2:450406/150‑1 (MQ=255)
                       ccAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCt                                                                                                                                >  1:558/1‑150 (MQ=255)
                        cAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTa                                                                                                                               >  1:749719/1‑150 (MQ=255)
                         aaTACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTaa                                                                                                                              <  1:1038442/150‑1 (MQ=255)
                             cTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCa                                                                                                                          >  1:444515/1‑150 (MQ=255)
                              ttCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCaa                                                                                                                         >  2:1683269/1‑150 (MQ=255)
                              ttCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCaa                                                                                                                         >  1:1767128/1‑150 (MQ=255)
                              ttCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCaa                                                                                                                         >  2:910644/1‑150 (MQ=255)
                               tCAAAATATGATTGATTACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACATCTGTTAAAATAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAat                                                                                                                        <  2:1011791/150‑1 (MQ=255)
                               tCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAat                                                                                                                        <  1:145768/150‑1 (MQ=255)
                                aaaaaaaTGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAata                                                                                                                       <  1:1383184/149‑1 (MQ=255)
                                  aaaaaTGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCATTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAatata                                                                                                                     >  2:36861/1‑150 (MQ=255)
                                    aaaTGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGGGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATAtt                                                                                                                   <  1:1542752/150‑1 (MQ=255)
                                            gatgaCCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAATTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAAc                                                                                                           <  2:934790/150‑1 (MQ=255)
                                            gatgaCCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAAc                                                                                                           >  2:246471/1‑150 (MQ=255)
                                              tgaCCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCt                                                                                                         >  1:750337/1‑150 (MQ=255)
                                                aCCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTg                                                                                                       <  1:995304/150‑1 (MQ=255)
                                                 ccTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGt                                                                                                      <  1:1195097/150‑1 (MQ=255)
                                                    tACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAgg                                                                                                   >  2:1534728/1‑150 (MQ=255)
                                                      cGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGcc                                                                                                 <  2:59488/150‑1 (MQ=255)
                                                           tttCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATa                                                                                            <  2:846656/150‑1 (MQ=255)
                                                           tttCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGATTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATa                                                                                            <  1:435045/150‑1 (MQ=255)
                                                             tCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAg                                                                                          >  2:1685995/1‑150 (MQ=255)
                                                             tCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAg                                                                                          >  2:485433/1‑150 (MQ=255)
                                                              ccTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGc                                                                                         <  2:1479961/150‑1 (MQ=255)
                                                              ccTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGc                                                                                         >  1:1872308/1‑150 (MQ=255)
                                                               cTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCg                                                                                        >  2:1426394/1‑150 (MQ=255)
                                                                 cACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTTCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTa                                                                                      >  1:403992/1‑150 (MQ=255)
                                                                  aCGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAg                                                                                     >  1:452025/1‑150 (MQ=255)
                                                                   cGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAgc                                                                                    <  1:271658/150‑1 (MQ=255)
                                                                   cGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAgc                                                                                    >  2:182462/1‑150 (MQ=255)
                                                                          tAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAg                                                                             <  1:584787/150‑1 (MQ=255)
                                                                           aaaTACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCAGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAgg                                                                            >  2:1741050/1‑150 (MQ=255)
                                                                                cTGCATTTGTCGGCAGCAACAACTGATAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGtt                                                                       >  1:1267860/1‑150 (MQ=255)
                                                                                  gCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGtttt                                                                     >  2:1470507/1‑150 (MQ=255)
                                                                                     tttGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCg                                                                  <  2:1063374/150‑1 (MQ=255)
                                                                                        gTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAATTCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGttt                                                               <  1:1153235/150‑1 (MQ=255)
                                                                                            gcagcaACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGtca                                                           >  1:892089/1‑150 (MQ=255)
                                                                                            gaagcaACAACTGTTAAATAAGGGCGCTTTGTTTATGCAGGATGCGGCGTAAACGCCTTATACGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATTAGCGTAGCGCATCAGGCAGTTTTGCGTTTGtca                                                           <  2:1208491/148‑1 (MQ=18)
                                                                                              agcaACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGtcatc                                                         <  2:1319132/150‑1 (MQ=255)
                                                                                              agcaACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCATACAAAATCGTGATATTTCAATATATTGCAGAAACCTTTTAGGCCTGATAAGCGTAGCGCATCAGGAAGTTTTGCGtttttaata                                                         >  1:956224/1‑144 (MQ=18)
                                                                                               gcaACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGtcatca                                                        <  2:1026912/150‑1 (MQ=255)
                                                                                                caacaaCTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAATCGTAGCGCATCAGGCAGTTTTGCTTTTGtcatcat                                                       >  1:770278/1‑149 (MQ=255)
                                                                                                  acaacTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGtc                                                     >  1:988522/1‑150 (MQ=255)
                                                                                                        gTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGGTAATTCAATATATTGCAGAAACCTGGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGAt                                               >  2:313346/1‑150 (MQ=255)
                                                                                                        gTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGAt                                               <  2:1153996/150‑1 (MQ=255)
                                                                                                          tAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGc                                             <  2:610362/150‑1 (MQ=255)
                                                                                                           aaaaaaGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCt                                            >  2:1726989/1‑150 (MQ=255)
                                                                                                            aaaaaGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTa                                           <  1:1121719/150‑1 (MQ=255)
                                                                                                             aaaaGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTAt                                          <  1:1470507/150‑1 (MQ=255)
                                                                                                              aaaGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTAtt                                         <  2:1855386/150‑1 (MQ=255)
                                                                                                               aaGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTa                                        <  1:1609003/150‑1 (MQ=255)
                                                                                                               aaGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCGGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCCGTCTCCGATGCTATTa                                        <  1:360057/150‑1 (MQ=255)
                                                                                                                   ccccTTTGTTTATGCCGGAGGCGGCGTAAACGCAGTAACCGGCATACAAAATCGTGCAAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGGAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATcc                                    <  2:177925/147‑1 (MQ=18)
                                                                                                                    cgcTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCt                                   <  2:749719/150‑1 (MQ=255)
                                                                                                                    cgcTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCt                                   <  1:1394774/150‑1 (MQ=255)
                                                                                                                     gcTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAATTCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCtt                                  >  1:1440211/1‑150 (MQ=255)
                                                                                                                         tGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTTATTCAATATATTGCAGAACCATTGTAGGCCTGCTAAGCGTCGCGCATCCGGCAGTTTTTCGGTTGTAATCAGTCTCCGATGCTTTTAATCCttctat                              >  1:1183167/1‑146 (MQ=18)
                                                                                                                            ttATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATccc                           >  2:773402/1‑150 (MQ=255)
                                                                                                                            ttATGCCGGATGCGGCGTAAAAGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATccc                           <  2:1667234/150‑1 (MQ=255)
                                                                                                                               ggCCGGATAAGGCGTTTACGCCGCATCGGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATGAGTCTCCGATGCTATTAATCCTTAAATCCACGc                        >  2:790911/2‑150 (MQ=21)
                                                                                                                                 ccGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGAATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGccc                      >  1:995990/1‑150 (MQ=255)
                                                                                                                                   ggataaggcgTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGccccc                    >  1:630728/7‑150 (MQ=21)
                                                                                                                                   ggataaggcgTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGccccc                    >  2:35381/7‑150 (MQ=21)
                                                                                                                                    gataaggcgTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCg                         <  2:15740/139‑1 (MQ=21)
                                                                                                                                    gataaggcgTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCg                         >  1:15740/6‑144 (MQ=21)
                                                                                                                                     aTGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAATTCCCCGCCCCCTg                  <  1:887838/150‑1 (MQ=255)
                                                                                                                                     aTGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTg                  <  1:1819042/150‑1 (MQ=255)
                                                                                                                                        cggcgTAAACGCCTTATCCGGCCTACAAAATTGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCt               >  2:799542/1‑150 (MQ=255)
                                                                                                                                          gcgTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTaa             >  2:523776/1‑150 (MQ=255)
                                                                                                                                          gcgTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCACCTGGCTaa             >  1:1299547/1‑150 (MQ=255)
                                                                                                                                           cgTAAACGCCTTATCCGGCCTACAAAATCGTGATAATTAAATATATTGCAGAAACCTAGTAGACCTGATAAGAGTAGCGCATCAGGCAGTTTTGCGTTTGTCATAAGTCTACGATGCTATTAATCCTTAAATCCCCGCCCCCTGGGTaaa            >  1:1275904/1‑150 (MQ=25)
                                                                                                                                             tAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAAt          >  2:485614/1‑150 (MQ=255)
                                                                                                                                             tAAACGCCTTATCCGGCCTACAAAATCGGGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAAt          <  1:87511/150‑1 (MQ=255)
                                                                                                                                               aaCGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGc        <  2:722610/150‑1 (MQ=255)
                                                                                                                                               aaCGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGc        <  2:853061/150‑1 (MQ=255)
                                                                                                                                                aCGCCTTATCCGGCCTACAAAATCGTTCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGct       <  1:60610/150‑1 (MQ=255)
                                                                                                                                                aCGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGct       >  1:1700144/1‑150 (MQ=255)
                                                                                                                                                aCGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGct       >  2:1832883/1‑150 (MQ=255)
                                                                                                                                                aCGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCACTGGCTAAATTGct       >  1:162775/1‑150 (MQ=255)
                                                                                                                                                   ccTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCtt    >  1:1595255/1‑150 (MQ=255)
                                                                                                                                                   ccTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCAGATGCTATTAATCCTTAAATCCCCGCACCCTGGCTAAAATGCTCtt    >  2:1685010/1‑150 (MQ=255)
                                                                                                                                                     ttATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTcc  <  1:1534728/150‑1 (MQ=255)
                                                                                                                                                     ttATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTcc  <  1:774168/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AATGGGGTGACGGTCAGGACTCACCAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCC  >  NC_000913/2538407‑2538705

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: