Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,491,171 C→A 6.8% intergenic (‑32/+35) 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,1710CA 6.8% 916.4 / 35.3 270intergenic (‑32/+35)idnR/idnTDNA‑binding transcriptional dual regulator IdnR/L‑idonate/5‑ketogluconate/gluconate transporter
Reads supporting (aligned to +/- strand):  ref base C (128/123);  new base A (8/11);  total (136/134)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

Only 100 of 299 total aligned reads displayed.

TGCGCTCGCCTGTTTCCTTTGCCACCTTTTTCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCAACGTTTCACCAACCGTCAGATTAAAATACCCT  >  NC_000913/4491022‑4491320
                                                                                                                                                     |                                                                                                                                                     
tGCGCTCGCCTGTTTCCTTTGCCACCTTTTTCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGaatcttcaatc                                                                                                                                                       <  2:1820658/150‑1 (MQ=255)
tGCGCTCGCCTGTTTCCTTTGCCACCTTTTTCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGaatcttcaatc                                                                                                                                                       >  2:153268/1‑150 (MQ=255)
 gcgcTCGCCTGTTTCCTTTGCCACCTTTTTCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGaatcttcaatct                                                                                                                                                      >  1:73678/1‑150 (MQ=255)
      cGCCTGTTTCCTTTGCCACGTTTTTCGGGGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTcccc                                                                                                                                                 <  2:1168504/150‑1 (MQ=255)
      cGCCTGTTTCCTTTGCCACCTTTTTCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTcccc                                                                                                                                                 <  2:1272560/150‑1 (MQ=255)
      cGCCTGTTTCCTTTGCCACCTTTTTCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTGCAATCTTcccc                                                                                                                                                 <  2:369326/150‑1 (MQ=255)
       gCCTGTTTCCTTTGCCACCTTTTTCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTccccc                                                                                                                                                <  2:629289/150‑1 (MQ=255)
          tGTTTCCTTTGCCACCTTTTTCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGaa                                                                                                                                             >  1:196477/1‑150 (MQ=255)
          tGTTTCCTTTGCCACCTTTTTCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAAGATTGAAGATTCCCCCGac                                                                                                                                             <  1:507596/150‑2 (MQ=255)
          tGTTTCCTTTGCCACCTTTTTCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAAGATTGAAGATTCCCCCGac                                                                                                                                             <  2:1344688/150‑2 (MQ=255)
           gTTTCCTTTGCCACCTTTTTCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGaaa                                                                                                                                            <  1:1882508/150‑1 (MQ=255)
           gTTTCCTTTGCCACCTTTTTCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAAGATTGAAGATTCCCCCGacg                                                                                                                                            <  2:1897117/150‑3 (MQ=255)
            tttCCTTTGCCACCTTTTTCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAAc                                                                                                                                           <  1:33755/150‑1 (MQ=255)
              tCCTTTGCCACCTTTTTCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAAGATTGAAGATTCCCCCGacgttg                                                                                                                                         <  1:1650814/150‑6 (MQ=255)
               ccTTTGCCACCTTTTTCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTAc                                                                                                                                        <  1:889077/150‑1 (MQ=255)
                  ttGCCACCTTTTTCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATc                                                                                                                                     <  2:1840236/150‑1 (MQ=255)
                          tttttCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTcccc                                                                                                                                                 <  2:1691601/130‑1 (MQ=255)
                           ttttCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATAttt                                                                                                                            >  2:1753383/1‑150 (MQ=255)
                           ttttCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATAttt                                                                                                                            <  1:132914/150‑1 (MQ=255)
                            tttCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATAtttt                                                                                                                           <  1:1089668/150‑1 (MQ=255)
                            tttCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAACGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATAtttt                                                                                                                           <  2:1665708/150‑1 (MQ=255)
                               cGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGaa                                                                                                                                             >  2:1734272/1‑129 (MQ=255)
                                ggCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTAtta                                                                                                                       >  1:1684208/1‑149 (MQ=255)
                                ggCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTc                                                                                                                       <  1:242231/150‑1 (MQ=255)
                                   gAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCttt                                                                                                                    <  2:1063526/150‑1 (MQ=255)
                                    aGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTc                                                                                                                   <  1:1762607/150‑1 (MQ=255)
                                       ggATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGt                                                                                                                >  2:337670/1‑150 (MQ=255)
                                        gATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTg                                                                                                               >  1:1278641/1‑150 (MQ=255)
                                        gATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTg                                                                                                               <  2:925143/150‑1 (MQ=255)
                                         atatAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCATTCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGc                                                                                                              <  2:1866190/150‑1 (MQ=255)
                                              aCGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAg                                                                                                                 <  2:507136/142‑1 (MQ=255)
                                               cGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACg                                                                                                        >  1:1404717/1‑150 (MQ=255)
                                                gACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACgg                                                                                                       <  1:1220687/150‑1 (MQ=255)
                                                gACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAAAACgg                                                                                                       <  1:694807/150‑1 (MQ=255)
                                                  cTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCg                                                                                                     >  1:518136/1‑150 (MQ=255)
                                                   tCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGt                                                                                                    <  2:892441/150‑1 (MQ=255)
                                                    cACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGtt                                                                                                   >  1:1746863/1‑150 (MQ=255)
                                                      cGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTaa                                                                                                 >  1:642700/1‑150 (MQ=255)
                                                      cGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTaa                                                                                                 >  1:106805/1‑150 (MQ=255)
                                                       ggTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAAt                                                                                                >  1:1569917/1‑150 (MQ=255)
                                                          cATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGc                                                                                             >  2:1458340/1‑150 (MQ=255)
                                                              tttGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAac                                                                                         <  1:1848182/150‑1 (MQ=255)
                                                                   tACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGcc                                                                                    >  2:1361534/1‑150 (MQ=255)
                                                                     cGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCa                                                                                  >  1:1415997/1‑150 (MQ=255)
                                                                       ccagccagCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGc                                                                                >  1:1515451/1‑150 (MQ=255)
                                                                        cagccagcGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCa                                                                               >  1:552405/1‑150 (MQ=255)
                                                                         agccagcGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCaa                                                                              <  1:350642/150‑1 (MQ=255)
                                                                           ccagcGTAGCGATATCCTGTATAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAAc                                                                            <  2:1090132/150‑1 (MQ=255)
                                                                                    cGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAg                                                                   <  1:1600008/150‑1 (MQ=255)
                                                                                         tCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATa                                                              >  2:970566/1‑150 (MQ=255)
                                                                                          ccTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTAc                                                                                                                                        <  1:1478233/75‑1 (MQ=255)
                                                                                          ccTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTAc                                                                                                                                        >  2:1478233/1‑75 (MQ=255)
                                                                                           cTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAg                                                            >  1:928552/1‑150 (MQ=255)
                                                                                           cTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAg                                                            >  1:377655/1‑150 (MQ=255)
                                                                                              tAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGt                                                         <  2:653526/150‑1 (MQ=255)
                                                                                              tAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGt                                                         >  2:884676/1‑150 (MQ=255)
                                                                                               aaaGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGtt                                                        >  2:490458/1‑150 (MQ=255)
                                                                                               aaaGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGtt                                                        <  2:552296/150‑1 (MQ=255)
                                                                                                aaGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGttt                                                       >  2:722031/1‑150 (MQ=255)
                                                                                                 aGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTc                                                      >  2:26973/1‑150 (MQ=255)
                                                                                                 aGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTc                                                      <  2:582190/150‑1 (MQ=255)
                                                                                                   aaaTTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCa                                                    >  1:608162/1‑150 (MQ=255)
                                                                                                      ttCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCa                                                 >  2:1308709/1‑150 (MQ=255)
                                                                                                          gtgATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGt                                             >  2:216636/1‑150 (MQ=255)
                                                                                                            gATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTcc                                           >  2:1808769/1‑150 (MQ=255)
                                                                                                            gATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTcc                                           >  1:1137146/1‑150 (MQ=255)
                                                                                                             aTTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTc                                                                                                                   <  1:1709212/77‑1 (MQ=255)
                                                                                                                 cgcATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTa                                      <  2:53642/150‑1 (MQ=255)
                                                                                                                   cATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACg                                    <  2:1972588/150‑1 (MQ=255)
                                                                                                                   cATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACg                                    <  1:1361534/150‑1 (MQ=255)
                                                                                                                   cATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACg                                    <  1:1914689/150‑1 (MQ=255)
                                                                                                                     tAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTATTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCa                                  >  2:6745/1‑150 (MQ=255)
                                                                                                                      aaCTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGt                                       >  2:765337/1‑144 (MQ=255)
                                                                                                                       aCTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCAAc                                >  1:1207659/1‑150 (MQ=255)
                                                                                                                       aCTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCAAc                                >  1:1988211/1‑150 (MQ=255)
                                                                                                                           ctcCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCAACGttt                            >  1:41738/1‑150 (MQ=255)
                                                                                                                             cccAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAAAAAGGGTTTCCATCACCGTCCAGGTACGCAACGTTTCa                          >  2:1637362/1‑150 (MQ=255)
                                                                                                                              ccAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACTCAACGTTTCAc                         >  2:1617678/1‑150 (MQ=255)
                                                                                                                              ccAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCAACGTTTCAc                         <  2:382309/150‑1 (MQ=255)
                                                                                                                              ccAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCAACGTTTCAc                         <  1:95815/150‑1 (MQ=255)
                                                                                                                              ccAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCAACGTTTCAc                         <  1:1589876/150‑1 (MQ=255)
                                                                                                                              ccAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCAACGTTTCAc                         <  2:86787/150‑1 (MQ=255)
                                                                                                                              acAACGTCGCGGGAATCTTCAATCTACCCCCTAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCAACGTTTCAc                         <  1:847016/149‑1 (MQ=255)
                                                                                                                               cAACGTCGGGGGAATATTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCAACGTTTCAcc                        <  1:181842/150‑1 (MQ=255)
                                                                                                                                  cGTCGGGGGAATCTTCAATCTTCCCCAGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCAACGTTTCACCAAc                     <  2:1949096/150‑1 (MQ=255)
                                                                                                                                   gTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCAACGTTTCACCAAcc                    <  2:237387/150‑1 (MQ=255)
                                                                                                                                    tCGGGGGAAGATTGAAGATTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCAACGTTTCACCAACCg                   >  2:315675/1‑150 (MQ=255)
                                                                                                                                    tCGGGGGAAGATTGAAGATTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCAACGTTTCACCAACCg                   <  1:315675/150‑1 (MQ=255)
                                                                                                                                     cGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCAACGTTTCACCAACCGt                  >  1:1936589/1‑150 (MQ=255)
                                                                                                                                     cGGGGGAAGATTGAAGATTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCAACGTTTCACCAACCGt                  >  2:1110184/1‑150 (MQ=255)
                                                                                                                                      gggggAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCAACGTTTCACCAACCGTc                 <  2:1569917/150‑1 (MQ=255)
                                                                                                                                      gggggAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCAACGTTTCACCAACCGTc                 <  2:1218483/150‑1 (MQ=255)
                                                                                                                                       ggggAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCAACGTTTCACCAACCGTCa                >  1:1439456/1‑150 (MQ=255)
                                                                                                                                          gAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCAACGTTTCACCAACCGTCAGAt             <  1:153268/150‑1 (MQ=255)
                                                                                                                                          gAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCAACGTTTCACCAACCGTCAGAt             >  2:1391717/1‑150 (MQ=255)
                                                                                                                                              cttcaatcttcCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCAACGTTTCACCAACCGTCAGATTaaa         >  1:1025704/1‑150 (MQ=255)
                                                                                                                                               ttcaatcttcCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCAACGTTTCACCAACCGTCAGATTaaaa        >  1:797749/1‑150 (MQ=255)
                                                                                                                                               ttcaatcttcCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCAACGTTTCACCAACCGTCAGATTaaaa        >  1:899551/1‑150 (MQ=255)
                                                                                                                                                 caatcttcCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCAACGTTTCACCAACCGTCAGATTAAAATa      >  1:1092648/1‑150 (MQ=255)
                                                                                                                                                     cttcCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCAACGTTTCACCAACCGTCAGATTAAAATACCCt  <  1:996298/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGCGCTCGCCTGTTTCCTTTGCCACCTTTTTCGGCGAGCGGATATAACGACTCACGGTCATTTTTGTTACGCCAGCCAGCGTAGCGATATCCTGTAAAGAAATTCTGTGATTGCGCATAACTTCTCCCAACGTCGGGGGAATCTTCAATCTTCCCCCGAAACTACATCACAATATTTTATTCTTTCAGTGCAATACGGCGTTAATGGCTAACACGCCCAGCAAACCCATAATAGAAATAAGGGTTTCCATCACCGTCCAGGTACGCAACGTTTCACCAACCGTCAGATTAAAATACCCT  >  NC_000913/4491022‑4491320

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: