Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,519,225 T→A 8.2% intergenic (+20/+32) yfeD → / ← gltX putative DNA‑binding transcriptional regulator YfeD/glutamate‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,519,2250TA 8.2% 453.7 / 25.0 158intergenic (+20/+32)yfeD/gltXputative DNA‑binding transcriptional regulator YfeD/glutamate‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base T (71/74);  new base A (2/11);  total (73/85)
Fisher's exact test for biased strand distribution p-value = 2.18e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.92e-01

Only 100 of 164 total aligned reads displayed.

GGCGTATCGACAAATCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTGGAGACTGCCCCGCAC  >  NC_000913/2519076‑2519372
                                                                                                                                                     |                                                                                                                                                   
ggCGTATCGACAAATCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAtatt                                                                                                                                                     <  1:276631/150‑1 (MQ=255)
    tATCGACAAATCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTAtct                                                                                                                                                 >  2:2248309/1‑150 (MQ=255)
     aTCGACAAATCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTTGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTAtctc                                                                                                                                                >  1:1874792/1‑150 (MQ=255)
       cGACAAATCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCt                                                                                                                                              <  1:975585/150‑1 (MQ=255)
        gACAAATCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTg                                                                                                                                             >  1:457297/1‑150 (MQ=255)
           agtacatcaACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCg                                                                                                                                          <  2:583460/146‑1 (MQ=255)
           aaaTCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCg                                                                                                                                          <  1:1973918/150‑1 (MQ=255)
           aaaTCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCg                                                                                                                                          >  1:1268928/1‑150 (MQ=255)
           aaaTCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCg                                                                                                                                          <  1:1120941/150‑1 (MQ=255)
            aaTCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGt                                                                                                                                         >  1:1586353/1‑150 (MQ=255)
            aaTCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGt                                                                                                                                         <  1:1028625/150‑1 (MQ=255)
               catcaACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTa                                                                                                                                      >  2:2017178/1‑150 (MQ=255)
               catcaACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTa                                                                                                                                      <  2:528528/150‑1 (MQ=255)
               catcaACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGtttt                                                                                                                                      <  1:1624536/150‑2 (MQ=255)
                atcaACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTAt                                                                                                                                     >  2:2097508/1‑150 (MQ=255)
                 tcaACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATc                                                                                                                                    <  2:326737/150‑1 (MQ=255)
                  caACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGtttttta                                                                                                                                   <  1:1002876/150‑5 (MQ=255)
                   aaCGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATATATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCtt                                                                                                                                  <  2:2037163/150‑1 (MQ=255)
                   aaCGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGttttttat                                                                                                                                  <  1:83096/150‑6 (MQ=255)
                    aCGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGttttttatt                                                                                                                                 <  2:1657883/150‑7 (MQ=255)
                    aCGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGttttttatt                                                                                                                                 <  1:1314072/150‑7 (MQ=255)
                     cgcgATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTTTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGttttttattt                                                                                                                                <  1:1936451/150‑8 (MQ=39)
                           agataATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACg                                                                                                                          <  2:2271449/150‑1 (MQ=255)
                           agataATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACg                                                                                                                          >  1:1710326/1‑150 (MQ=255)
                           agataATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACg                                                                                                                          <  1:1971532/150‑1 (MQ=255)
                               aataTGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTaa                                                                                                                      >  2:492110/1‑150 (MQ=255)
                                atatGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAAt                                                                                                                     >  2:384242/1‑150 (MQ=255)
                                      cagacaTTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgc                                                                                                               >  1:789384/1‑150 (MQ=255)
                                       agacaTTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgct                                                                                                              >  1:1357472/1‑150 (MQ=255)
                                        gacaTTGAGCAGGAAATAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgctg                                                                                                             <  2:1281947/150‑1 (MQ=255)
                                        gacaTTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgctg                                                                                                             <  2:1316176/150‑1 (MQ=255)
                                        gacaTTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgctg                                                                                                             >  1:1516441/1‑150 (MQ=255)
                                          caTTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGAt                                                                                                           >  1:519254/1‑150 (MQ=255)
                                           aTTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACTGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGAtt                                                                                                          <  2:1145119/150‑1 (MQ=255)
                                              gAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTc                                                                                                       >  2:2155659/1‑150 (MQ=255)
                                               aaaatgAAAAATTCGCCAAAATTTTGTCACATGAAAGAATTAGCAATTACCTCTCAACTATCGAAAATAACGAATCAGAATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTAGACGGTAATTACTGGTGATTTTcg                                                                                                      >  2:2143721/6‑150 (MQ=255)
                                                gCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTcgc                                                                                                     <  1:2248309/150‑1 (MQ=255)
                                                  aGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGt                                                                                                   <  1:2042940/150‑1 (MQ=255)
                                                     aaaaaGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCATTCa                                                                                                <  2:815564/150‑1 (MQ=255)
                                                       aaaGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGc                                                                                              <  2:734327/150‑1 (MQ=255)
                                                          gTCGCTCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGACCTTCTCGATATTGACGAATCAGCATAAATAAAATACGGTAGGTTATTATTTCGTGCAGTTTATCTTTTTACACGCTAATTACTTCTGATTTTCGAGTTCTGTAAt                                                                                           <  2:455609/150‑1 (MQ=255)
                                                            cgcgCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATaa                                                                                         >  2:150659/1‑150 (MQ=255)
                                                              cgcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATaaaa                                                                                       >  2:1541948/1‑150 (MQ=255)
                                                                cAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATc                                                                                     >  2:1197612/1‑150 (MQ=255)
                                                                 gaTTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATcc                                                                                    <  2:123670/149‑1 (MQ=255)
                                                                  aTTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCACGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCa                                                                                   <  2:1952807/150‑1 (MQ=255)
                                                                        gTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCtt                                                                             >  1:2284526/1‑150 (MQ=255)
                                                                          cACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttg                                                                           >  1:1668731/1‑150 (MQ=255)
                                                                           aCGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgt                                                                          <  1:1003239/150‑1 (MQ=255)
                                                                           aCGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgt                                                                          <  1:819868/150‑1 (MQ=255)
                                                                             gTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgttg                                                                        >  2:342463/1‑150 (MQ=255)
                                                                             gTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgttg                                                                        <  2:1947126/150‑1 (MQ=255)
                                                                              tGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGa                                                                       <  2:1847721/150‑1 (MQ=255)
                                                                               gAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGAt                                                                      <  2:1752684/150‑1 (MQ=255)
                                                                                aaGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATa                                                                     <  2:650369/150‑1 (MQ=255)
                                                                                  ggAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCAGTTTATCTTATTACACGCTAATTACTGCTGATTTTCGATTTCAGCAATAAAATCCAGCGCTTTGTTGATaag                                                                   >  1:662758/1‑148 (MQ=255)
                                                                                    aaTTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCt                                                                 >  1:1305268/1‑150 (MQ=255)
                                                                                     aTTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTc                                                                <  1:1729650/150‑1 (MQ=255)
                                                                                     aTTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTc                                                                <  2:173445/150‑1 (MQ=255)
                                                                                      ttCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCg                                                               <  2:727728/150‑1 (MQ=255)
                                                                                      ttCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCg                                                               <  2:1722643/150‑1 (MQ=255)
                                                                                        cGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGAt                                                             <  1:2474971/150‑1 (MQ=255)
                                                                                        cGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGAt                                                             >  1:345355/1‑150 (MQ=255)
                                                                                                ccTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCggg                                                     <  1:1348730/150‑1 (MQ=255)
                                                                                                  tCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTc                                                   <  2:659191/150‑1 (MQ=255)
                                                                                                  tCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTc                                                   >  1:103424/1‑150 (MQ=255)
                                                                                                       cGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTAcc                                              <  2:103424/150‑1 (MQ=255)
                                                                                                       cGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTAcc                                              >  2:130568/1‑150 (MQ=255)
                                                                                                            cGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTg                                         >  2:1609777/1‑150 (MQ=255)
                                                                                                            cGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTg                                         <  2:519254/150‑1 (MQ=255)
                                                                                                            cGATATTGACGAATCAGCATAAATAAAAAACCGCAGGATATTATATCATGCCGTTAATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCACCAATAAAATCCAGACATTTGTTGATACGCTCGATACTCCGGCTCTTACCAATTg                                         >  2:947619/1‑150 (MQ=255)
                                                                                                             gATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGc                                        >  2:126331/1‑150 (MQ=255)
                                                                                                              atatTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCg                                       <  1:140747/150‑1 (MQ=255)
                                                                                                              atatTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCg                                       <  1:59982/150‑1 (MQ=255)
                                                                                                                atTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTg                                     <  1:2130994/150‑1 (MQ=255)
                                                                                                                   gACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAAc                                  >  2:1643688/1‑150 (MQ=255)
                                                                                                                   gACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAAc                                  <  2:2064788/150‑1 (MQ=255)
                                                                                                                    aCGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACg                                 >  1:1259776/1‑150 (MQ=255)
                                                                                                                     cGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACgg                                >  1:82761/1‑150 (MQ=255)
                                                                                                                     cGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACgg                                <  2:259051/150‑1 (MQ=255)
                                                                                                                       aaTCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTGAGCAATAATATCCTGCGCTTTGTTTTTTCGCTCGATACTTCGTGTCTTACCAATTTCGTGAACGGTa                              <  2:656270/150‑1 (MQ=255)
                                                                                                                       aaTCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTa                              <  1:2226719/150‑1 (MQ=255)
                                                                                                                           aGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACAt                          >  2:809476/1‑150 (MQ=255)
                                                                                                                           aGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACAt                          >  2:1153699/1‑150 (MQ=255)
                                                                                                                            gCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATc                         >  2:1759762/1‑150 (MQ=255)
                                                                                                                            gCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATc                         >  1:2064602/1‑150 (MQ=255)
                                                                                                                               taaataaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAg                      >  1:2097954/1‑150 (MQ=255)
                                                                                                                                aaataaaAAACGGCAGGATATTAACTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAAAGGTACCATCCAGt                     >  2:1724971/1‑150 (MQ=255)
                                                                                                                                  ataaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGc                   >  1:1480972/1‑150 (MQ=255)
                                                                                                                                  ataaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGc                   <  1:513984/150‑1 (MQ=255)
                                                                                                                                   taaaAAACGGCAGGATATTATCTCCTGCCGTTTATGTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGGTTTGTTGATACGCTCGATACTGAGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCt                  >  1:2408575/1‑150 (MQ=255)
                                                                                                                                   taaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCt                  >  2:78819/1‑150 (MQ=255)
                                                                                                                                   taaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTATCGCGTTCAGCAATAAAATCCAGCGCTATGTTGATACGCTCGATACTGCGGGTCTTACCAATTCCGTGAACGCTACCATCCAGTggt                  >  2:1078157/1‑148 (MQ=255)
                                                                                                                                    gataaaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTg                 >  1:116474/4‑150 (MQ=255)
                                                                                                                                       aaaCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTGgag              <  1:492110/150‑1 (MQ=255)
                                                                                                                                            gCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTGGAGACTGc         <  1:1566136/150‑1 (MQ=255)
                                                                                                                                            gCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTAACAATTGCGTGAACGGTAACATCCAGTGCTGGAGACTGc         >  1:1053000/1‑150 (MQ=255)
                                                                                                                                              aGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTGGAGACTGccc       <  2:1445905/150‑1 (MQ=255)
                                                                                                                                               ggATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTGGAGACTGcccc      <  1:1735874/150‑1 (MQ=255)
                                                                                                                                                   attatCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTGGAGACTGCCCCGCAc  <  1:2115038/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GGCGTATCGACAAATCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTGGAGACTGCCCCGCAC  >  NC_000913/2519076‑2519372

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