Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,987,105 A→G 5.6% intergenic (+29/‑431) yqeG → / → yqeH putative transporter/putative LuxR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,987,1050AG 5.6% 436.3 / 11.5 143intergenic (+29/‑431)yqeG/yqeHputative transporter/putative LuxR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (61/73);  new base G (2/6);  total (64/79)
Fisher's exact test for biased strand distribution p-value = 3.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.60e-01

Only 100 of 144 total aligned reads displayed.

ATCATGCCTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGATGCATAACTTAG  >  NC_000913/2986960‑2987254
                                                                                                                                                 |                                                                                                                                                     
atcatGCCTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAAt                                                                                                                                                   <  1:1550885/150‑1 (MQ=255)
 tcatGCCTACGCTGTCAGCGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACTCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTGTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAAtt                                                                                                                                                  >  1:92940/1‑150 (MQ=255)
 tcatGCCTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAAtt                                                                                                                                                  <  2:1421645/150‑1 (MQ=255)
      ccTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGaa                                                                                                                                             >  1:388286/1‑150 (MQ=255)
      ccTACGCTGTCAACGAATCTCAACCCGGCGCTTAAACCCTGGAGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGaa                                                                                                                                             <  1:930362/150‑1 (MQ=255)
       cTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGTAGCAATTTTGAAc                                                                                                                                            <  1:847988/150‑1 (MQ=255)
       cTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGTAGCAATTTTGAAc                                                                                                                                            <  2:1574394/150‑1 (MQ=255)
       cTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGTAGCAATTTTGAAc                                                                                                                                            <  2:1671640/150‑1 (MQ=255)
       cTACGCTGTCAACGTATCTCATCCCGGCGCTAAACCCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTAACCTGTGCGTAGCGGTAATGTTCTTTAGCTAATACCAGTTGTGGATGGCTCAAAATTGTAGCAATTTTGAAc                                                                                                                                            <  2:1444696/150‑1 (MQ=255)
        tACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGTAGCAATTTTGAAcc                                                                                                                                           <  2:1760989/150‑1 (MQ=255)
        tACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGAGTATCCTGTGCGGATCGGTATTGTTCTTTAGCTAATATCAGTTGTTGATGGTTCAAAATTGATACAATTTTGAAcc                                                                                                                                           >  2:88953/1‑150 (MQ=255)
         aCGCTGTCAACGTTTCTCATCCCGGTGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAAcca                                                                                                                                          >  2:1211617/1‑150 (MQ=255)
            cTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAAccattc                                                                                                                                       <  1:1825962/150‑1 (MQ=255)
               tCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAAccattccat                                                                                                                                    >  1:1324287/1‑150 (MQ=255)
                 aaCGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCAttt                                                                                                                                  >  2:977462/1‑150 (MQ=255)
                  aCGTATCTCATCCCGGCGCTTATACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCAtttt                                                                                                                                 <  1:1294388/150‑1 (MQ=255)
                    gTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGAGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTc                                                                                                                               <  1:1573540/150‑1 (MQ=255)
                      aTCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCat                                                                                                                             <  2:514803/150‑1 (MQ=255)
                       tctcATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCata                                                                                                                            >  2:477372/1‑150 (MQ=255)
                        ctcATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGGTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCataa                                                                                                                           <  1:335223/150‑1 (MQ=255)
                         tcATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCataat                                                                                                                          >  2:917250/1‑150 (MQ=255)
                         tcATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCataat                                                                                                                          >  1:1304783/1‑150 (MQ=255)
                              ccGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCa                                                                                                                     >  2:1497385/1‑150 (MQ=255)
                              ccGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCa                                                                                                                     <  1:1347809/150‑1 (MQ=255)
                               cGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCAt                                                                                                                    >  2:465348/1‑150 (MQ=255)
                                 tggcTTAAACCCTGGCGTTACATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCAttt                                                                                                                  <  1:651836/148‑1 (MQ=255)
                                    cTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGtt                                                                                                               >  1:28815/1‑150 (MQ=255)
                                        aaCCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGttatta                                                                                                           <  1:1404238/150‑1 (MQ=255)
                                         aCCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTATCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTtattat                                                                                                          <  1:1803315/150‑1 (MQ=255)
                                          cccTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTattatt                                                                                                         >  1:892126/1‑150 (MQ=255)
                                              ggCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCAc                                                                                                     <  1:1069509/150‑1 (MQ=255)
                                              ggCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATAGGTAATGTTCTTTAGCTAATATCAGTTTTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCAc                                                                                                     <  1:30953/150‑1 (MQ=255)
                                                  ttCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCTGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTa                                                                                                 <  1:1065694/150‑1 (MQ=255)
                                                      aTCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCAc                                                                                             >  2:1028389/1‑150 (MQ=255)
                                                      aTCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCAc                                                                                             <  1:363315/150‑1 (MQ=255)
                                                      aTCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATAAc                                                                                             >  1:736769/1‑150 (MQ=255)
                                                      aTCGGAAATCTGATTACGCTGATCGGGGGTATCCTGTGCGGATCGGTAATGTTCTTAAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCAc                                                                                             <  2:301707/150‑1 (MQ=255)
                                                         ggAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTtc                                                                                          >  1:157749/1‑150 (MQ=255)
                                                         ggAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTtc                                                                                          >  2:1433242/1‑150 (MQ=255)
                                                          gAAATCTGATTACGCTGATCGTGGGTATACTGTGCGTATCGGTAATGTTCTTTAGCTAATATCTGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATTACTtct                                                                                         <  1:118272/150‑1 (MQ=255)
                                                           aaaTCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTtctc                                                                                        <  2:792412/150‑1 (MQ=255)
                                                           aaaTCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTtctc                                                                                        <  2:538527/150‑1 (MQ=255)
                                                            aaTCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTtctct                                                                                       >  2:420105/1‑150 (MQ=255)
                                                              tCTGATTACGCTGATAGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCg                                                                                     <  2:1591832/150‑1 (MQ=255)
                                                               cTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATAGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGt                                                                                    <  1:1140436/150‑1 (MQ=255)
                                                                tGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTc                                                                                   >  2:833960/1‑150 (MQ=255)
                                                                 gATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCa                                                                                  >  2:176214/1‑150 (MQ=255)
                                                                  tttACGCTTATAGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCat                                                                                 <  1:327823/149‑1 (MQ=255)
                                                                  aTTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCat                                                                                 >  2:783860/1‑150 (MQ=255)
                                                                      cGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCAACCCTATCCATTATCTCGTCATATcc                                                                             >  1:1604136/1‑150 (MQ=255)
                                                                        cTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTa                                                                           >  2:1445181/1‑150 (MQ=255)
                                                                           aTCGTGTGTATACTGTGAGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTTGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATATTATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGta                                                                        <  1:1541245/150‑1 (MQ=255)
                                                                             cGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTATCTCGTCATATCCTAGtatt                                                                      >  1:1598714/1‑150 (MQ=255)
                                                                              gTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTACACCCTATCACTTCTCTCGTCATATCCTAGtattt                                                                     >  1:973809/1‑149 (MQ=255)
                                                                                    aTCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCGCCCTATCACTTCTCTCGTCATATCCTAGTATTATCagtt                                                               <  2:1750155/150‑1 (MQ=255)
                                                                                     tcatctgCGTAACGGTAATGTTCTGTAGCTAATCTCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATATATCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCagtta                                                              <  2:1206064/145‑1 (MQ=255)
                                                                                        tgtgCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTtc                                                                                          <  2:386650/119‑1 (MQ=255)
                                                                                         gtgCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAgttagtta                                                          <  2:1803121/150‑1 (MQ=255)
                                                                                          tgCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAATATTGCTACAATTTTGAACCATTCCATTTTTCATATTATTCATTTGATATTATTCCACCCTATCACTTCTCTCGTCATATGCTTGTTTTATCAGTTAGTTaa                                                         <  2:660109/150‑1 (MQ=255)
                                                                                          tgCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTaa                                                         <  2:1211635/150‑1 (MQ=255)
                                                                                           gCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTaaa                                                        >  2:1167400/1‑150 (MQ=255)
                                                                                            cGTATCGGTAATGTTCTTTAGCTAATATCAGTTTTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATATTATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAac                                                       <  2:419507/150‑1 (MQ=255)
                                                                                            cGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCGACCCTATCACTTCTATCGTCATATCCTAGTATTATCAGTTAGTTaaaa                                                       >  1:142329/1‑149 (MQ=255)
                                                                                              tATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTATCGTCATATCCTATTATTATCAGTTAGTTAAacac                                                     >  1:59658/1‑150 (MQ=255)
                                                                                                tCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAacacac                                                   <  1:1109752/150‑1 (MQ=255)
                                                                                                    tAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAag                                               >  2:319903/1‑150 (MQ=255)
                                                                                                    tAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAAAACACAAag                                               >  2:1820119/1‑150 (MQ=255)
                                                                                                       tGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAgaga                                            <  1:977462/150‑1 (MQ=255)
                                                                                                       tGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAgaga                                            <  1:144281/150‑1 (MQ=255)
                                                                                                         ttCTTTAGCTAATATCAGTTGTGGATTGTTCAAAATTGCTAGAATTTTGAACCATTCCATTTTTCATATTATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATc                                          <  2:1475563/150‑1 (MQ=255)
                                                                                                           cTTTAGCTTATATCAGTTGTGGATTGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACTCACAAAGAGATCTg                                        <  1:1281844/150‑1 (MQ=255)
                                                                                                           cTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTg                                        >  2:388092/1‑150 (MQ=255)
                                                                                                            tttAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGa                                       <  1:1715769/150‑1 (MQ=255)
                                                                                                              tAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAAc                                     <  2:923323/150‑1 (MQ=255)
                                                                                                               aGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACa                                    >  1:248120/1‑150 (MQ=255)
                                                                                                                   aaTATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCtt                                >  2:1508747/1‑150 (MQ=255)
                                                                                                                   aaTATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCtt                                <  1:1001355/150‑1 (MQ=255)
                                                                                                                     tatCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTaa                              <  2:1120884/150‑1 (MQ=255)
                                                                                                                      atCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTaaa                             >  1:1795192/1‑150 (MQ=255)
                                                                                                                       tCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAAt                            >  2:105067/1‑150 (MQ=255)
                                                                                                                          gTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAg                         <  2:1324287/150‑1 (MQ=255)
                                                                                                                          gTTGTGGATGGTTCAAAATTGCTACAAATTATAACCATTCCATTTTTCATAATATTAATTTGTTATTATTACACCCTATAACTTCTCTCGTCATATACTAGTATTATCAGTTAGTTAAACACAAAAAGAGATCTGAAAATAATaaaaaaa                         >  2:1640119/1‑149 (MQ=255)
                                                                                                                            tgtgGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGAc                       >  2:429229/1‑150 (MQ=255)
                                                                                                                            tgtgGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGATATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGAc                       <  1:17406/150‑1 (MQ=255)
                                                                                                                              tgGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTa                     <  1:571760/150‑1 (MQ=255)
                                                                                                                              tgGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTa                     >  2:825261/1‑150 (MQ=255)
                                                                                                                               ggATGGTTCAAAATTGCTTCAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCAACCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTAAGTTAAACACTCAAAGAGATCTGAACATCTTAAATATGACTAt                    >  1:103305/1‑150 (MQ=255)
                                                                                                                               ggATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTAt                    >  2:889871/1‑150 (MQ=255)
                                                                                                                                    ggtCAAAATTGCTACAATTTTGAACAATTCCATTTTTCATAATATTCATTTGTTATTATTACACCCTATGACTTCTCTCGTCATATCGTAGTTTTATCAGTTAGTTTAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAg               <  1:637743/148‑1 (MQ=255)
                                                                                                                                    gTTCAAAATTGTAGCAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAg               >  1:95018/1‑150 (MQ=255)
                                                                                                                                    gTTCAAAATTGTAGCAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAg               >  1:251954/1‑150 (MQ=255)
                                                                                                                                    gTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAg               >  2:1306257/1‑150 (MQ=255)
                                                                                                                                       cAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGAATATTTAAGATg            >  2:983857/1‑150 (MQ=255)
                                                                                                                                       cAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCAACCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATATTAAATAAGACTATTTAAGATg            >  1:891147/1‑150 (MQ=255)
                                                                                                                                           aTTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGATGCATa        <  1:422968/150‑1 (MQ=255)
                                                                                                                                           aTTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGATGCATa        <  1:388092/150‑1 (MQ=255)
                                                                                                                                            ttGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAGTCTTATTTAAGATgttcag       <  1:420105/150‑6 (MQ=255)
                                                                                                                                              gCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTAATATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAGTCTTATTTAAGATgttcagat     <  1:1670435/150‑8 (MQ=255)
                                                                                                                                               cTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGATGCATAACtt    >  1:782197/1‑150 (MQ=255)
                                                                                                                                                tACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGATGCATAACTTa   <  1:1445181/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
ATCATGCCTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGATGGTTCAAAATTGCTACAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGATGCATAACTTAG  >  NC_000913/2986960‑2987254

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