Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,672,639 A→C 6.0% Q198P (CAG→CCG) ‡ 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,6390AC 6.0% 578.4 / 26.4 166Q198P (CAG→CCG) ‡hcaDphenylpropionate dioxygenase, ferredoxin reductase subunit
Reads supporting (aligned to +/- strand):  ref base A (79/77);  new base C (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 = 1.07e-01

Only 100 of 170 total aligned reads displayed.

ATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCAACTGGC  >  NC_000913/2672491‑2672783
                                                                                                                                                     |                                                                                                                                                
aTTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACcagca                                                                                                                                                  <  1:32221/149‑1 (MQ=255)
      ggTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTgg                                                                                                                                            >  3:790756/1‑149 (MQ=255)
          ccGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGtgtg                                                                                                                                        <  1:1609777/149‑1 (MQ=255)
           cGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTgc                                                                                                                                       <  3:826302/149‑1 (MQ=255)
            ggAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTgcg                                                                                                                                      >  1:200296/1‑149 (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)
                    tGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGc                                                                                                                              <  4:790756/149‑1 (MQ=255)
                    tGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGc                                                                                                                              >  2:1629111/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                                                                                                <  3:1676149/149‑1 (MQ=255)
                                                  gcgcagATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATg                                                                                                <  2:1168212/149‑1 (MQ=255)
                                                  gcgcagATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATg                                                                                                <  4:711071/149‑1 (MQ=255)
                                                    gcagATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGt                                                                                              <  3:915661/149‑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                                                                                  <  3:1671868/149‑1 (MQ=255)
                                                                tGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGaa                                                                                  >  2:1710981/1‑149 (MQ=255)
                                                                tGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGaa                                                                                  <  2:1099800/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)
                                                                         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)
                                                                                 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                                                             >  3:748630/1‑145 (MQ=255)
                                                                                         cGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGaa                                                             <  4:748630/145‑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)
                                                                                               gggCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAgg                                                   >  1:1592518/1‑149 (MQ=255)
                                                                                               gggCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAgg                                                   >  1:1668798/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                                         <  1:868920/149‑1 (MQ=255)
                                                                                                         gCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATgtggt                                         <  2:1592518/149‑1 (MQ=255)
                                                                                                           accaccGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGa                                       >  2:1613704/1‑149 (MQ=255)
                                                                                                              accGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGAttt                                    >  2:887783/1‑149 (MQ=255)
                                                                                                                cGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTAc                                  <  1:1291031/149‑1 (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)
                                                                                                                         aaCGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGt                         >  1:1351569/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)
                                                                                                                            gCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATc                      <  1:179115/149‑1 (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)
                                                                                                                                           gcacaccaGCCTGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCa                                                     >  3:496154/4‑104 (MQ=255)
                                                                                                                                           gcacaccaGCCTGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCa        >  4:1474700/4‑149 (MQ=255)
                                                                                                                                           gcacaccaGCCTGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCa        >  4:134130/4‑149 (MQ=255)
                                                                                                                                           gcacaccaGCCTGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCa        >  1:279057/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)
                                                                                                                                                     |                                                                                                                                                
ATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCAACTGGC  >  NC_000913/2672491‑2672783

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