Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,672,640 G→T 6.0% 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 6.0% 588.2 / 19.4 166Q198H (CAG→CAT) ‡hcaDphenylpropionate dioxygenase, ferredoxin reductase subunit
Reads supporting (aligned to +/- strand):  ref base G (79/77);  new base T (6/4);  total (85/81)
Fisher's exact test for biased strand distribution p-value = 7.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.66e-01

Only 100 of 170 total aligned reads displayed.

GGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCAACTGGCTCGCG  >  NC_000913/2672497‑2672788
                                                                                                                                                |                                                                                                                                                    
ggTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTgg                                                                                                                                                 >  3:790756/1‑149 (MQ=255)
     cGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTgc                                                                                                                                            <  3:826302/149‑1 (MQ=255)
      ggAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTgcg                                                                                                                                           >  2:48120/1‑149 (MQ=255)
       gAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCACACCAGCCTGCTGgtggcgc                                                                                                                                           <  2:256469/149‑5 (MQ=255)
           tATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTc                                                                                                                                      >  1:358762/1‑149 (MQ=255)
             ttGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTgc                                                                                                                                            >  1:1626471/1‑141 (MQ=255)
              tGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGc                                                                                                                                   >  2:1629111/1‑149 (MQ=255)
              tGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGc                                                                                                                                   <  4:790756/149‑1 (MQ=255)
                 tCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCa                                                                                                                                >  1:1603887/1‑149 (MQ=255)
                     gAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGATCaataa                                                                                                                            <  2:919123/149‑1 (MQ=255)
                      aaCTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCaataat                                                                                                                           <  3:1033460/149‑1 (MQ=255)
                         tGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGcc                                                                                                                        <  1:48120/149‑1 (MQ=255)
                          ggCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCa                                                                                                                       >  2:1291031/1‑149 (MQ=255)
                           gCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCAt                                                                                                                      <  3:326572/149‑1 (MQ=255)
                           gCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCGCAATAATGCCAt                                                                                                                      <  2:586355/149‑1 (MQ=255)
                               ccAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGaa                                                                                                                  >  3:1542617/1‑149 (MQ=255)
                               ccAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGaa                                                                                                                  <  2:1537973/149‑1 (MQ=255)
                                 aGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAAc                                                                                                                 >  1:439408/1‑148 (MQ=255)
                                       aCGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGACAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGt                                                                                                          >  2:427545/1‑149 (MQ=255)
                                            gcgcagATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATg                                                                                                     <  1:1669847/149‑1 (MQ=255)
                                            gcgcagATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATg                                                                                                     <  4:711071/149‑1 (MQ=255)
                                            gcgcagATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATg                                                                                                     <  2:1168212/149‑1 (MQ=255)
                                            gcgcagATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATg                                                                                                     <  3:1676149/149‑1 (MQ=255)
                                              gcagATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGt                                                                                                   <  3:915661/149‑1 (MQ=255)
                                              gcagATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGt                                                                                                   >  1:1463973/1‑149 (MQ=255)
                                                   tGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGt                                                                                                   <  1:1329776/144‑1 (MQ=255)
                                                   tGTAAGGTGACAGTGATTGAACTCGCGGCAACCGTCATGGGCCGTAATGCACCCACGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTATGCTCAATAATGCCATTGAACATGTGGTCGCTGGTGaaaa                                                                                              >  3:1447786/1‑149 (MQ=255)
                                                    gTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGaaaaa                                                                                             <  1:1274613/149‑1 (MQ=255)
                                                    gTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGGGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGACGGTGaaaaa                                                                                             <  2:1653606/149‑1 (MQ=255)
                                                     tAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAg                                                                                            >  1:1694951/1‑149 (MQ=255)
                                                          tGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGaa                                                                                       >  2:1710981/1‑149 (MQ=255)
                                                          tGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGaa                                                                                       <  2:1099800/149‑1 (MQ=255)
                                                          tGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGaa                                                                                       <  3:1671868/149‑1 (MQ=255)
                                                              aGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCGATCTTTTACAGCGC‑CACCAGCAGGCTTGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGa                                                                                   >  2:985574/1‑149 (MQ=255)
                                                               gTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGAc                                                                                  <  4:1542617/149‑1 (MQ=255)
                                                                tGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACg                                                                                 >  4:981105/1‑149 (MQ=255)
                                                                 gATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGc                                                                                <  3:317373/149‑1 (MQ=255)
                                                                  aTTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACAt                                                                                                               <  2:226363/117‑1 (MQ=255)
                                                                  aTTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCt                                                                               >  2:1044838/1‑149 (MQ=255)
                                                                  aTTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCt                                                                               <  1:855705/149‑1 (MQ=255)
                                                                  aTTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATGGAACAt                                                                                                               >  1:226363/1‑117 (MQ=255)
                                                                   ttGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCaataa                                                                                                                            >  3:1317689/1‑103 (MQ=255)
                                                                   ttGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCaataa                                                                                                                            <  4:1317689/103‑1 (MQ=255)
                                                                   ttGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTg                                                                              >  2:1572052/1‑149 (MQ=255)
                                                                   ttGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTg                                                                              >  2:910193/1‑149 (MQ=255)
                                                                      aaCTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCaa                                                                           >  3:549075/1‑149 (MQ=255)
                                                                       aCTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAGAGTAGAACTGACGCTGCaaa                                                                          <  1:186356/149‑1 (MQ=255)
                                                                           gcggcAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTgg                                                                      <  1:1044838/149‑1 (MQ=255)
                                                                            cggcAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTggg                                                                     <  3:1638983/149‑1 (MQ=255)
                                                                             ggcAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTgggg                                                                    >  4:1067500/1‑149 (MQ=255)
                                                                                aaCCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCCTTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGaaa                                                                 >  4:28439/1‑149 (MQ=255)
                                                                                 aCCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAAc                                                                >  1:408628/1‑149 (MQ=255)
                                                                                 aCCGTCACGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAAc                                                                <  1:26298/149‑1 (MQ=255)
                                                                                   cGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGaa                                                                  <  4:748630/145‑1 (MQ=255)
                                                                                   cGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGaa                                                                  >  3:748630/1‑145 (MQ=255)
                                                                                   cGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGc                                                              <  1:294708/149‑1 (MQ=255)
                                                                                   cGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGc                                                              >  4:934542/1‑149 (MQ=255)
                                                                                       aTGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACTTGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCa                                                          <  1:1104126/149‑1 (MQ=255)
                                                                                         gggCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAgg                                                        >  1:1141857/1‑149 (MQ=255)
                                                                                          ggCCGTAATGCACCACCGCCCGTGCAACGGTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGc                                                       >  3:1007027/1‑149 (MQ=255)
                                                                                          ggCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGc                                                       >  1:1585287/1‑149 (MQ=255)
                                                                                             cGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGa                                                    <  4:549075/149‑1 (MQ=255)
                                                                                             cGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGa                                                    <  1:879852/149‑1 (MQ=255)
                                                                                             cGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGa                                                    >  2:1316736/1‑149 (MQ=255)
                                                                                                aaTGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGAtgt                                                 >  3:257473/1‑149 (MQ=255)
                                                                                                   gCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATgtggt                                              <  1:529937/149‑1 (MQ=255)
                                                                                                   gCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATgtggt                                              <  4:620686/149‑1 (MQ=255)
                                                                                                   gCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATgtggt                                              <  2:1592518/149‑1 (MQ=255)
                                                                                                   gCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATgtggt                                              <  1:868920/149‑1 (MQ=255)
                                                                                                     accaccGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGa                                            >  2:1613704/1‑149 (MQ=255)
                                                                                                        accGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGAttt                                         >  2:887783/1‑149 (MQ=255)
                                                                                                             ccGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGc                                                                                                                                                    >  3:1560065/1‑37 (MQ=255)
                                                                                                             ccGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGc                                                                                                                                                    <  4:1560065/37‑1 (MQ=255)
                                                                                                              cGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTa                                   <  2:358762/149‑1 (MQ=255)
                                                                                                              cGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTa                                   <  4:1642185/149‑1 (MQ=255)
                                                                                                               gTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTAt                                  <  4:252954/149‑1 (MQ=255)
                                                                                                                tGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTAtt                                 <  2:861469/149‑1 (MQ=255)
                                                                                                                tGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTAtt                                 <  4:385653/149‑1 (MQ=255)
                                                                                                                  cAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTgg                               <  2:774707/149‑1 (MQ=255)
                                                                                                                   aaCGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGt                              >  1:88389/1‑149 (MQ=255)
                                                                                                                     cGCTATCTTTTTCAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTAt                            <  1:985574/149‑1 (MQ=255)
                                                                                                                     cGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATg                                                                                                                          <  3:730411/55‑1 (MQ=255)
                                                                                                                     cGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTAt                            >  4:629610/1‑149 (MQ=255)
                                                                                                                       cTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCa                          >  4:1143409/1‑149 (MQ=255)
                                                                                                                         aTCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAgc                        <  3:629610/149‑1 (MQ=255)
                                                                                                                          tCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAgcg                       >  1:623245/1‑149 (MQ=255)
                                                                                                                            ttttACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGcc                     >  2:11035/1‑149 (MQ=255)
                                                                                                                                 cAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGa                <  1:385481/149‑1 (MQ=255)
                                                                                                                                   gcgc‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGc              >  1:1191192/1‑149 (MQ=255)
                                                                                                                                     gcacaccaGCCTGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCa                                                          >  3:496154/4‑104 (MQ=255)
                                                                                                                                     gcacaccaGCCTGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCa             >  4:1474700/4‑149 (MQ=255)
                                                                                                                                     gcacaccaGCCTGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCa             >  4:134130/4‑149 (MQ=255)
                                                                                                                                     gcacaccaGCCTGCTGGTG‑GCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGAtgtg                                                <  3:918786/110‑1 (MQ=255)
                                                                                                                                      cacaccaGCCTGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTAtt                                 <  4:837946/126‑1 (MQ=255)
                                                                                                                                        caccaGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGccagt          <  4:72497/149‑5 (MQ=255)
                                                                                                                                        caccaGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCAACt          >  2:144261/1‑149 (MQ=255)
                                                                                                                                          ccaGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCAACTgg        <  2:408628/149‑1 (MQ=255)
                                                                                                                                           cagcagGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCAACTGGc       >  1:740592/1‑149 (MQ=255)
                                                                                                                                            agcagGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATgtgg                                               <  4:356232/108‑1 (MQ=255)
                                                                                                                                            agcagGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATgtgg                                               >  3:356232/1‑108 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCAACTGGCTCGCG  >  NC_000913/2672497‑2672788

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