Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,668,977 C→T intergenic (‑81/‑55) hcaR ← / → hcaE DNA‑binding transcriptional dual regulator HcaR/putative 3‑phenylpropionate/cinnamate dioxygenase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,668,9770CT98.0% 167.1 / NA 51intergenic (‑81/‑55)hcaR/hcaEDNA‑binding transcriptional dual regulator HcaR/putative 3‑phenylpropionate/cinnamate dioxygenase subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (24/25);  minor base G (0/1);  total (24/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCGCTTATCCCTCACCAGTAACGGAACACCGACACAGTTTTCAAGATCGCGGATCTGGCTGCTTAACGAAGGCTGTGAGGTATGCAGTTTTTCCGCCGCACGGGTAAAGTTCAGTGCCTGCGCCACTGCGACGAAATAGCGTAAATGCCGTAGTTCCATCACCTTCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTATTTCACATATTAGCAACCAACCAGCAACATCCTTATGGCACAAAAATAGAAGGTCAATACATCTTATCTTTCAGGATTAAAAAATATGACCACACCCTCAGATTTGAACATTTACCAACTGATTGACACCCAAAATGGTCGGGTCACTCCGCGTATTTATACCGACCCGGATATTTACCAACTGGAGCTTGAGCGTATTTTCGGTCGTTGCTGGCTATTTCTCGCCCACGAAAGCCAGATCCCAAAACCCGGTGATTTCTTTAACACCTAC  >  NC_000913/2668737‑2669217
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                
tGCGCTTATCCCGAACCAGTAACGGAACACCGACACAGTTTTCAAGATCGCGGATCTGGCTGCTTAACGAAGGCTGTGAGGTATGCAGTTTTTCCGCCGCACGGGTAAAGTTCAGTGCCTGCGCCACTGCGACGAAATAGCGTAAATGCCGTAGTTCCATCACCTTCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTATTTCACATATTAGCAACCAACCAGCAACATCCTTATGGca                                                                                                                                                                                                                                         <  2:235435/250‑1 (MQ=255)
      tATCCCTCACCAGTAACGGAACACCGACACAGTTTTCAAGATCGCGGATCTGGCTGCTTAACGAAGGCTGTGAGGTATGCAGTTTTTCCGCCGCACGGGTAAAGTTCAGTGCCTGCGCCACTGCGACGAAATAGCGTAAATGCCGTAGTTCCATCACCTTCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTATTTCACATATTAGCAACCAACCAGCAACATTCTTATGGCACaaaa                                                                                                                                                                                                                                    <  1:101311/249‑1 (MQ=255)
       aTCCCGCACCAGTAGCGGAACACCGACACAGGGTTCAAGATCGCGGATCTGGCTGCGTAACGAAGGCTGTGAGGTATGCAGTTTTTCCGCCGCACGGGTAAAGTTCAGTGCCTGCGCCACTGCGACGAAATAGCGTAAATGCCGTAGTTCCATCACCTTCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTATTTCACATATTAGCAACCAACCAGCAACATTCTTATGGCACaaaaaga                                                                                                                                                                                                                                 <  2:326868/251‑3 (MQ=255)
        tCCCTCACCAGTAACGGAACACCGACACAGTTTTCAAGATCGCGGATCTGGCTGCTTAACGAAGGCTGTGAGGTATGCAGTTTTTCCGCCGCACGGGTAAAGTTCAGTGCCTGCGCCACTGCGACGAAATAGCGTAAATGCCGTAGTTCCATCACCTTCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTATTTCACATATTAGCAACCAACCAGCAACATTCTTATGGCACAAAAAt                                                                                                                                                                                                                                  >  1:655257/1‑249 (MQ=255)
                cAGTAACGGACCACCGACACAGTTTTCAAGATCGCGGATCTGGCTGCTTAACGAAGGCTGTGAGGTATGCAGTTTTTCCGCCGCACGGGTAAAGTTCAGTGCCTGCGCCACTGCGACGAAATAGCGTAAATGCCGTAGTTCCATCACCTTCCCCTTGTTATCGAAAACACGTCTCAGCTGGTAGAAATTAAATATTTCACATATTAGCAACCAACCAGCAACATTATTATGGCACAAAAATAGCAGGTCa                                                                                                                                                                                                                         >  2:117016/1‑250 (MQ=255)
                        gAACACCGACACAGTTTTCAAGATCGCGGATCTGGCTGCTTAACGAAGGCTGTGAGGTATGCAGTTTTTCCGCCGCACGGGTAAAGTTCAGTGCCTGCGCCACTGCGACGAAATAGCGTAAATGCCGTAGTTCCATCACCTTCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTATTTCACATATTAGCAACCAACCAGCAACATTCTTATGGCACAAAAATAGAAGGTCAATACatctt                                                                                                                                                                                                                >  1:689526/1‑251 (MQ=255)
                            acCGACACAGTTTTCAAGATCGCGGATCTGGCTGCTTAACGAAGGCTGTGAGGTATGCAGTTTTTCCGCCGCACGGGTAAAGTTCAGTGCCTGCGCCACTGCGACGAAATAGCGTAAATGCCGTAGTTCCATCACCTTCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTATTTCACATATTAGCAACCAACCAGCAACATTCTTATGGCACAAAAATAGAAGGTCAATACatcttatct                                                                                                                                                                                                            >  1:249114/1‑251 (MQ=255)
                                acacAGTTTTCAAGATCGCGGATCTGGCTGCTTAACGAAGGCTGTGAGGTATGCAGTTTTTCCGCCGCACGGGTAAAGTTCAGTGCCTGCGCCACTGCGACGAAATAGCGTAAATGCCGTAGTTCCATCACCTTCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTATTTCACATATTAGCAACCAACCAGCAACATTCTTATGGCACAAAAATAGAAGGTCAATACATCTTATCTTTCa                                                                                                                                                                                                        <  1:550618/251‑1 (MQ=255)
                                    aGTTTTCCAGATCGCGGATCTGGCTGCTTAACGAAGGCTGTGAGGTATGCAGTTTTTCCGCCGCACGGGTAAAGTTCAGTGCCTGCGCCACTGCGACGAAATAGCGTAAATGCCGTAGTTCCATCACCTTCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTATTTCACATATTAGCAACCAACCAGCAACATTCTTATGGCACAAAAATAGAAGGTCAATACATCTTATCTTTCAGGAt                                                                                                                                                                                                    <  2:7984/251‑1 (MQ=255)
                                        ttCAAGATCGCGGATCTGGCTGCTTAACGAAGGCTGTGAGGTATGCAGTTTTTCCGCCGCACGGGTAAAGTTCAGTGCCTGCGCCACTGCGACGAAATAGCGTAAATGCCGTAGTTCCATCACCTTCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTATTTCACATATTAGCAACCAACCAGCAACATTCTTATGGCACAAAAATAGAAGGTCAATACATCTTATCTTTCAGGATTaa                                                                                                                                                                                                 <  2:249114/250‑1 (MQ=255)
                                             gATCGCGGATCTGGCTGCTTAACGAAGGCTGTGAGGTATGCAGTTTTTCCGCCGCACGGGTAAAGTTCAGTGCCTGCGCCACTGCGACGAAATAGCGTAAATGCCGTAGTTCCATCACCTTCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTATTTCACATATTAGCAACCAACCAGCAACATTCTTATGGCACAAAAATAGAAGGTCAATACATCTGATCTTTCAGGATTAAAAAata                                                                                                                                                                                           <  1:435024/251‑1 (MQ=255)
                                              aTCGCGGATCTGGCTGTATAACGAAGGCTGTGAGGTAGGCAGTTTTTCCGCCGCACGGGTAAAGTGCAGTGCCTGCGCCACTGCGACGAAATAGCGTAAATGCCGTAGTTCCATCACCTTCCCCTTGTTAACGAAAAAACGTCGCAGCTGGTAGAAATTAACTATTGCACATATTAGCAACCAACCAGCAACATTCGTATGGAACAAAAATAGAAGGTCAATACATCGTATCTTTCAGGATGAAAAAatat                                                                                                                                                                                          <  2:416068/251‑1 (MQ=255)
                                              aTCGCGGATCTGGCTGCTTAACGAAGGCTGTGAGGTATGCAGTTTTTCCGCCGCACGGGTAAAGTTCAGTGCCTGCGCCACTGCGACGAAATAGCGTAAATGCCGTAGTTCCATCACCTTCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTATTTCACATATTAGCAACCAACCAGCAACATTCTTATGGCACAAAAATAGAAGGTCAATACATCTTATCTTTCAGGATTAAAAAatat                                                                                                                                                                                          <  1:604400/251‑1 (MQ=255)
                                                cgcgGATCTGGCTGCTTAACGAAGGCTGTGAGGTATGCAGTTTTTCCGCCGCACGGGTAAAGTTCAGTGCCTGCGCCACTGCGACGAAATAGCGTAAATGCCGTAGTTCCATCACCTTCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTATTTCACATATTAGCAACCAACCAGCAACATTCTTATGGCACAAAAATAGAAGGTCAATACATCTTATCTTTCAGGATTAAAAAATATg                                                                                                                                                                                         >  2:282471/1‑250 (MQ=255)
                                                  cgGATCTGGCTGCTTAACGAAGGCTGTGAGGTATGCAGTTTTTCCGCCGCACGGGTAAAGTTCAGTGCCTGCGCCACTGCGACGAAATAGCGTAAATGCCGTAGTTCCATCACCTTCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTATTTCACATATTAGCAACCAACCAGCAACATTCTTATGGCACAAAAATAGAAGGTCAATACATCTTATCTTTCAGGATTAAAAAATATGAcc                                                                                                                                                                                      <  1:219284/251‑1 (MQ=255)
                                                         ggCTGCGTAACGAAGGCTGTGAGGTATGCAGTTTTTCCGCCGCACGGGTAAAGTTCAGTGCCTGCGCCACTGCGACGAAATAGCGTAAATGCCGTAGTTCCATCACCTTCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTATGTCACATATTAGCAACCAACCAGCAACATTCTTATGGCACAAAAATAGAAGGTCAATACATCTTATCTTTCAGGATTAAAAAATATGACCACAccc                                                                                                                                                                                <  1:406030/250‑1 (MQ=255)
                                                             gctTAACGAAGGCTGTGAGGTATGCAGTTTTTTAGCCGCACGGGTAAAGTTCTATGCCTGCGCCACTGAGACGAAATAGCGTAAATGCCGTAGTTCCATTACCTTCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTTGAAATTAACTCTTTCACATATTAGCAACCAACCAGCAACATTCTTATGGCACAAACATAGAAGGACAATACATCTTATCTTACAGGATTAAAAAATATGACCACACACTc                                                                                                                                                                              >  2:269608/1‑248 (MQ=255)
                                                             gctTAACGAAGGCTGTGAGGTATGCAGTTTTTCCGCCGCACGGGTAAAGTTCAGTGCCTGCGCCACTGCGACGAAATAGCGTAAATGCCGTAGTTCCATCACCTTCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTATTTCACATATTAGCAACCAACCAGCAACATTCTTATGGCACAAAAATAGAAGGTCAATACATCTTATCTTTCAGGATTAAAAAATATGACCACACCCTCAg                                                                                                                                                                            <  1:409845/250‑1 (MQ=255)
                                                                   cGAAGGCTGTGAGGTATGCAGTTTTTCCGCCGCACGGGTAAAGTTCAGTGCCTGCGCCACTGCGACGAAAAAGCGTAAATGCCGTAGTTCCATCACCTTCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTATTTCACATATTAGCAACCAACCAGCAACATTCTTATGGCACAAAAATAGAAGGTCAATACATCTTATCTTTCAGGATTAAAAAATATGACCACACCCTCAGATTTGaa                                                                                                                                                                     <  1:508763/251‑1 (MQ=255)
                                                                        gCTGTGAGGTATGCAGTTTTTCCGCCGCACGGGTAAAGTTCAGTGCCTGCGCCACTGCGACGAAATAGCGTAAATGCCGTAGTTCCATCACCTTCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTATTTCACATATTAGCAACCAACCAGCAACATTCTTATGGCACAAAAATAGAAGGTCAATACATCTTATCTTTCAGGATTAAAAAATATGACCACACCCTCAGATGTGAACAttt                                                                                                                                                                <  1:167766/251‑1 (MQ=255)
                                                                         cTGTGTGGTATGCCGTTTTTCCGCCGCACGGGTAACGTTCAGAGCCTGCACAACTGCGACGAAAGAGCGTACATGCAGTAGTTCCATCACCTTCACCATGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTATTTCACATAGTAGAAACAAACCAGCGACATTATAATGGCACAAAGATAGAAGGACACTACAACTTATCTTTCAGGATTAAAAATTATGAcagtgccctacagtttgaccatgt                                                                                                                                                                >  2:482427/1‑227 (MQ=255)
                                                                                          tttCCGCCGCACGGGTAAAGTTCAGTGCCTGCGCCACTGCGACGAAATAGCGTAAATGCCGTAGTTCCATCACCTTCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTATTTCACATATTAGCAACCAACCAGCAACATTCTTATGGCACAAAAATAGAAGGTCAATACATCTTATCTTTCAGGATTAAAAAATATGACCACACCCTCAGATTTGAACATTTACCAACTGATTGACAcc                                                                                                                                               <  1:307742/250‑1 (MQ=255)
                                                                                            tCCGCCGCACGGGTAAAGTTCAGTGCCTGCGCCACTGCGACGAAATAGCGTAAATGCCGTAGTTCCATCACCTTCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTATTTCACATATTAGCAACCAACCAGCAACATTCTTATGGCACAAAAATAGAAGGTCAATACATCTTATCTTTCAGGATTAAAAAATATGACCACACCCTCAGATTTGAACATTTACCAACTGATTGACACCCa                                                                                                                                             <  1:317826/250‑1 (MQ=255)
                                                                                            tCCGCCGCACGGGTAAAGTTCAGTGCCTGCGCCACTGCGACGAAATAGCGTAAATGCCGTAGTTCCATCACCTTCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTATTTCACATATTAGCAACCAACCAGCAACATTCTTATGGCACAAAAATAGAAGGTCAATACATCTTATCTTTCAGGATTAAAAAATATGACCACACCCTCAGATTTGAACATTTACCAACTGATTGACACCCa                                                                                                                                             <  1:801626/250‑1 (MQ=255)
                                                                                                        gaaatgttcagagcaagcgcCACTGCTACGAAATAGCGTGTATGCCGTAGTTCCATCACCTTCCCCTTGTTATCGAAAAAACGGCTCAGCTGGTCGAAAATAACTATTTCACATATTAGCAACCAACCAGCAACCTTCTTATGGCCCCAAAATAGAAGGTCAATACATCTTATCTTTCAGGTTTAAAAAATAGGACCACAGCATCAGATTTGAACATTTACCAACAGATTGAGACCAAAAATGGTAGGGTc                                                                                                                                >  2:659750/17‑251 (MQ=255)
                                                                                                                     ccTGCGCCACTGCGACGAAATAGCGTAAATGCCGTAGTTCCAACACCTTCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTATTTCACATATTAGCAACCAACCAGCAACATTCTTATGGCACAAAAATAGAAGGTCAATACATCTTATCTTTCAGGATTAAAAAATATGACCACACCCTCAGATTTGAACATTTACCAACTGATTGACACCCAAAATGGTCGGGTCACTCCGCGTAttt                                                                                                                   <  2:689526/251‑1 (MQ=255)
                                                                                                                            cACTGCGACGAAATAGCGTAAATGCCGTAGTTCCATCACCTTCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTATTTCACATATTAGCAACCAACCAGCAACATTCTTATGGCACAAAAATAGAAGGTCAATACATCTTATCTTTCAGGATTAAAAAATATGACCACACCCTCAGATTTGAACATTTACCAACTGATTGACACCCAAAATGGTCGGGTCACTCCGCGTATTTATACCg                                                                                                             >  2:456772/1‑250 (MQ=255)
                                                                                                                                gCGACGAAATAGCGTAAATGCCGTAGTTCCATCACCTTCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTATTTCACATATTAGCAACCAACCAGCAACATTCTTATGGCACAAAAATAGAAGGTCAATACATCTTATCTTTCAGGATTAAAAAATATGACCACACCCTCAGATTTGAACATTTACCAACTGATTGACACCCAAAATGGTCGGGTCACTCCGCGGATTTATACCGAcc                                                                                                          <  1:669045/249‑1 (MQ=255)
                                                                                                                                 cgacgaAATAGCGTAAATGCCGTAGTTCCATCACCTTCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTATTTCACATATTAGCAACCAACCAGCAACATTCTTATGGCACAAAAATAGAAGGTCAATACATCTTATCTTTCAGGATTAAAAAATATGACCACACCCTCAGATTTGAACATTTACCAACTGATTGACACCCAAAATGGTCGGGTCACTCCGCGTATTTATACCGACCCgg                                                                                                       >  2:71457/1‑251 (MQ=255)
                                                                                                                                            cGTAAATGCCGTAGTTCCATCACCTTCCCCTCGTTATCGAAACAACGTCTCAGCTGGTAGAAATTAACTATTTCACATATTAGCAACCAACCAGCAACATTCTTATGGCACAAAAATAGAAGGTCAATACATCTTATCTTTCAGGATTAAAAAATATGACCACACCCTCAGATTTGAACATTTACCAACTGATTGACACCCAAAATGGTCGGGTCACTCCGCGTATTTATACCGACCCGGATCTTTACCa                                                                                             >  1:187659/1‑250 (MQ=255)
                                                                                                                                                     cGTAGTTCCATCACCTTCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTATTTCACATATTAGCAACCAACCAGCAACATTCTTATGGCACAAAAATAGAAGGTCAATACATCTTATCTTTCAGGATTAAAAAATATGACCACACCCTCAGATTTGAACATTTACCAACTGATTGACACCCAAAATGGTCGGGTCACTCCGCGTATTTATACAGACCCGGTTATTTACCAACTGGAGCtt                                                                                   >  1:85266/1‑251 (MQ=255)
                                                                                                                                                                             ttATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTATTTCACATATTAGCAACCAACCAGCAACATTCTTATGGCACAAAAATAGAAGGTCAATACATCTTATCTTTCAGGATTAAAAAATATGACCACACCCTCAGATTTGAACATTTACCAACTGATTGACACCCAAAATGGTCGGGTCACTCCGCGTATTTATACCGACCCGGATATTTACCAACTGGAGCTTGAGCGTATTTTCGGTCGTTGCTgg                                                           >  1:111750/1‑251 (MQ=255)
                                                                                                                                                                              tATCGAAAAAACGTCTCAGCTGGTAGCAATTAACTATTTCACATATTAGCAACCAACCAGCAACATTCTTATGGCACAAAAATAGAAGGTCAATACATCTTATCTTTCAGGATTAAAAAATATGACCACACCCTCAGATTTGAACATTTACCAACTGATTGACACCCAAAATGGTCGGGTCACTCCGCGTATTTATACCGACCCGGATATTTACCAACTGGAGCTTGAGCGTATTTTCGGTCGTTGCTGGc                                                          >  2:756814/1‑251 (MQ=255)
                                                                                                                                                                                   aaaaaaCGTCTCAGCTGGTAGAAATTAACTATTTCACATATTAGCAACCAACCAGCAACATTCTTATGGCACAAAAATAGAAGGTCAATACATCTTATCTTTCAGGATTAAAAAATATGACCACACCCTCAGATTTGAACATTTACCAACTGATTGACACCCAAAATGGTCGGGTCACTCCGCGTATTTATACCGACCCGGATATTTACCAACTGGAGCTTGAGCGTATTTTCGGTCGTTGCTGGCTAttt                                                     >  1:646794/1‑251 (MQ=255)
                                                                                                                                                                                    aaaaaCGTCTCATCTGGTAGAAATTAACTATTTCACATATTAGCAACCAACCAGCAACATTCTTATGGCACAAAAATAGAAGGTCAATACATCTTATCTTTCAGGATTAAAAAATATGACCTCACCCTCAGATTTGAACATTTACCAACTGATTGACACCCAAAATGGTCGGGTCACTCCGCGTATTTATACCGACCCGGATATTTTCCAACTGGAGGTTGAGCGTATTTTCGGTCGTTGCTGGCTATTtc                                                    <  2:775116/251‑1 (MQ=255)
                                                                                                                                                                                                gCTGGTAGAAATTAACTATTTCACATATTAGCAACCAACCAGCAACATTCTTATGGCACAAAAATAGAAGGTCAATACATCTTATCTTTCAGGATTAAAAAATATGACCACACCCTCAGATTTGAACATTTACCAACTGATTGACACCCAAAATGGTCGGGTCACTCCGCGTATTTATACCGACCCGGATATTTACCAACTGGAGCTTGAGCGTATTTTCGGTCGTTGCTGGCTATTTCTCGCCCACGaaa                                        <  1:350261/251‑1 (MQ=255)
                                                                                                                                                                                                  tGGTAGAAATTAACTATTTCACATATTAGCAACCAACCAGCAACATTCTTATGGCACAAAAATAGAAGGTCAATACATCTTATCTTTCAGGATTAAAAAATATGACCACACCCTCAGATTTGAACATTTACCAACTGATTGACACCCAAAATGGTCGGGTCACTCCGCGTATTTATACCGACCCGGat                                                                                                     <  2:624836/188‑1 (MQ=255)
                                                                                                                                                                                                  tGGTAGAAATTAACTATTTCACATATTAGCAACCAACCAGCAACATTCTTATGGCACAAAAATAGAAGGTCAATACATCTTATCTTTCAGGATTAAAAAATATGACCACACCCTCAGATTTGAACATTTACCAACTGATTGACACCCAAAATGGTCGGGTCACTCCGCGTATTTATACCGACCCGGat                                                                                                     >  1:624836/1‑188 (MQ=255)
                                                                                                                                                                                                    gTAGAAATTAACTATTTCACATATTAGCAACCATCCAACAACATTCTTATGGCACAAAAATAGAAGGTCAATACATCTTATCTTTCAGGATTAAAAAATATGACCACACCCTCAGATTTGAACATTTACAAACTGATTGACACCCAAAATGGTCGGGTCACTCCGCGTATTTATACCGACCCGGATATTTACCAACTGGAGCTTGAGCGTATTTTCGGTCGTTGCTGTCTATTTCTCGCCCACGAAAGcc                                     >  2:226312/1‑250 (MQ=255)
                                                                                                                                                                                                            tAACTATTTCACATATTAGCAACCAACCAGCAACATTCTTATGGCACAAAAATAGAAGGTCAATACATCTTATCTTTCAGGATTAAAAAATATGACCACACCCTCAGATTTGAACATTTACCAACTGATTGACACCCAAAATGGTCGGGTCACTCCGCGTATTTATACCGACCCGGATATTTACCAACTGGAGCTTGAGCGTATTTTCGGTCGTTGCTGGCTATTTCTCGCCCACGAAAGCCAGATCCCa                             >  1:758963/1‑250 (MQ=255)
                                                                                                                                                                                                                  tttCACATATTAGCAACCAACCAGCAACATTCTTATGGCACAAAAATAGAAGGACAATACATCTTATCTTTCAGGATTAAAAAATATGACCACACCCTCAGATTTGAACATTTACCAACTGATTGACACCCAAAATGGTCGGGTCACTCCGCGTATTTATACCGACCCGGATATTTACCCACTGTAGCTTGAGCGTATTTTCGGTCGTTGCTGGCTATTTCTCGCCCACGAAAGTCAGATCCCaaaaatcg                      >  1:159982/1‑247 (MQ=255)
                                                                                                                                                                                                                   ttCACATATTAGCAACCAAACAGCAACATTCTTATGGCACAAAAATAGAAGGTCAATACATCTTATCTTTCAGGATTAAAAAATATGACCACACCCTCAGATTTGAACATTTACCAACTGATTGACACCCAAAATGGTCGGGTCACTCCGCGTATTTATACCGACCCGGCTATTTACCAACTGGAGCTTGAGCGTATTTTCGGTCGTTGCTGGCTATTTCTCGCCCACGAAAGCCAGATCCCAAAACCCgg                     >  2:486806/1‑251 (MQ=255)
                                                                                                                                                                                                                    tCACATATTAGCAACCAACCAGCAACATTCTTATGGCACAAAAATAGAAGGTCAATACATCTTATCTTTCAGGATTAAAAAATATGACCACACCCTCAGATTTGAACATTTACCAACTGATTGACACCCAAAATGGTCGGGTCACTCCGCGTATTTATACCGACCCGGATATTTCCCGACTGGAGCTTGAGCGTATTTTCGGTCGTTGCt                                                             >  1:634604/1‑210 (MQ=255)
                                                                                                                                                                                                                    tCACATATTAGCAACAAACCAGCTACATTCTTAGTGCATAAACATAGAATGTCAATACATCTTATCTTTCAGGATTAAAAAATATGATCACACTCTCAGATTTGAACATTTACCAACTGATTGACTCCCAAATGGGTCGGGTCACTCCGCGTATTTATACCGACCCGGATATTTACCAACTGGAGCTTGAGCGTAGTTTCGGTCGGTGCt                                                             <  2:634604/210‑1 (MQ=255)
                                                                                                                                                                                                                     cacaTATTAGCAACCAACCAGCAACATTCTTATTGCACAAAAATAGACAGTCAATACATCTTAGCTTTCAGGATTAAAAAATATGACCACACCCTCATATTTGAACATTTACCAACTGACTGACACCCAAAATGGTCGGGTCACTCCGCGTATTTATACCGACCCGGATATTTACCAACTTGAGCTTGATCGTATTTTCGGCCGTTGGTGGCTATTTCGCGCCCACGAAAGCCAGATCCCAAACCCCGGt                    <  2:232872/250‑1 (MQ=255)
                                                                                                                                                                                                                       caTATTAGCAACCAACCAGCAACATTCTTATGGCACAAAAATAGAAGGTCAATACATCTTATCTTTCAGGATTAAAAAATATGACCACACCCTCAGATTTGAACATTTACCAACTGATTGACACCCAAAATGGTCGGGTCACTCCGCGTATTTATACCGACCCGGATATTTACCAACTGGAGCTTGAGCGTATTTTCGGTCGTTGCTGGCTATTTCTCGCCCACGAAAGCCAGATCCCAAAACCCGGTGAt                 >  1:676556/1‑251 (MQ=255)
                                                                                                                                                                                                                             aGCAACCAACCAGCAACATTCTTATGGCACAAAAATAGAAGGTCAATACATCTTATCTTTCAGGATTAAAAAATATGACCACACCCTCAGATTTGAACATTTACCAACTGATTGACACCCAAAATGGTCGGGTCACTCCGCGTATTTATACCGACCCGGATATTTACCAACTGGAGCTTGAGCGTATTTTCGGTCGTTGCTGGCTATTTCTCGCCCACGAAAGCCAGATCCCAAAACCCGGTGATTTCttt           <  2:184398/251‑1 (MQ=255)
                                                                                                                                                                                                                             aGCAACCAACCAGCAACATGCTTATGGCACAAAAATAGAAGGTCAATACATCTTATCTTTCAGGATTAAAAAATATGACCACACCCTCAGATTTGAACATTTACCAACTGATTGACACCCAAAATGGTCGGGTCACTCCGCGTATTTATACCGACCCGGATATTTACCAACTGGAGCTTGAGCGTATTTTCGGTCGTTGCTGGCTATTTCTCGCCCACGAAAGCCAGATCCCAAAACCCGGTGATTTCtt            <  1:542583/250‑1 (MQ=255)
                                                                                                                                                                                                                                aaccaaccAGCAACATTCTTATGGCACAAAAATAGAAGGTCAATACATCTTATCTTTCAGGATTAAAAAATATGACCACACCCTCAGATTTGAACATTTACCAACTGATTGACACCCAAAATGGTCGGGTCACTCCGCGTATTTATACCGACCCGGATATTTACCAACTGGAGCTTGAGCGTATTTTCGGTCGTTGCTGGCTATTTCTCGCCCACGAAAGCCAGATCCCAAAACCCGGTGATTTCTTTAac        >  2:742917/1‑251 (MQ=255)
                                                                                                                                                                                                                                aaccaaccAGCAACATTCTTATGGCACAAAAATAGAAGGTCAATACATCTTATCTTTCAGGATTAAAAAATATGACCACACCCTCAGATTTGAACATTTACCAACTGATTGACACCCAAAATGGTCGGGTCACTCCGCGTATTTATACCGACCCGGATATTTACCAACTGGAGCTTGAGCGTATTTTCGGTCGTTGCTGGCTATTTCTCGCCCACGAAAGCCAGATCCCAAAACCCGGTGATTTCTTTAac        <  1:371496/251‑1 (MQ=255)
                                                                                                                                                                                                                                       caGCAACATTCTTATGGCACAAAAATAGAAGGTCAATACATCTTATCTTTCAGGATTAAAAAATATGACCACACCCTCAGATTTGAACATTTACCAACTGATTGACACCCAAAATGGTCGGGTCACTCCGCGTATTTATACCGACCCGGATATTTACCAACTGGAGCTTGAGCGTATTTTCGGTCGTTGCTGGCTATTTCTCGCCCACGAAAGCCAGATCCCAAAACCCGGTGATTTCTTTAACACCTAc  >  1:129226/1‑250 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                
TGCGCTTATCCCTCACCAGTAACGGAACACCGACACAGTTTTCAAGATCGCGGATCTGGCTGCTTAACGAAGGCTGTGAGGTATGCAGTTTTTCCGCCGCACGGGTAAAGTTCAGTGCCTGCGCCACTGCGACGAAATAGCGTAAATGCCGTAGTTCCATCACCTTCCCCTTGTTATCGAAAAAACGTCTCAGCTGGTAGAAATTAACTATTTCACATATTAGCAACCAACCAGCAACATCCTTATGGCACAAAAATAGAAGGTCAATACATCTTATCTTTCAGGATTAAAAAATATGACCACACCCTCAGATTTGAACATTTACCAACTGATTGACACCCAAAATGGTCGGGTCACTCCGCGTATTTATACCGACCCGGATATTTACCAACTGGAGCTTGAGCGTATTTTCGGTCGTTGCTGGCTATTTCTCGCCCACGAAAGCCAGATCCCAAAACCCGGTGATTTCTTTAACACCTAC  >  NC_000913/2668737‑2669217

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