Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,538,557 T→C 6.3% intergenic (+74/+115) 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,5570TC 6.3% 285.4 / 7.6 95intergenic (+74/+115)yfeS/cysMWGR domain protein/cysteine synthase B (O‑acetylserine sulfhydrolase B)
Reads supporting (aligned to +/- strand):  ref base T (43/46);  new base C (4/2);  total (47/48)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.03e-01

GGGGTGACGGTCAGGACTCACCAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCC  >  NC_000913/2538410‑2538706
                                                                                                                                                   |                                                                                                                                                     
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:814765/150‑1 (MQ=255)
                  cACCAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTg                                                                                                                                   >  1:1855386/1‑150 (MQ=255)
                   aCCAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGGTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGc                                                                                                                                  <  2:1315893/150‑1 (MQ=255)
                   aCCAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGc                                                                                                                                  <  1:304224/150‑1 (MQ=255)
                   aCCAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGc                                                                                                                                  <  1:1136273/150‑1 (MQ=255)
                   aCCAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGc                                                                                                                                  <  2:1404601/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                                                                                                                          >  2:910644/1‑150 (MQ=255)
                           ttCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCaa                                                                                                                          >  1:1767128/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                                                                                     >  2:182462/1‑150 (MQ=255)
                                                                cGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAgc                                                                                     <  1:271658/150‑1 (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:853061/150‑1 (MQ=255)
                                                                                                                                            aaCGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGc         <  2:722610/150‑1 (MQ=255)
                                                                                                                                             aCGCCTTATCCGGCCTACAAAATCGTTCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGct        <  1:60610/150‑1 (MQ=255)
                                                                                                                                             aCGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGct        >  2:1832883/1‑150 (MQ=255)
                                                                                                                                             aCGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGct        >  1:1700144/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:774168/150‑1 (MQ=255)
                                                                                                                                                  ttATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTcc   <  1:1534728/150‑1 (MQ=255)
                                                                                                                                                   tATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTccc  <  1:522349/150‑1 (MQ=255)
                                                                                                                                                   tATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTccc  <  1:485433/150‑1 (MQ=255)
                                                                                                                                                   tATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGATATTAATCCTTAAATCCCCGCCCCCTGCCTAAAATGCTCTTccc  >  2:444562/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GGGGTGACGGTCAGGACTCACCAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGCTCTTCCC  >  NC_000913/2538410‑2538706

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: