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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,672,6400GT 5.8% 363.2 / 10.2 104Q198H (CAG→CAT) ‡hcaDphenylpropionate dioxygenase, ferredoxin reductase subunit
Reads supporting (aligned to +/- strand):  ref base G (46/52);  new base T (2/4);  total (48/56)
Fisher's exact test for biased strand distribution p-value = 6.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.66e-01

Only 100 of 106 total aligned reads displayed.

GTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCAACTGGCTCGC  >  NC_000913/2672498‑2672787
                                                                                                                                              |                                                                                                                                                   
gTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTgg                                                                                                                                                >  3:854153/1‑148 (MQ=255)
gTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGgt                                                                                                                                               <  1:970235/149‑1 (MQ=255)
  gCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGgtgt                                                                                                                                             <  1:868174/149‑1 (MQ=255)
  gCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGgtgt                                                                                                                                             <  1:616842/149‑1 (MQ=255)
  gCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGgtg                                                                                                                                              <  4:599473/148‑1 (MQ=255)
   ccGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGtgtg                                                                                                                                            >  1:572859/1‑149 (MQ=255)
     ggAACTATTGGTCTGGAACTGGCTGCCAGCGCCATGCTGCGCAGATGTGAGGTGACAGTGATTGGACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTgcg                                                                                                                                          >  1:541995/1‑149 (MQ=255)
         cTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCAtt                                                                                                                                      <  1:612404/149‑1 (MQ=255)
          tATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCATCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTc                                                                                                                                     >  1:186145/1‑149 (MQ=255)
          tATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTc                                                                                                                                     >  2:389442/1‑149 (MQ=255)
           aTTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTc                                                                                                                                     >  4:39625/1‑148 (MQ=255)
              ggTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCt                                                                                                                                 <  2:248685/149‑1 (MQ=255)
               gTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCt                                                                                                                                 <  4:447078/148‑1 (MQ=255)
                 cTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCa                                                                                                                               <  3:999929/148‑1 (MQ=255)
                 cTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCa                                                                                                                               <  4:539703/148‑1 (MQ=255)
                  tGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCaat                                                                                                                             >  1:418282/1‑149 (MQ=255)
                    gagcTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCaataa                                                                                                                           <  2:394717/146‑1 (MQ=255)
                      aCTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCaataat                                                                                                                          >  4:1039279/1‑148 (MQ=255)
                        tGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGcc                                                                                                                       <  1:328588/149‑1 (MQ=255)
                         ggCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGcc                                                                                                                       >  3:761813/1‑148 (MQ=255)
                           cTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCAt                                                                                                                     <  4:593916/148‑1 (MQ=255)
                            tGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCAtt                                                                                                                    >  4:735958/1‑148 (MQ=255)
                             gCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTg                                                                                                                   >  4:115827/1‑148 (MQ=255)
                             gCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGa                                                                                                                  <  1:76041/149‑1 (MQ=255)
                               cAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGaa                                                                                                                 >  3:1041540/1‑148 (MQ=255)
                                aGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCCTGCTGgtgg                                                                                                                                             >  2:143490/1‑118 (MQ=255)
                                aGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCCTGCTGgtag                                                                                                                                             <  1:143490/119‑3 (MQ=255)
                                aGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAAc                                                                                                                <  3:300340/148‑1 (MQ=255)
                                aGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACa                                                                                                               <  1:12122/149‑1 (MQ=255)
                                   gcCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACAtg                                                                                                             <  4:828746/148‑1 (MQ=255)
                                    ccACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACAtgtg                                                                                                           >  1:754968/1‑149 (MQ=255)
                                       cgcagcgcagATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTc                                                                                                        <  2:145454/149‑1 (MQ=255)
                                        gcagcgcagATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCg                                                                                                       >  2:398007/1‑149 (MQ=255)
                                           gcgcagATGTAAGGTGACAGTGATTGAACTGGCGGTAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGATGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGAt                                                                                                     <  3:723418/148‑1 (MQ=255)
                                           gcgcagATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGAt                                                                                                     <  3:64139/148‑1 (MQ=255)
                                              cagATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTg                                                                                                 <  2:418282/149‑1 (MQ=255)
                                                gATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGaa                                                                                               >  2:445939/1‑149 (MQ=255)
                                                  tGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCTGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGaaaa                                                                                             <  2:560764/149‑1 (MQ=255)
                                                   gTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGaaaa                                                                                             <  4:954264/148‑1 (MQ=255)
                                                    tAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGaaaaa                                                                                            >  3:531412/1‑148 (MQ=255)
                                                         tGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGaa                                                                                      <  1:780086/149‑1 (MQ=255)
                                                          gACAGGGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGaa                                                                                      >  4:16139/1‑148 (MQ=255)
                                                             aGTGATTGAACTGGCGTCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTg                                                                                   <  4:761813/148‑1 (MQ=255)
                                                             aGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGa                                                                                  >  2:720991/1‑149 (MQ=255)
                                                                   tGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTg                                                                             >  3:459122/1‑148 (MQ=255)
                                                                        tGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAg                                                                        <  3:735958/148‑1 (MQ=255)
                                                                        tGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGt                                                                       <  1:398007/149‑1 (MQ=255)
                                                                        gtccatcAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAg                                                                        <  4:1023295/142‑1 (MQ=255)
                                                                             gcAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGa                                                                  >  2:650798/1‑149 (MQ=255)
                                                                               aaCCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTTCTCAATAATGCCATTGAACATGTGGTCGATGGTTAAAAAGTAGAACTGACGCTGCAAAGTGGGGaaa                                                                >  2:391356/1‑149 (MQ=255)
                                                                               aaCCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGaaa                                                                >  1:900051/1‑149 (MQ=255)
                                                                               aaCCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGaa                                                                 <  4:1041540/148‑1 (MQ=255)
                                                                                aCCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGaaa                                                                >  4:479733/1‑148 (MQ=255)
                                                                                   gTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCCGCAAAGTGGGGAAACGCt                                                            >  2:427891/1‑149 (MQ=255)
                                                                                    tCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCt                                                            >  3:607744/1‑148 (MQ=255)
                                                                                     cATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTc                                                          <  2:572859/149‑1 (MQ=255)
                                                                                       tGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCa                                                         <  4:828750/148‑1 (MQ=255)
                                                                                         ggCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTg                                                                      <  4:886502/133‑1 (MQ=255)
                                                                                         ggCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTg                                                                      >  3:886502/1‑133 (MQ=255)
                                                                                             gTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGa                                                   <  4:531412/148‑1 (MQ=255)
                                                                                               aaTGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGAtgt                                                <  2:1066506/149‑1 (MQ=255)
                                                                                                  gCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATgtggt                                             <  1:826405/149‑1 (MQ=255)
                                                                                                   caccacCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCAtt                                                                                                                                      >  4:1041511/1‑59 (MQ=255)
                                                                                                   caccacCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCAtt                                                                                                                                      <  3:1041511/59‑1 (MQ=255)
                                                                                                   caccacCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATgtggtg                                            >  1:635015/1‑149 (MQ=255)
                                                                                                   caccacCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATgtggt                                             <  3:165484/148‑1 (MQ=255)
                                                                                                    accaccGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTa                                                                                         <  1:614895/103‑1 (MQ=255)
                                                                                                    accaccGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTa                                                                                         >  2:614895/1‑103 (MQ=255)
                                                                                                    accaccGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGa                                           >  1:830516/1‑149 (MQ=255)
                                                                                                    accaccGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGa                                           <  2:119124/149‑1 (MQ=255)
                                                                                                    accaccGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGa                                           >  2:48264/1‑149 (MQ=255)
                                                                                                      caccGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGAtt                                         >  2:685307/1‑149 (MQ=255)
                                                                                                      caccGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGAtt                                         >  1:715193/1‑149 (MQ=255)
                                                                                                        ccGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTa                                       >  2:472601/1‑149 (MQ=255)
                                                                                                        ccGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTa                                       >  1:307904/1‑149 (MQ=255)
                                                                                                               tGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTAt                                 <  3:1039279/148‑1 (MQ=255)
                                                                                                                gCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTAtt                                <  3:224019/148‑1 (MQ=255)
                                                                                                                   aCGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGt                             >  3:349612/1‑148 (MQ=255)
                                                                                                                   aCGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTa                            <  2:397631/149‑1 (MQ=255)
                                                                                                                    cGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTAt                           <  1:48264/149‑1 (MQ=255)
                                                                                                                      cTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACgg                                    <  1:76844/138‑1 (MQ=255)
                                                                                                                      cTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACgg                                    >  2:76844/1‑138 (MQ=255)
                                                                                                                         tCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGa                                                                  >  2:508361/1‑105 (MQ=255)
                                                                                                                         tCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGa                                                                  <  1:508361/105‑1 (MQ=255)
                                                                                                                          cTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAgcgc                     >  2:828540/1‑149 (MQ=255)
                                                                                                                                 aGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAg              <  1:673747/149‑1 (MQ=255)
                                                                                                                                  gcgcCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTAc                                      <  4:331846/124‑1 (MQ=255)
                                                                                                                                  gcgcCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTAc                                      >  3:331846/1‑124 (MQ=255)
                                                                                                                                  gcgcCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAg                        >  4:978518/1‑138 (MQ=255)
                                                                                                                                  gcgcCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAg                        <  3:978518/138‑1 (MQ=255)
                                                                                                                                  gcgcCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAg                        <  3:613341/138‑1 (MQ=255)
                                                                                                                                  gcgcCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAAAGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAg                        >  4:613341/1‑138 (MQ=255)
                                                                                                                                   gcacaccaGCCTGCTGGTGGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGaa                                                                                      >  1:616190/4‑74 (MQ=37)
                                                                                                                                   gcacaccaGCCTGCTGGTGGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGaa                                                                                      <  2:616190/71‑1 (MQ=37)
                                                                                                                                     ccaccaGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCAAc          >  2:893224/1‑149 (MQ=255)
                                                                                                                                     acaccaGCCTGCTGGTGGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCAACt         <  2:305159/148‑1 (MQ=255)
                                                                                                                                       accaGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCAACTg        >  2:1039522/1‑149 (MQ=255)
                                                                                                                                        ccaGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCAACTgg       <  2:878321/149‑1 (MQ=255)
                                                                                                                                             agGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCAACTGGCTcgc  >  2:479050/1‑149 (MQ=255)
                                                                                                                                              ggCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCAACTGGCTcgc  >  3:211391/1‑148 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
GTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCAACTGGCTCGC  >  NC_000913/2672498‑2672787

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