Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,672,640 G→T 5.3% Q198H (CAG→CAT) ‡ hcaD → phenylpropionate dioxygenase, ferredoxin reductase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,672,6400GT 5.3% 340.3 / 7.6 94Q198H (CAG→CAT) ‡hcaDphenylpropionate dioxygenase, ferredoxin reductase subunit
Reads supporting (aligned to +/- strand):  ref base G (48/41);  new base T (3/2);  total (51/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.08e-01

TTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCAACTGGC  >  NC_000913/2672492‑2672783
                                                                                                                                                    |                                                                                                                                               
ttATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACcagcag                                                                                                                                                 >  3:651319/1‑149 (MQ=255)
  aTCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGc                                                                                                                                               >  1:921656/1‑149 (MQ=255)
   tCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCt                                                                                                                                              >  3:678635/1‑149 (MQ=255)
      gTGCCGGAACGATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGgt                                                                                                                                           <  4:981156/149‑1 (MQ=255)
         ccGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGtgtg                                                                                                                                        <  1:851253/149‑1 (MQ=255)
         ccGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGtgtg                                                                                                                                        >  1:653144/1‑149 (MQ=255)
           ggAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTgcg                                                                                                                                      >  3:699860/1‑149 (MQ=255)
                tATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTc                                                                                                                                 >  1:349571/1‑149 (MQ=255)
                      tCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCa                                                                                                                           >  4:1021260/1‑149 (MQ=255)
                      tCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCa                                                                                                                           >  3:897660/1‑149 (MQ=255)
                      gcTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCa                                                                                                                           >  1:881963/2‑149 (MQ=255)
                          gAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCaataa                                                                                                                       >  3:687207/1‑149 (MQ=255)
                           aaCTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCaataat                                                                                                                      <  4:897660/149‑1 (MQ=255)
                           aaCTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTAATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCaataat                                                                                                                      >  1:571917/1‑149 (MQ=255)
                                gCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCAt                                                                                                                 >  2:162339/1‑149 (MQ=255)
                                  tGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTg                                                                                                               >  4:907074/1‑149 (MQ=255)
                                   gCCAGCGCCACACAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGa                                                                                                              >  3:652640/1‑149 (MQ=255)
                                    ccAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGaa                                                                                                             >  4:646691/1‑149 (MQ=255)
                                      aGCGCCACGCAGCGCCGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACa                                                                                                           >  4:641190/1‑149 (MQ=255)
                                      aGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACa                                                                                                           <  3:1072047/149‑1 (MQ=255)
                                         gcCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACAtgt                                                                                                        <  4:699860/149‑1 (MQ=255)
                                           cACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTgg                                                                                                      <  2:349571/149‑1 (MQ=255)
                                            aCGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGt                                                                                                     <  3:1057646/149‑1 (MQ=255)
                                                 gcgccgATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATg                                                                                                <  3:641190/149‑1 (MQ=255)
                                                   gcagATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGt                                                                                              >  2:539372/1‑149 (MQ=255)
                                                        gatgaGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGaaaa                                                                                         <  4:1112536/145‑1 (MQ=255)
                                                         gTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCaa                                                                                                                          >  4:771436/1‑115 (MQ=255)
                                                         gTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCaa                                                                                                                          <  3:771436/115‑1 (MQ=255)
                                                          tAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAg                                                                                       <  2:809414/149‑1 (MQ=255)
                                                           aaGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGt                                                                                      >  3:991296/1‑149 (MQ=255)
                                                           aaGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGt                                                                                      >  2:362687/1‑149 (MQ=255)
                                                             ggTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATg                                                                                                <  4:523847/137‑1 (MQ=255)
                                                             ggTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATg                                                                                                >  3:523847/1‑137 (MQ=255)
                                                               tGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAATGTAGaa                                                                                  <  2:816019/149‑1 (MQ=255)
                                                                     tGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACg                                                                            >  1:733064/1‑149 (MQ=255)
                                                                        ttGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTg                                                                         >  3:410698/1‑149 (MQ=255)
                                                                                gcggcAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTgg                                                                 <  1:593080/149‑1 (MQ=255)
                                                                                gcggcAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTgg                                                                 >  1:575695/1‑149 (MQ=255)
                                                                                gcggcAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTgg                                                                 >  1:197081/1‑149 (MQ=255)
                                                                                 cggcAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTggg                                                                <  2:571917/149‑1 (MQ=255)
                                                                                  ggcAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTgggg                                                               >  1:515653/1‑149 (MQ=255)
                                                                                  ggcAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTgggg                                                               >  3:331367/1‑149 (MQ=255)
                                                                                  ggcAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTgggg                                                               >  4:924406/1‑149 (MQ=255)
                                                                                  ggcAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTgggg                                                               <  1:162339/149‑1 (MQ=255)
                                                                                     aaCCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGaaa                                                            <  4:678635/149‑1 (MQ=255)
                                                                                     aaCCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGaaa                                                            >  3:21420/1‑149 (MQ=255)
                                                                                      aCCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAAc                                                           <  3:924406/149‑1 (MQ=255)
                                                                                      aCCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAAc                                                           <  2:733064/149‑1 (MQ=255)
                                                                                         gTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCt                                                        >  2:868147/1‑149 (MQ=255)
                                                                                          tCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGt                                                                                              <  1:1142609/110‑1 (MQ=255)
                                                                                          tCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGt                                                                                              >  2:1142609/1‑110 (MQ=255)
                                                                                           cATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTc                                                      <  4:991296/149‑1 (MQ=255)
                                                                                             tGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAg                                                    <  4:331367/149‑1 (MQ=255)
                                                                                             tGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAg                                                    >  3:910077/1‑149 (MQ=255)
                                                                                                gCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGACACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCt                                                 >  3:833450/1‑149 (MQ=255)
                                                                                                   gTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACAtgt                                                                                                        <  3:512604/91‑1 (MQ=255)
                                                                                                   gTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACAtgt                                                                                                        >  4:512604/1‑91 (MQ=255)
                                                                                                   gTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGTGGAAACGCTTCAGGCTGAt                                              <  3:907074/149‑1 (MQ=255)
                                                                                                     aaTGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTg                                                                         <  4:822556/120‑1 (MQ=255)
                                                                                                     aaTGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTg                                                                         >  3:822556/1‑120 (MQ=255)
                                                                                                         caccacCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATgtggtg                                        >  2:969220/1‑149 (MQ=255)
                                                                                                          accaccGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTTGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGa                                       >  2:918452/1‑149 (MQ=255)
                                                                                                                gCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGGCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACg                                 <  4:652640/149‑1 (MQ=255)
                                                                                                                gCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATgtggtg                                        >  3:356665/1‑142 (MQ=255)
                                                                                                                gCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATgtggtg                                        <  4:356665/142‑1 (MQ=255)
                                                                                                                 cccGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCAtt                                                                                                                <  2:345994/69‑1 (MQ=255)
                                                                                                                 cccGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCAtt                                                                                                                >  1:345994/1‑69 (MQ=255)
                                                                                                                    gTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGAt                                      >  4:637453/1‑140 (MQ=255)
                                                                                                                    gTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGAt                                      <  3:637453/140‑1 (MQ=255)
                                                                                                                      gCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTg                           <  4:651319/149‑1 (MQ=255)
                                                                                                                       cAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATgtgg                                          <  2:222643/133‑1 (MQ=255)
                                                                                                                       cAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATgtgg                                          >  1:222643/1‑133 (MQ=255)
                                                                                                                        aaCGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGt                         >  3:709480/1‑149 (MQ=255)
                                                                                                                        aaCGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGt                         >  3:689887/1‑149 (MQ=255)
                                                                                                                         aCGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATgtgg                                          <  4:434494/131‑1 (MQ=255)
                                                                                                                         aCGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATgtgg                                          >  3:434494/1‑131 (MQ=255)
                                                                                                                            cTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCa                     <  2:653144/149‑1 (MQ=255)
                                                                                                                             tATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGaa                                                                                  <  1:309684/87‑1 (MQ=255)
                                                                                                                             tATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGaa                                                                                  >  2:309684/1‑87 (MQ=255)
                                                                                                                               tCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTg                           <  1:175120/140‑1 (MQ=255)
                                                                                                                               tCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTg                           >  2:175120/1‑140 (MQ=255)
                                                                                                                               tCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAgcg                  >  3:997361/1‑149 (MQ=255)
                                                                                                                                  tttACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTAtt                            >  2:135597/1‑136 (MQ=255)
                                                                                                                                  tttACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTAtt                            <  1:135597/136‑1 (MQ=255)
                                                                                                                                      cAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGa           <  3:646691/149‑1 (MQ=255)
                                                                                                                                      cAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGa           <  4:997361/149‑1 (MQ=255)
                                                                                                                                       aGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGc                                                         >  1:712892/1‑102 (MQ=255)
                                                                                                                                       aGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGc                                                         <  2:712892/102‑1 (MQ=255)
                                                                                                                                        gcgcCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGc         <  1:539372/149‑1 (MQ=255)
                                                                                                                                         gcacaccaGCCTGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGt                         <  2:337266/129‑1 (MQ=255)
                                                                                                                                         gcacaccaGCCTGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGt                         >  1:337266/4‑132 (MQ=255)
                                                                                                                                         gcacaccaGCCTGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCa        >  1:1019143/4‑149 (MQ=255)
                                                                                                                                          gcCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCaa       <  1:570443/149‑1 (MQ=255)
                                                                                                                                           acaccaGCCTGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCAAc      >  2:722838/2‑149 (MQ=255)
                                                                                                                                            caccaGCCTGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCAACt     <  2:1019143/149‑1 (MQ=255)
                                                                                                                                               cagcagGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCAACTGGc  >  4:570637/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
TTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCAACTGGC  >  NC_000913/2672492‑2672783

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