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

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

Only 100 of 101 total aligned reads displayed.

TTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCAACTGGCTCGCGA  >  NC_000913/2672492‑2672789
                                                                                                                                                     |                                                                                                                                                     
ttATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAgg                                                                                                                                                      >  2:1368903/1‑150 (MQ=255)
ttATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTCAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAgg                                                                                                                                                      <  1:975037/150‑1 (MQ=255)
     ggTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGgt                                                                                                                                                 >  2:160134/1‑150 (MQ=255)
      gTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGgtg                                                                                                                                                <  1:593070/150‑1 (MQ=255)
      gTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGgtg                                                                                                                                                >  2:94751/1‑150 (MQ=255)
      gTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGgtg                                                                                                                                                <  2:701665/150‑1 (MQ=255)
      gTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGgtg                                                                                                                                                <  2:614685/150‑1 (MQ=255)
        gCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTGAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCTGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGtgtg                                                                                                                                              <  2:386804/150‑1 (MQ=255)
        gCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGtgtg                                                                                                                                              <  2:1682996/150‑1 (MQ=255)
         ccGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTgc                                                                                                                                             >  1:149788/1‑150 (MQ=255)
          cGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGAGTgcg                                                                                                                                            >  2:376106/1‑150 (MQ=255)
            gAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCa                                                                                                                                          >  1:207349/1‑150 (MQ=255)
            gAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCa                                                                                                                                          <  2:701498/150‑1 (MQ=255)
              aCTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCACACCAGCCTGCTGgtggcgctgt                                                                                                                                         <  2:345498/150‑8 (MQ=255)
               cTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCCACGCTATCTTTTACCGCGC‑CACCAGCAGGCTGGTGTGCGCAttg                                                                                                                                       >  2:1076930/1‑149 (MQ=255)
                   tGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCt                                                                                                                                   <  2:1660280/150‑1 (MQ=255)
                    gggcTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTc                                                                                                                                  <  1:771722/147‑1 (MQ=255)
                     gTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCa                                                                                                                                 <  2:1345355/150‑1 (MQ=255)
                     gTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGAGACAGTGACTGAACTGGCGGCAACCGTCATGGGCCGTACTGCACCACCGCCCGTGCAACGCTTTCTTTTACAGCGC‑AACCCCCAGGCTGCGGTGCGCATTCTGCtaa                                                                                                                                 >  1:1045055/1‑148 (MQ=255)
                        tGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCaata                                                                                                                              <  2:1234910/150‑1 (MQ=255)
                            aCTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGc                                                                                                                          <  2:976929/150‑1 (MQ=255)
                                gCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGACGTAATGCACCACCGCACGTGCAACGCAAGCATTTACAGCGC‑CACCAGCATGCTGGTGTGCGCATTCTGCTCAATAAAGCCAtg                                                                                                                      >  2:1656788/1‑149 (MQ=255)
                                  tGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCGGCTCAATAATGCCATTGa                                                                                                                    >  2:1371940/1‑150 (MQ=255)
                                   gCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCAGTGaa                                                                                                                   >  1:1003157/1‑150 (MQ=255)
                                     cAGCGCCTCGCAGCGCAGATGGAAGGTGAGAGTGATTGAACTGGCCGCAACCGTCATGGGCCGTAATGCACCACCGCCGGTGCATCGCTATCTTTTACAGCGC‑CACCAGCATGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAAct                                                                                                                 <  2:291559/150‑2 (MQ=255)
                                      aGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACAt                                                                                                                <  1:1619880/150‑1 (MQ=255)
                                      aGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACAt                                                                                                                >  2:750186/1‑150 (MQ=255)
                                         gcCACGCAGCGCAGATGTGAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACAtgtg                                                                                                             <  2:1265034/150‑1 (MQ=255)
                                         gcCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACAtgtg                                                                                                             <  2:191583/150‑1 (MQ=255)
                                         gcCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACAtgtg                                                                                                             >  2:1035273/1‑150 (MQ=255)
                                           cACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGGTCAATAATGCCATTGATCATGTGGt                                                                                                           >  2:70156/1‑150 (MQ=255)
                                           cACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTg                                                                                                                     >  1:530611/1‑140 (MQ=255)
                                           cACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTg                                                                                                                     <  2:530611/140‑1 (MQ=255)
                                             cgcagcgcagATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCg                                                                                                         >  2:186579/1‑150 (MQ=255)
                                                agcgcagATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATg                                                                                                      >  1:1448294/1‑150 (MQ=255)
                                                       aTGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCTC‑CACCAGCAGGCTGTTGTGCGCATTCTGCTCAATAATGCCATCGACCACATGTTCAATggcattat                                                                                               <  2:344272/150‑7 (MQ=255)
                                                         gTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAg                                                                                             <  2:712361/150‑1 (MQ=255)
                                                          tAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGt                                                                                            >  1:38890/1‑150 (MQ=255)
                                                            aGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAg                                                                                          >  2:1413721/1‑150 (MQ=255)
                                                            aGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGAGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAg                                                                                          <  1:1267597/150‑1 (MQ=255)
                                                             ggTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGa                                                                                         <  1:1333472/150‑1 (MQ=255)
                                                                gACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACt                                                                                      <  1:379792/150‑1 (MQ=255)
                                                                  cAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGa                                                                                    >  1:207271/1‑150 (MQ=255)
                                                                    gTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACg                                                                                  <  2:789182/150‑1 (MQ=255)
                                                                       aTTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTg                                                                               <  2:1066413/150‑1 (MQ=255)
                                                                       aTTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTg                                                                               >  2:597933/1‑150 (MQ=255)
                                                                       aTTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTg                                                                               <  1:169847/150‑1 (MQ=255)
                                                                         tGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCa                                                                             >  2:646890/1‑150 (MQ=255)
                                                                             cTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGt                                                                         >  1:1312630/1‑150 (MQ=255)
                                                                             cTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGt                                                                         >  1:785483/1‑150 (MQ=255)
                                                                             cTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGt                                                                         >  1:792070/1‑150 (MQ=255)
                                                                               ggcggcAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTgg                                                                       <  1:806545/150‑1 (MQ=255)
                                                                                 cggcAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTgggg                                                                     <  1:94751/150‑1 (MQ=255)
                                                                                    catccGTCATGGGCCGTAATGCACCACCGCCCGTGGAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCTTTTAACATGTGGTGTATGGTGAATAAGTAGAACTGACGCTGCAAAGTGGGGata                                                                  <  2:527504/147‑3 (MQ=255)
                                                                                     aaCCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAAc                                                                 >  1:1549355/1‑150 (MQ=255)
                                                                                     aaCCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAAc                                                                 >  1:53688/1‑150 (MQ=255)
                                                                                     aaCCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAAc                                                                 <  2:149788/150‑1 (MQ=255)
                                                                                      aCCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGAGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACg                                                                <  1:376106/150‑1 (MQ=255)
                                                                                       ccGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGc                                                               <  1:505658/150‑1 (MQ=255)
                                                                                          tCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATACTGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGCCGCTACAAAGTGGGGAAACGCTTc                                                            <  1:1191774/150‑1 (MQ=255)
                                                                                           cATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCa                                                           >  2:1416807/1‑150 (MQ=255)
                                                                                           cATGGGCCGTAAAGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCAGCACAATAATGCCTTTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTCCAAAGTGGGGAAACGCTTCa                                                           >  2:405566/1‑150 (MQ=255)
                                                                                             tGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAgg                                                         >  2:450135/1‑150 (MQ=255)
                                                                                                gCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTg                                                      >  1:1402669/1‑150 (MQ=255)
                                                                                                gCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTg                                                      >  2:262616/1‑150 (MQ=255)
                                                                                                   gTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGAtg                                                   >  2:386984/1‑150 (MQ=255)
                                                                                                    tAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGACTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGAtgt                                                  >  2:1517237/1‑150 (MQ=255)
                                                                                                      aTGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATgtgg                                                >  2:739323/1‑150 (MQ=255)
                                                                                                      aTGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTCAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATgtgg                                                >  1:319293/1‑150 (MQ=255)
                                                                                                       tGCACCACCGCCCGTGCATCGCTATCTTTTACAGCGC‑CACCCGCAGGCTGGTGTGCGCATTCAGCTCAATACTGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATgtggt                                               <  1:252588/150‑1 (MQ=255)
                                                                                                          accaccGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGAt                                            <  2:1428452/150‑1 (MQ=255)
                                                                                                           ccaccGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGAtt                                           >  2:1608934/1‑150 (MQ=255)
                                                                                                            caccGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGAttt                                          >  1:1498267/1‑150 (MQ=255)
                                                                                                               cGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACg                                       <  2:1258928/150‑1 (MQ=255)
                                                                                                               cGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACg                                       <  1:927143/150‑1 (MQ=255)
                                                                                                                gCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACgg                                      <  1:160134/150‑1 (MQ=255)
                                                                                                                 cccGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGt                                     <  2:1459750/150‑1 (MQ=255)
                                                                                                                  ccGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTa                                    >  1:1232765/1‑150 (MQ=255)
                                                                                                                   cGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTAt                                   >  1:640680/1‑150 (MQ=255)
                                                                                                                      gCATCGCTTTCTTTTACAGCGC‑CACCAGCAGGCTGGGGTGCGCATTCTGCTCAATAATGCCCTTGAACATGTGGTCGAGGGGGAAACAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTTTTgg                                <  2:17959/150‑1 (MQ=255)
                                                                                                                      gCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTgg                                <  2:1189237/150‑1 (MQ=255)
                                                                                                                        aaCGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTa                              <  2:430037/150‑1 (MQ=255)
                                                                                                                          cGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATc                            <  2:207271/150‑1 (MQ=255)
                                                                                                                            cTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAg                          >  1:1189236/1‑150 (MQ=255)
                                                                                                                             tATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAgc                         >  2:1085378/1‑150 (MQ=255)
                                                                                                                              aTCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAgcg                        <  2:1362301/150‑1 (MQ=255)
                                                                                                                                cTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGcc                      >  2:424395/1‑150 (MQ=255)
                                                                                                                                    tACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACg                  <  1:1386871/150‑1 (MQ=255)
                                                                                                                                       aGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGc               >  2:468946/1‑150 (MQ=255)
                                                                                                                                       aGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGc               <  2:1448294/150‑1 (MQ=255)
                                                                                                                                         cgcacaccaGCCTGCTGGTGTGCGCATTCTGCTCAATAAAGCCATTGAACAAGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAAAGAGCa              >  1:144970/5‑150 (MQ=255)
                                                                                                                                           cacaccaGCCTGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCAAc            >  1:1040482/3‑150 (MQ=255)
                                                                                                                                           cacaccaGCCTGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCAAc            >  2:307597/3‑150 (MQ=255)
                                                                                                                                           cacaccaGCCTGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCAAc            >  1:1398139/3‑150 (MQ=255)
                                                                                                                                             caccaGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCAACTg          >  1:43154/1‑150 (MQ=255)
                                                                                                                                                cagcagGCTGGTGTGCGCATTCTGCTCAATACTTCCATTGAACATGTGGTCGATGGGGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTAGCAGCGCCAACGAGCAACTGGCt       <  1:270134/150‑1 (MQ=255)
                                                                                                                                                   cagGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCAACTGGCTcgc    >  1:1074304/1‑150 (MQ=255)
                                                                                                                                                     ggCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCAACTGGCTCGCGa  >  2:651099/1‑150 (MQ=255)
                                                                                                                                                     ggCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCAACTGGCTCGCGa  >  2:1352346/1‑150 (MQ=255)
                                                                                                                                                     ggCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCAACTGGCTCGCGa  >  2:1059506/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGC‑CACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCAACTGGCTCGCGA  >  NC_000913/2672492‑2672789

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