Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,672,639 A→C 6.9% 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.9% 333.4 / 15.9 101Q198P (CAG→CCG) ‡hcaDphenylpropionate dioxygenase, ferredoxin reductase subunit
Reads supporting (aligned to +/- strand):  ref base A (48/46);  new base C (3/4);  total (51/50)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.67e-01

Only 100 of 103 total aligned reads displayed.

ATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCAACTGGCTC  >  NC_000913/2672491‑2672785
                                                                                                                                                     |                                                                                                                                                  
aTTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACcagca                                                                                                                                                    <  3:82906/149‑1 (MQ=255)
  tATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAgg                                                                                                                                                  <  2:728972/149‑1 (MQ=255)
         gCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCCTGCTGgtgg                                                                                                                                           <  4:150828/149‑2 (MQ=255)
            ggAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCACACCAGCCTGCTGgtggcg                                                                                                                                         <  4:449278/149‑4 (MQ=255)
              aaCTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTCTACAGCGC‑CACCAGCAGGCTGGTGTGCGCa                                                                                                                                      <  1:514118/149‑1 (MQ=255)
               aCTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCAt                                                                                                                                     >  3:682172/1‑149 (MQ=255)
                cTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCAtt                                                                                                                                    >  1:347146/1‑149 (MQ=255)
                 tATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTc                                                                                                                                   >  1:801731/1‑149 (MQ=255)
                  aTTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCt                                                                                                                                  <  1:834810/149‑1 (MQ=255)
                   ttGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTg                                                                                                                                 >  3:794274/1‑149 (MQ=255)
                       tCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCa                                                                                                                             >  4:207688/1‑149 (MQ=255)
                        cTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATTCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCaa                                                                                                                            <  2:284472/149‑1 (MQ=255)
                         tGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCaat                                                                                                                           >  1:289297/1‑149 (MQ=255)
                         tGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCaat                                                                                                                           >  1:706636/1‑149 (MQ=255)
                         tGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGAGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGCTGTGCGCATTCTGCTCaat                                                                                                                           >  4:92801/1‑149 (MQ=255)
                           gAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCaataa                                                                                                                         >  2:420486/1‑149 (MQ=255)
                           gAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCaataa                                                                                                                         >  2:546741/1‑149 (MQ=255)
                           gAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCaataa                                                                                                                         >  3:16002/1‑149 (MQ=255)
                               tGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGcc                                                                                                                     >  4:697484/1‑149 (MQ=255)
                                 gCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCAt                                                                                                                   <  1:756490/149‑1 (MQ=255)
                                     ccAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGaa                                                                                                               >  2:796479/1‑149 (MQ=255)
                                     ccAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGaa                                                                                                               <  4:347879/149‑1 (MQ=255)
                                     ccAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGaa                                                                                                               <  2:4347/149‑1 (MQ=255)
                                           ccACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACAtgtg                                                                                                         >  1:184338/1‑149 (MQ=255)
                                            cACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTgg                                                                                                        >  4:191875/1‑149 (MQ=255)
                                            cACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTgg                                                                                                        >  4:643674/1‑149 (MQ=255)
                                              cgcagcgcagATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTc                                                                                                      <  4:794274/149‑1 (MQ=255)
                                                cagcgcagATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGa                                                                                                    <  1:483289/149‑1 (MQ=255)
                                                 agcgcagATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGAt                                                                                                   <  2:555731/149‑1 (MQ=255)
                                                 agcgcagATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGAt                                                                                                   <  3:86350/149‑1 (MQ=255)
                                                   cgcagATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATgg                                                                                                 <  2:289297/149‑1 (MQ=255)
                                                    gcagATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGt                                                                                                >  1:504317/1‑149 (MQ=255)
                                                     cagATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTg                                                                                               <  4:191779/149‑1 (MQ=255)
                                                        aTGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGaaa                                                                                            >  4:522976/1‑149 (MQ=255)
                                                          gTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑GACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGaaaaa                                                                                          <  4:559037/149‑1 (MQ=255)
                                                            aaGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGt                                                                                        >  1:374161/1‑149 (MQ=255)
                                                             aGGTGATAGTGATTGACCTGGCGGCAACCGTCATGGGCCGTCATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCCGCAGGCTGGTGTGCGAATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAGAAGTa                                                                                       <  4:156818/149‑1 (MQ=255)
                                                             aGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGc                                                                                                                                >  2:391918/1‑108 (MQ=255)
                                                             aGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGc                                                                                                                                <  1:391918/108‑1 (MQ=255)
                                                             aGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTa                                                                                       <  2:801731/149‑1 (MQ=255)
                                                             aGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCATGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTa                                                                                       <  2:252656/149‑1 (MQ=255)
                                                               gTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGa                                                                                     >  1:829273/1‑149 (MQ=255)
                                                                 gACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAAc                                                                                   <  3:196368/149‑1 (MQ=255)
                                                                 gACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAAc                                                                                   <  3:161480/149‑1 (MQ=255)
                                                                    aGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGa                                                                                >  2:106163/1‑149 (MQ=255)
                                                                       gATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGc                                                                             <  3:655953/149‑1 (MQ=255)
                                                                         ttGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTg                                                                           >  4:669865/1‑149 (MQ=255)
                                                                             aCTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCaaa                                                                       >  1:305566/1‑149 (MQ=255)
                                                                             aCTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCaaa                                                                       <  4:447758/149‑1 (MQ=255)
                                                                               tGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGt                                                                     <  1:546741/149‑1 (MQ=255)
                                                                                ggcggcAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTg                                                                    <  3:218204/149‑1 (MQ=255)
                                                                                 gcggcAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTgg                                                                   >  1:533615/1‑149 (MQ=255)
                                                                                  cggcAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTggg                                                                  <  1:219554/149‑1 (MQ=255)
                                                                                      aaCCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGaaa                                                              >  1:108500/1‑149 (MQ=255)
                                                                                      aaCCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGaaa                                                              >  3:141252/1‑149 (MQ=255)
                                                                                         cGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGaaa                                                                                            <  2:606289/116‑1 (MQ=255)
                                                                                         cGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGaaa                                                                                            >  1:606289/1‑116 (MQ=255)
                                                                                             atcggCCGGAATGCACCACCGCGCGAGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTTGTGTGCGCCTTCTGCTCAGTAATGCCGTGGAACATGTGGGCGATGGTGAAAAGGTAGAACGGAGGCTGCTAGGAGGTGAAAGGCttta                                                       >  4:685274/4‑147 (MQ=255)
                                                                                             aTGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCa                                                       >  3:813771/1‑149 (MQ=255)
                                                                                             aTGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCa                                                       <  1:420486/149‑1 (MQ=255)
                                                                                             aTGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCa                                                       <  1:443911/149‑1 (MQ=255)
                                                                                             aTGGGCCGTAATGCACCACCGCCCGTGCAACGCTACCTTCTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGTCGCTGCAAAGTGGGGAAACGCTTCa                                                       <  2:675509/149‑1 (MQ=255)
                                                                                               gggCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAgg                                                     >  2:406847/1‑149 (MQ=255)
                                                                                                 gCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCt                                                   <  4:525569/149‑1 (MQ=255)
                                                                                                 gCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCt                                                   >  3:249650/1‑149 (MQ=255)
                                                                                                    gTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGAt                                                >  3:8125/1‑149 (MQ=255)
                                                                                                        tGCACCACCGGCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATgtgg                                            <  1:538594/149‑1 (MQ=255)
                                                                                                        tGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATgtgg                                            <  2:41922/149‑1 (MQ=255)
                                                                                                        tGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATgtgg                                            <  2:376226/149‑1 (MQ=255)
                                                                                                         gCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATgtggt                                           <  4:229805/149‑1 (MQ=255)
                                                                                                           accaccGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGa                                         >  3:405194/1‑149 (MQ=255)
                                                                                                            ccaccGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGAt                                        <  1:796479/149‑1 (MQ=255)
                                                                                                             caccGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGAtt                                       >  1:433132/1‑149 (MQ=255)
                                                                                                              accGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGAttt                                      >  1:146234/1‑149 (MQ=255)
                                                                                                              accGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGAttt                                      >  1:560934/1‑149 (MQ=255)
                                                                                                                cGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTAc                                    <  3:669865/149‑1 (MQ=255)
                                                                                                                 gCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACg                                   <  4:405194/149‑1 (MQ=255)
                                                                                                                     gTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTAt                               <  1:58413/149‑1 (MQ=255)
                                                                                                                      tGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTAtt                              >  1:670879/1‑149 (MQ=255)
                                                                                                                       gCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTg                             <  2:533615/149‑1 (MQ=255)
                                                                                                                         aaCGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACg                                   >  2:368420/1‑141 (MQ=255)
                                                                                                                         aaCGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACg                                   <  1:368420/141‑1 (MQ=255)
                                                                                                                         aaCGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGt                           >  2:829491/1‑149 (MQ=255)
                                                                                                                         aaCGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGt                           >  3:800376/1‑149 (MQ=255)
                                                                                                                          aCGCTATCTTTTGCAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTAACGGTATTGGTa                          <  2:132714/149‑1 (MQ=255)
                                                                                                                                tCTTTTACAGCGC‑CACCAGCAGGCCGGTGTGCGCATGCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAgcg                    >  4:228992/1‑149 (MQ=255)
                                                                                                                                  ttttACAGCGC‑CACCAGCAGGCTGGTGTGTGAATTCTGCTCAATAATGCCATTGAAAATGTCGACGATTGTGAAAAAGTAGAACTGACACCGCAAAGGGGCGCAACGCTTAAAGCTGATGTGGTGCTTTAAGGCATTCGTATCAGCGcc                  >  1:105887/1‑149 (MQ=255)
                                                                                                                                         gcgc‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGc           >  2:671212/1‑149 (MQ=255)
                                                                                                                                         gcgc‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGc           >  1:227021/1‑149 (MQ=255)
                                                                                                                                           gcacaccaGCCTGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTg                                                                    <  3:47836/88‑1 (MQ=255)
                                                                                                                                           gcacaccaGCCTGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTg                                                                    >  4:47836/4‑91 (MQ=255)
                                                                                                                                           gcacaccaGCCTGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGa                                         <  1:383451/115‑1 (MQ=255)
                                                                                                                                           gcacaccaGCCTGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGa                                         >  2:383451/4‑118 (MQ=255)
                                                                                                                                             acaccaGCCTGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCAAc        >  2:261995/2‑149 (MQ=255)
                                                                                                                                              caccaGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCAACt       >  3:507916/1‑149 (MQ=255)
                                                                                                                                               accaGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGc                                                    <  3:256771/103‑1 (MQ=255)
                                                                                                                                               accaGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGc                                                    >  4:256771/1‑103 (MQ=255)
                                                                                                                                                ccaGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCAACTgg     <  4:698970/149‑1 (MQ=255)
                                                                                                                                                ccaGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGCTCGATGGTGAAAAAGTAGAACTGACCCTGCAAAGCGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCAACTgg     >  4:620868/1‑149 (MQ=255)
                                                                                                                                                   gcagGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCAACTGGCTc  <  3:272455/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
ATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCAACTGGCTC  >  NC_000913/2672491‑2672785

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