Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,491,164 C→A 7.2% intergenic (‑25/+42) idnR ← / ← idnT DNA‑binding transcriptional dual regulator IdnR/L‑idonate/5‑ketogluconate/gluconate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,491,1640CA 7.2% 919.0 / 37.9 270intergenic (‑25/+42)idnR/idnTDNA‑binding transcriptional dual regulator IdnR/L‑idonate/5‑ketogluconate/gluconate transporter
Reads supporting (aligned to +/- strand):  ref base C (128/122);  new base A (8/12);  total (136/134)
Fisher's exact test for biased strand distribution p-value = 3.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 294 total aligned reads displayed.

GGCGATGCGCTCGCCTGTTTCCTTTGCCACCTTTTTCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCAACGTTTCACCAACCGTCAGATTAAA  >  NC_000913/4491017‑4491313
                                                                                                                                                   |                                                                                                                                                     
ggCGATGCGCTCGCCTGTTTCCTTTGCCACCTTTTTCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGaatctt                                                                                                                                                     <  2:535238/150‑1 (MQ=255)
    aTGCGCTCGCCTGTTTCCTTTGCCACCTTTTTCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGaatcttcaat                                                                                                                                                 <  1:903633/150‑1 (MQ=255)
     tGCGCTCGCCTGTTTCCTTTGCCACCTTTTTCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGaatcttcaata                                                                                                                                                >  2:1533889/1‑149 (MQ=255)
       cgcTCGCCTGTTTCCTTTGCCACCTTTTTCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGaatcttcaatctt                                                                                                                                              >  1:1201699/1‑150 (MQ=255)
           cGCCTGTTTCCTTTGCCACGTTTTTCGGGGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTcccc                                                                                                                                          <  2:1168504/150‑1 (MQ=255)
           cGCCTGTTTCCTTTGCCACCTTTTTCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTcccc                                                                                                                                          <  2:1272560/150‑1 (MQ=255)
           cGCCTGTTTCCTTTGCCACCTTTTTCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTcccc                                                                                                                                          <  1:1329654/150‑1 (MQ=255)
            gCCTGTTTCCTTTGCCACCTTTTTCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAAGATTGAAGATTccccc                                                                                                                                         <  2:1769634/150‑1 (MQ=255)
               tGTTTCCTTTGCCACCTTTTTCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGaa                                                                                                                                      >  2:71174/1‑150 (MQ=255)
               tGTTTCCTTTGCCACCTTTTTCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAAGATTGAAGATTCCCCCGac                                                                                                                                      <  2:1344688/150‑2 (MQ=255)
                gTTTCCTTTGCCACCTTTTTCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAAGATTGAAGATTCCCCCGacg                                                                                                                                     <  1:1869086/150‑3 (MQ=255)
                gTTTCCTTTGCCACCTTTTTCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAAGATTGAAGATCCCCCCGacg                                                                                                                                     <  1:1591269/150‑3 (MQ=255)
                 tttCCTTTGCCACCTTTTTCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAAc                                                                                                                                    <  1:33755/150‑1 (MQ=255)
                   tCCTTTGCCACCTTTTTCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTa                                                                                                                                  >  1:1972588/1‑150 (MQ=255)
                   tCCTTTGCCACCTTTTTCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAAGATTGAAGATTCCCCCGacgttg                                                                                                                                  <  1:1650814/150‑6 (MQ=255)
                    ccTTTGCCACCTTTTTCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTAc                                                                                                                                 <  1:889077/150‑1 (MQ=255)
                     cTTTGCCACCTTTTTCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACa                                                                                                                                >  1:1387877/1‑150 (MQ=255)
                     cTTTGCCACCTTTTTCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACa                                                                                                                                >  2:95815/1‑150 (MQ=255)
                     cTTTGCCACCTTTTTCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAAGATTGAAGATTCCCCCGacgttggg                                                                                                                                <  1:1475641/150‑8 (MQ=255)
                       ttGCCACCTTTTTCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATc                                                                                                                              <  1:1513491/150‑1 (MQ=255)
                        tGCCACCTTTTTCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATca                                                                                                                             >  1:640897/1‑150 (MQ=255)
                         gCCACCTTTTTCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATcac                                                                                                                            <  1:1616403/150‑1 (MQ=255)
                          ccACCTTTTTCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATcaca                                                                                                                           <  1:1421890/150‑1 (MQ=255)
                               tttttCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTcccc                                                                                                                                          >  1:1691601/1‑130 (MQ=255)
                               tttttCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTcccc                                                                                                                                          <  2:1691601/130‑1 (MQ=255)
                                 tttCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATAtttt                                                                                                                    <  2:1858199/150‑1 (MQ=255)
                                 tttCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATAtttt                                                                                                                    <  1:1089668/150‑1 (MQ=255)
                                 tttCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAACGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATAtttt                                                                                                                    <  2:1665708/150‑1 (MQ=255)
                                  ttCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATAttttt                                                                                                                   >  2:1580531/1‑149 (MQ=255)
                                   tCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTAt                                                                                                                  <  1:423641/150‑1 (MQ=255)
                                    cGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGaa                                                                                                                                      >  2:1734272/1‑129 (MQ=255)
                                     ggCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTAttg                                                                                                                >  2:373469/1‑149 (MQ=255)
                                     ggCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTc                                                                                                                <  1:242231/150‑1 (MQ=255)
                                      gCGAGCGGATATAACGACTCACGGTCGTTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCt                                                                                                               <  2:253947/150‑1 (MQ=255)
                                        gAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCttt                                                                                                             <  2:1063526/150‑1 (MQ=255)
                                         aGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTc                                                                                                            <  1:1762607/150‑1 (MQ=255)
                                               tataACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCa                                                                                                      >  2:150005/1‑150 (MQ=255)
                                                 taACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAAt                                                                                                    >  1:1561083/1‑150 (MQ=255)
                                                 taACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAAt                                                                                                    >  2:577333/1‑150 (MQ=255)
                                                  aaCGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATa                                                                                                   >  2:1619329/1‑150 (MQ=255)
                                                   aCGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAg                                                                                                          <  2:507136/142‑1 (MQ=255)
                                                         cACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGtt                                                                                            >  1:1746863/1‑150 (MQ=255)
                                                           cGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTaa                                                                                          >  1:106805/1‑150 (MQ=255)
                                                            ggTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGt                                                                                             <  1:538971/146‑1 (MQ=255)
                                                            ggTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAAt                                                                                         >  1:37347/1‑150 (MQ=255)
                                                             gTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATg                                                                                        <  1:1892477/150‑1 (MQ=255)
                                                              tCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGTGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATgg                                                                                       <  2:1824204/150‑1 (MQ=255)
                                                                 tttttGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTa                                                                                    >  2:1914689/1‑150 (MQ=255)
                                                                   tttGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAac                                                                                  <  2:889771/150‑1 (MQ=255)
                                                                        tACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGcc                                                                             >  2:1361534/1‑150 (MQ=255)
                                                                         aCGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGcac                                                                            >  1:1247304/1‑148 (MQ=255)
                                                                          cGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCa                                                                           >  2:665206/1‑150 (MQ=255)
                                                                          cGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCa                                                                           >  2:1875098/1‑150 (MQ=255)
                                                                          cGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCa                                                                           >  1:110191/1‑150 (MQ=255)
                                                                          cGCCAGACAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGAGTTAATGGCTAACACGcaca                                                                           >  1:1949096/1‑150 (MQ=255)
                                                                           gccagccaGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACCTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAg                                                                          >  2:1461127/1‑150 (MQ=255)
                                                                             cagccagcGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCa                                                                        <  1:1933298/150‑1 (MQ=255)
                                                                                ccagcGTAGCGATATCCTGTATAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAAc                                                                     <  2:1090132/150‑1 (MQ=255)
                                                                                ccagcGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAAc                                                                     >  1:1578060/1‑150 (MQ=255)
                                                                                 cagcGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAAcc                                                                    <  1:66057/150‑1 (MQ=255)
                                                                                      tAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCataa                                                               >  1:1527790/1‑150 (MQ=255)
                                                                                               ccTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATaa                                                      >  2:1962553/1‑150 (MQ=255)
                                                                                                cTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAg                                                     >  2:346739/1‑150 (MQ=255)
                                                                                                cTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAg                                                     >  1:928552/1‑150 (MQ=255)
                                                                                                cTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCCTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAg                                                     >  1:511252/1‑150 (MQ=255)
                                                                                                 tGTAAAGAAATTCTGTGATTGCTCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAgg                                                    <  2:547317/150‑1 (MQ=255)
                                                                                                 tGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCACAGCAAACCCATAAAAGAAAAAAgg                                                    >  2:1986872/1‑150 (MQ=255)
                                                                                                   tAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGt                                                  <  2:1640702/150‑1 (MQ=255)
                                                                                                    aaaGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGtt                                                 >  2:490458/1‑150 (MQ=255)
                                                                                                      aGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTc                                               <  1:1407317/150‑1 (MQ=255)
                                                                                                      aGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTc                                               <  2:844408/150‑1 (MQ=255)
                                                                                                        aaaTTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCa                                             <  1:577333/150‑1 (MQ=255)
                                                                                                           ttCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCa                                          >  2:1308709/1‑150 (MQ=255)
                                                                                                                  aTTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTc                                                                                                            <  1:1709212/77‑1 (MQ=255)
                                                                                                                  aTTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCa                                   <  1:1309191/150‑1 (MQ=255)
                                                                                                                   ttGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAg                                  <  1:346739/150‑1 (MQ=255)
                                                                                                                    tGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAgg                                 <  1:722031/150‑1 (MQ=255)
                                                                                                                           aaCTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGt                                <  1:765337/144‑1 (MQ=255)
                                                                                                                            aCTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCAAc                         >  1:1988211/1‑150 (MQ=255)
                                                                                                                             cTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCAACg                        >  1:1210554/1‑150 (MQ=255)
                                                                                                                             cTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCAGAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCAACg                        >  2:1019777/1‑150 (MQ=255)
                                                                                                                                  cccAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCAACGTTTCa                   >  2:216665/1‑150 (MQ=255)
                                                                                                                                   ccAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCAACGTTTCAc                  <  2:382309/150‑1 (MQ=255)
                                                                                                                                   ccAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCAACGTTTCAc                  <  1:95815/150‑1 (MQ=255)
                                                                                                                                   acAACGTCGCGGGAATCTTCAATCTACCCCCTAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCAACGTTTCAc                  <  1:847016/149‑1 (MQ=255)
                                                                                                                                    cAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCAACGTTTCAcc                 <  2:608960/150‑1 (MQ=255)
                                                                                                                                     aaCGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCAACGttt                     <  2:908275/145‑1 (MQ=255)
                                                                                                                                      gttttgggggAAGATTGAAGATTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCAACGTTTCACCaa               >  2:741040/6‑150 (MQ=255)
                                                                                                                                        ttCGGGGGAAGATTGAAGATTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCAACGTTTCACCAAcc             >  2:870806/2‑150 (MQ=255)
                                                                                                                                         tCGGGGGAAGATTGAAGATTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCAACGTTTCACCAACCg            <  1:315675/150‑1 (MQ=255)
                                                                                                                                          cGGGGGAAGATTGAAGATTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCAACGTTTCACCAACCGt           >  2:1110184/1‑150 (MQ=255)
                                                                                                                                          cGGGGGAAGATTGAAGATTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCAACGTTTCACCAACCGt           >  2:283354/1‑150 (MQ=255)
                                                                                                                                           gggggAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCAACGTTTCACCAACCGTc          <  2:1218483/150‑1 (MQ=255)
                                                                                                                                              ggAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCAACGTTTCACCAACCGTCAGa       <  1:1829511/150‑1 (MQ=255)
                                                                                                                                               gAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCAACGTTTCACCAACCGTCAGAt      >  2:1391717/1‑150 (MQ=255)
                                                                                                                                                 atcttcaatcttcCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCAACGTTTCACCAACCGTCAGATTa    <  2:642700/150‑1 (MQ=255)
                                                                                                                                                  tcttcaatcttcCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCAACGTTTCACCAACCGTCAGATTaa   >  1:386676/1‑150 (MQ=255)
                                                                                                                                                   cttcaatcttcCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCAACGTTTCACCAACCGTCAGATTaaa  >  1:1025704/1‑150 (MQ=255)
                                                                                                                                                   cttcaatcttcCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCAACGTTTCACCAACCGTCAGATTaaa  >  2:328713/1‑150 (MQ=255)
                                                                                                                                                   cttcaatcttcCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCAACGTTTCACCAACCGTCAGATTaaa  >  2:1809165/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GGCGATGCGCTCGCCTGTTTCCTTTGCCACCTTTTTCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCAACGTTTCACCAACCGTCAGATTAAA  >  NC_000913/4491017‑4491313

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