Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,538,550 A→T 6.7% intergenic (+67/+122) 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,5500AT 6.7% 268.2 / 7.5 89intergenic (+67/+122)yfeS/cysMWGR domain protein/cysteine synthase B (O‑acetylserine sulfhydrolase B)
Reads supporting (aligned to +/- strand):  ref base A (39/44);  new base T (4/2);  total (43/46)
Fisher's exact test for biased strand distribution p-value = 4.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.56e-01

CAATGGGGTGACGGTCAGGACTCACCAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGC  >  NC_000913/2538406‑2538699
                                                                                                                                                |                                                                                                                                                     
cAATGGGGTGACGGTCAGGACTCACCAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGcc                                                                                                                                                  >  2:1819042/1‑150 (MQ=255)
 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:1136273/150‑1 (MQ=255)
                       aCCAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGc                                                                                                                           <  1:304224/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                                                                                                                   >  1:1767128/1‑150 (MQ=255)
                               ttCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCaa                                                                                                                   >  2:1683269/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                                                                                   >  1:1872308/1‑150 (MQ=255)
                                                               ccTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGc                                                                                   <  2:1479961/150‑1 (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                             <  1:1394774/150‑1 (MQ=255)
                                                                                                                     cgcTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCt                             <  2:749719/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)
                                                                                                                                                |                                                                                                                                                     
CAATGGGGTGACGGTCAGGACTCACCAATACTTCAAAAAATGATTGATGACCTTACGGCGTTTCCTCACGAAAATTAAATACTGCATTTGTCGGCAGCAACAACTGTTAAAAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCGTTTGTCATCAGTCTCCGATGCTATTAATCCTTAAATCCCCGCCCCCTGGCTAAAATGC  >  NC_000913/2538406‑2538699

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: