Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,906,828 A→T 5.8% intergenic (‑261/+25) bglG ← / ← phoU transcriptional antiterminator of the bgl operon/negative regulator of PhoR/PhoB two‑component regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,906,8280AT 5.8% 399.0 / 9.0 156intergenic (‑261/+25)bglG/phoUtranscriptional antiterminator of the bgl operon/negative regulator of PhoR/PhoB two‑component regulator
Reads supporting (aligned to +/- strand):  ref base A (66/80);  new base T (4/6);  total (70/86)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.85e-01

Only 100 of 159 total aligned reads displayed.

ATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATA  >  NC_000913/3906679‑3906976
                                                                                                                                                     |                                                                                                                                                    
ataACGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGATATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCACAATAATTa                                                                                                                                                      >  1:1437962/1‑150 (MQ=255)
  aaCGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTaat                                                                                                                                                    >  2:823996/1‑149 (MQ=255)
   aCGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGt                                                                                                                                                   <  1:1150837/150‑1 (MQ=255)
    cgcaCATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGtt                                                                                                                                                  >  1:155336/4‑150 (MQ=255)
    cGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGtt                                                                                                                                                  <  1:511390/150‑1 (MQ=255)
    cGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGtt                                                                                                                                                  >  1:1702073/1‑150 (MQ=255)
     gAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTa                                                                                                                                                 <  2:1690411/150‑1 (MQ=255)
         aTCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGAGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTgg                                                                                                                                             <  1:1675197/150‑1 (MQ=255)
           ccAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAACTTAATTATTGGGa                                                                                                                                           <  1:168801/150‑1 (MQ=255)
            cAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGAt                                                                                                                                          >  2:649530/1‑150 (MQ=255)
            cAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAACTTAATTATTGGGAt                                                                                                                                          <  2:1601251/150‑1 (MQ=255)
             aGGTTCGGAAATTAATTTAATTGCGGGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGTAATTTATAAAGTTATATATAACAAATCCCAATAACTTAATTATTGGGAtt                                                                                                                                         <  2:313695/150‑1 (MQ=255)
              ggTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGAttt                                                                                                                                        >  1:590390/1‑150 (MQ=255)
              ggTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGAttt                                                                                                                                        >  2:635943/1‑150 (MQ=255)
               tttCGGAAATTAATTTAATTTCGTGCAGAGAAAATAGCAATGCGCTAATGATAAAAATATGACCATGCTCGCAGTTATTTACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATACTTAAGTTATTGGGATTTg                                                                                                                                       <  1:1457638/149‑1 (MQ=255)
                 tCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAACTTAATTATTGGGATTTGtt                                                                                                                                     <  2:302153/150‑2 (MQ=255)
                  cGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCt                                                                                                                                    >  1:989088/1‑150 (MQ=255)
                    gAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTgg                                                                                                                                  >  1:1666352/1‑150 (MQ=255)
                    gAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTAGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATACAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTgg                                                                                                                                  <  2:1112776/150‑1 (MQ=255)
                      aattaattTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTg                                                                                                                                >  1:1297940/1‑150 (MQ=255)
                         taatttaattGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAat                                                                                                                             <  1:1528263/150‑1 (MQ=255)
                          aatttaattGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAatt                                                                                                                            <  2:346311/150‑1 (MQ=255)
                             ttaattGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAattatt                                                                                                                         <  2:1356139/150‑1 (MQ=255)
                              taattGCGTGCAGAGAAAATAGCAATGAGCTTTTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTTAGTTATTGGGATTTGTCTGGTGAATTAttt                                                                                                                        <  2:384665/150‑1 (MQ=255)
                               aattGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTg                                                                                                                       <  2:414216/150‑1 (MQ=255)
                               aattGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTg                                                                                                                       <  1:1066183/150‑1 (MQ=255)
                                 ttGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTc                                                                                                                     <  2:1139027/150‑1 (MQ=255)
                                   gCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGc                                                                                                                   >  2:507253/1‑150 (MQ=255)
                                   gCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGc                                                                                                                   <  2:1030992/150‑1 (MQ=255)
                                    cGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCt                                                                                                                  >  2:1631519/1‑150 (MQ=255)
                                    cGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTAAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGTTGAATTATTTGTCGCt                                                                                                                  >  2:774417/1‑150 (MQ=255)
                                     gTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTa                                                                                                                 >  1:1238845/1‑150 (MQ=255)
                                       gCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATc                                                                                                               >  1:929261/1‑150 (MQ=255)
                                          gagaAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCttt                                                                                                            >  1:952068/1‑150 (MQ=255)
                                              aaaTAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTcccc                                                                                                        <  1:641649/150‑1 (MQ=255)
                                                 tAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGcc                                                                                                     <  2:106882/150‑1 (MQ=255)
                                                   gCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTccc                                                                                                         <  2:382020/144‑1 (MQ=255)
                                                   gCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTccc                                                                                                         >  1:382020/1‑144 (MQ=255)
                                                   gCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATAAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCcag                                                                                                   >  2:699088/1‑150 (MQ=255)
                                                    cAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCcagc                                                                                                  >  1:651435/1‑150 (MQ=255)
                                                    cAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCcagc                                                                                                  <  1:684283/150‑1 (MQ=255)
                                                    aaaTTCGCTATTGATAAAAATATGACCATGCTCGCATTTATTAACTTTGTGTAATTTTAGTAATTTAAAAACTTATATAAAACAAATACCAATAATTAAGTAATTTGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCcagccagc                                                                                                  <  2:1183580/149‑1 (MQ=255)
                                                      aTGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCcagcag                                                                                                <  2:636641/150‑1 (MQ=255)
                                                        gcgcTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGtt                                                                                              >  2:658844/1‑150 (MQ=255)
                                                          tcTATTGATAAAAATATGAACATGCTCGCAGTTATTAAGTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATATTTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTa                                                                                            <  2:1214511/149‑1 (MQ=255)
                                                          gcTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGAATTTTa                                                                                            >  1:1320966/1‑150 (MQ=255)
                                                            tATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATc                                                                                          >  1:887179/1‑150 (MQ=255)
                                                             aTTGATAAAAATATGACCATGCTCTCAGTTATTAACTTTGTGTTATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATcc                                                                                         <  2:683287/150‑1 (MQ=255)
                                                             aTTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATcc                                                                                         <  2:338965/150‑1 (MQ=255)
                                                                 aTAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCt                                                                                     <  1:777578/150‑1 (MQ=255)
                                                                     aaaTATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTAATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGATATCTTTCCCCGCCAGCAGATTATCCAGatcatc                                                                                 >  1:821293/1‑150 (MQ=255)
                                                                        tatGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGcc                                                                              <  1:665763/150‑1 (MQ=255)
                                                                          tGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCAc                                                                            <  2:228517/150‑1 (MQ=255)
                                                                             ccATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACAGa                                                                         >  2:1060657/1‑150 (MQ=255)
                                                                                 gCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTg                                                                     >  2:922115/1‑150 (MQ=255)
                                                                                 gCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTg                                                                     >  2:361425/1‑150 (MQ=255)
                                                                                  cTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGa                                                                    >  2:488584/1‑150 (MQ=255)
                                                                                     gCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACgg                                                                 <  2:1016435/150‑1 (MQ=255)
                                                                                         ttattaACTTTGTGTAATTTTAGGATTTTTTATAGTTATATATAAAAAATGCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTAGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAAt                                                             <  1:1714153/150‑1 (MQ=255)
                                                                                              aaCTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGc                                                        >  1:1601188/1‑150 (MQ=255)
                                                                                               aCTTTGTGTAATTTTAGGAATTTCTAAAGTTATATATAACAAATCCCAAAAATTAAGTTATTGGGATTTGTCTGCTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGcc                                                       <  1:254974/150‑1 (MQ=255)
                                                                                               aCTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGc                                                                                      >  1:473129/1‑119 (MQ=255)
                                                                                                   tgtgTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCtt                                                   >  1:720334/1‑150 (MQ=255)
                                                                                                    gtgtAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTc                                                  >  1:142973/1‑150 (MQ=255)
                                                                                                    gtgtAATTTTAGGAATTTATAAAGTTATATATAACAAATCACAATAATTAAGTTATTGGGATTTGTCTGTTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCtta                                                  >  1:1072278/1‑149 (MQ=255)
                                                                                                      gtAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCAc                                                >  2:1314941/1‑150 (MQ=255)
                                                                                                      gtAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCAc                                                >  2:1687287/1‑150 (MQ=255)
                                                                                                       tAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCCTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACg                                               <  1:729546/150‑1 (MQ=255)
                                                                                                        aaTTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACgt                                              <  1:1505406/150‑1 (MQ=255)
                                                                                                             tAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGtagtag                                         <  2:590390/150‑1 (MQ=255)
                                                                                                             tAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGtagtag                                         >  2:893157/1‑150 (MQ=255)
                                                                                                               tgAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCAGCGACGTGACGGAAATCCTGCCCCTTCACGTAGTagaa                                       <  1:43894/149‑1 (MQ=255)
                                                                                                               tgAATTTATAAAGTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTAACTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTagaag                                      <  1:1130067/149‑1 (MQ=255)
                                                                                                               ggAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTTCCCCTTCACGTAGTagaa                                       >  1:553726/1‑150 (MQ=255)
                                                                                                               ggAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTagaa                                       >  2:123625/1‑150 (MQ=255)
                                                                                                               ggAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTagaa                                       >  2:479074/1‑150 (MQ=255)
                                                                                                                 aaTTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTagaaga                                     <  2:722412/150‑1 (MQ=255)
                                                                                                                   tttATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATa                                   <  1:1080590/150‑1 (MQ=255)
                                                                                                                    ttATAAAGTTATATATAACAAATCCCAATAATTAAGTTATAGGGATTTGTCTGGTGAATTATTTGTCGATATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGAACCTTCACGTAGTAGAAGATaa                                  >  1:526211/1‑150 (MQ=255)
                                                                                                                      ataAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAAc                                >  2:1481011/1‑150 (MQ=255)
                                                                                                                      ataAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAAc                                <  2:458511/150‑1 (MQ=255)
                                                                                                                           gTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTcaca                           <  2:1702073/150‑1 (MQ=255)
                                                                                                                           gTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTcaca                           >  1:435220/1‑150 (MQ=255)
                                                                                                                            ttATATATAACACATCCAAATAATTAAGTTATTGTGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCATCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACaa                          <  2:1530919/150‑1 (MQ=255)
                                                                                                                              atatataACAAATCCCAATAATTAAGTTATTGGGATTTGTCAGGTGAATTATATGTCGCTATCTTTCCCCGCCAGAATTTTATCCAGATCATAGCCACCGACGAGACGGACATCCTGCCCCTTCACGTAGTAGAAGATAAACTCTCAAat                        <  1:280027/150‑1 (MQ=255)
                                                                                                                                 tataACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATTGCCACCGACGTGACGGAAATCCTGCCGCTTCACGTAGAAGAAGATAAACTCACAAATAtt                     >  1:256418/1‑150 (MQ=255)
                                                                                                                                  ataACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTc                    >  2:992669/1‑150 (MQ=255)
                                                                                                                                   agaCAAATCCCAATAACTTAATTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAATACCTGCCCCTTCACGTTTTATAATATAAAATCACAAATTTTcg                   >  1:877855/3‑149 (MQ=255)
                                                                                                                                   agaCAAATCCCAATAACTTAATTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCt                   >  1:1424540/3‑150 (MQ=255)
                                                                                                                                    aaCAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTg                  <  2:1666352/150‑1 (MQ=255)
                                                                                                                                    aaCAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTg                  <  2:1593099/150‑1 (MQ=255)
                                                                                                                                    aaCAAATCCCAATAATTAAGTTATTGGGATTTGGCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTg                  <  2:1376036/150‑1 (MQ=255)
                                                                                                                                      cAATTCACAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCTCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCTCGTAGTAGAAGATAAACTCACAAATATTCTGGc                <  1:1051658/150‑1 (MQ=255)
                                                                                                                                        aaTCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAg              >  1:228630/1‑150 (MQ=255)
                                                                                                                                        aaTCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAg              >  1:788907/1‑150 (MQ=255)
                                                                                                                                         aTCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGc             >  1:1658831/1‑150 (MQ=255)
                                                                                                                                         aTCCCAATAACTTAATTATTGGGATTTGTCTGGTGAATTATTTGTCTCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCAATCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGc             >  1:481761/1‑150 (MQ=255)
                                                                                                                                             cAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCAGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCACCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTc         >  2:881681/1‑150 (MQ=255)
                                                                                                                                                taatTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGcc      <  2:582824/150‑1 (MQ=255)
                                                                                                                                                aaatTAAGTAATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCTCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGcc      <  1:291662/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ATAACGAACATCCAGGTTCGGAAATTAATTTAATTGCGTGCAGAGAAAATAGCAATGCGCTATTGATAAAAATATGACCATGCTCGCAGTTATTAACTTTGTGTAATTTTAGGAATTTATAAAGTTATATATAACAAATCCCAATAATTAAGTTATTGGGATTTGTCTGGTGAATTATTTGTCGCTATCTTTCCCCGCCAGCAGTTTATCCAGCTCATCGCCACCGACGTGACGGAAATCCTGCCCCTTCACGTAGTAGAAGATAAACTCACAAATATTCTGGCAGCGGTCGCCAATA  >  NC_000913/3906679‑3906976

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: