Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,449,050 T→A 7.8% intergenic (+32/+26) feaB → / ← tynA phenylacetaldehyde dehydrogenase/copper‑containing amine oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,449,0500TA 7.8% 255.1 / 13.2 95intergenic (+32/+26)feaB/tynAphenylacetaldehyde dehydrogenase/copper‑containing amine oxidase
Reads supporting (aligned to +/- strand):  ref base T (57/28);  new base A (3/4);  total (63/32)
Fisher's exact test for biased strand distribution p-value = 2.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.45e-03

ACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCTTCGGCGCGGGC  >  NC_000913/1448909‑1449189
                                                                                                                                             |                                                                                                                                            
aCGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTg                                                                                                                                             >  2:195025/1‑143 (MQ=255)
  gttAACTTACTGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCGGGAGGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGt                                                                                                                                           <  1:254031/141‑1 (MQ=255)
  gCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGAGGGTTGGTGTGAAACTAAGTCGGTGTGCGTACGGGAGTAATCTGGGTCGCTCATAAGTAACAAAGCGCATCTgat                                                                                                                                           >  2:801917/1‑141 (MQ=255)
  gCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCACCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTGATAAGTAAAAACCGGCACCTGGt                                                                                                                                           >  1:491362/1‑143 (MQ=255)
   cTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTg                                                                                                                                          <  1:665605/143‑1 (MQ=255)
   cTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTAAGGGATTAATCTGGTTCGCTCATAAGTAAAAAACAGCACATGGTg                                                                                                                                          >  1:833031/1‑143 (MQ=255)
        ttACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCGCCGACTGGCTGGACGGGTGGTGGGAAACTAAGTCGGTGAGAGGACGGTATTAATGTGGTTCGCTCATAAGTAAAAAACGGCACGTGGTGCCggt                                                                                                                                     >  2:436741/1‑141 (MQ=255)
        ttACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGtt                                                                                                                                    >  1:58262/1‑143 (MQ=255)
        ttACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGtt                                                                                                                                    >  2:789829/1‑143 (MQ=255)
         tACCGTTTGGTTGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCAGGTGCCGttt                                                                                                                                   <  1:67608/143‑1 (MQ=255)
         tACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCAGGTGCCGttt                                                                                                                                   <  1:407038/143‑1 (MQ=255)
          atcGTTTGGGGGGATCAAGCAGTCAGGAACGGACCGTGATTTTGGCCCCTACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCAGGTGCCGtttt                                                                                                                                  <  2:402055/141‑1 (MQ=255)
             gTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAACAAACGGCACCTGGTGCCGttttttt                                                                                                                               >  1:127525/1‑143 (MQ=255)
             gTTCGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGttttttt                                                                                                                               <  1:19580/143‑1 (MQ=255)
                  gtggGATGAGGCAGGCAGGAGCGGGCCGTGATTTTGGCCCCCACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCAGGTGCCGttttttacttat                                                                                                                          <  1:740374/143‑7 (MQ=255)
                  gtggGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATGAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTg                                                                                                                          >  1:355211/1‑143 (MQ=255)
                   tggGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTGGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGa                                                                                                                         >  2:188449/1‑143 (MQ=255)
                    gggATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGaa                                                                                                                        <  1:423384/143‑1 (MQ=255)
                     tgATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGaaa                                                                                                                       <  2:758196/142‑1 (MQ=255)
                     ggATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGGGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGaaa                                                                                                                       <  2:496225/143‑1 (MQ=255)
                      gATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCCACGGGAGCCGTTGTTTTGTCTGCAAc                                                                                                                      >  2:32563/1‑143 (MQ=255)
                       aTGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGGTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACa                                                                                                                     >  1:741849/1‑143 (MQ=255)
                           aGCAGTCAGGAACGGGCCGTGCTTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTAGCGCAGACGTAAGCGGAGGCGACCGGTGCCGTTTTTTTGTGTGAAACAATCa                                                                                                                 >  2:816449/1‑143 (MQ=255)
                                 cAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTAtc                                                                                                           <  1:528695/143‑1 (MQ=255)
                                    gAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTCTTTTGTCTGAAACAATCACTTAtcttt                                                                                                        >  2:553931/1‑143 (MQ=255)
                                      aCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTCATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGCGCCGTTTTTTTGGCTGGGACAATCACTTAtctacct                                                                                                      >  1:279599/1‑139 (MQ=255)
                                       cGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTAtctttctt                                                                                                     <  1:565967/143‑1 (MQ=255)
                                        gggCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCAACTGGTGCCGTTTGTTTGTGTGAAACAATCCACTATATttcttc                                                                                                    >  2:332586/1‑143 (MQ=255)
                                           ccGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCCCCTGGTGCCGCTTTTTTATCTGAAACAATCACTTATCTTTCTTCCgc                                                                                                 >  2:645988/1‑143 (MQ=255)
                                             gTGATTTTGTCTCCGACTGGCTGGATGGTTCGTGTGATACTAAGTCGGTGTGTGGACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAgcgc                                                                                               <  1:249552/143‑1 (MQ=255)
                                             gTGATTTTGGCCCCGACTGGCTGTACGGTTGGTTGGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACACTCACTTATCTTTCTTCAgcga                                                                                               >  1:189176/1‑142 (MQ=255)
                                             gTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCAACCGGTGCGGCGGCATGGTCTGAAACAATCCCTTATCTTTCTTCAgagc                                                                                               >  2:860485/1‑143 (MQ=255)
                                              tGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGcc                                                                                              >  2:715666/1‑143 (MQ=255)
                                              tGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTAACTTTCTTAAGCGcc                                                                                              >  1:415828/1‑143 (MQ=255)
                                                  tttGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCCCCCCTa                                                                                          >  1:760355/1‑143 (MQ=255)
                                                       ccccGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGtt                                                                                     >  2:278374/1‑143 (MQ=255)
                                                        cccGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCCCTTATCTTTCTTCAGCCCCCCTAGCGTTg                                                                                    >  2:117119/1‑143 (MQ=255)
                                                        cccGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGGTTTTTTGTCTGAAACAATCGATTATCTTTCGTCAGCGCCCCTAggcttc                                                                                    >  1:395169/1‑138 (MQ=255)
                                                            aCTGTCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTATCCTGGTTCGCGCATAAGTAAAAAACGACACCTTGTGCCGTTTTTTTATCTGAAACAATCACTTATCTTTGTTCAGCGTCCCTAGCGTTGGCga                                                                                >  2:848800/1‑142 (MQ=255)
                                                                gctggACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCg                                                                            >  2:685176/1‑143 (MQ=255)
                                                                 ctggACGGTGGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGt                                                                           <  2:235745/143‑1 (MQ=255)
                                                                    gACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAAGAATCACTTATCTTTCTTGAGCGCCCGTAGCGTTGGCGTGTGGTCaa                                                                        >  1:448419/1‑143 (MQ=255)
                                                                     ccgcttGGTGTGAAACTATGTCGGTGTGTGTACGGTATTAATCTGCTTCGCTCATAAGTAAAAAACGGCACCTGCTGCCGTTTTTTTGTCTTAAACAATCACTTATCTTTCTTCAGCGCCTCTAGCGTTGGCGTTTCGTCaaa                                                                       <  2:77694/139‑1 (MQ=255)
                                                                     aCGGTTGGTGTGAAACTAAGTCGGGGTGTGTACGGTATTAATCTCGTTCGCTCATAAGTAAAAAACGGCCCCTGGTGCCGTTTCTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCaaa                                                                       <  2:481712/143‑1 (MQ=255)
                                                                       ggTTGGTGTGAAACTAAGTCGGAGTGTGTACGGTATTAATCTGGTTCGCACATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAaagc                                                                     >  2:584158/1‑142 (MQ=255)
                                                                          tGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGt                                                                  <  2:721529/143‑1 (MQ=255)
                                                                          tGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTACCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGt                                                                  >  1:149127/1‑143 (MQ=255)
                                                                             tgtgAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTcc                                                               >  1:379640/1‑143 (MQ=255)
                                                                             tgtgAAACTAAGTCGGTGTGCGTGCTGCGTTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCCAAGACGTTcc                                                               >  2:538953/1‑143 (MQ=255)
                                                                              gtgAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCa                                                              >  1:562644/1‑143 (MQ=255)
                                                                               tgAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCCAAGAAGTTCCAt                                                             >  2:838282/1‑143 (MQ=255)
                                                                               tgAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCAt                                                             >  2:790046/1‑143 (MQ=255)
                                                                               tgAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGCCGTTTCGTCAAAGAAGTTCCAt                                                             >  2:231879/1‑143 (MQ=255)
                                                                               tgAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTGCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCAt                                                             <  2:690437/143‑1 (MQ=255)
                                                                               tgAAACTAAGTCGGTGTGTGTACGGAATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGGTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCAt                                                             >  1:534085/1‑143 (MQ=255)
                                                                                   aCTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGtt                                                         >  2:475094/1‑143 (MQ=255)
                                                                                    cTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAACGAAGTTCCATGGttt                                                        >  1:137410/1‑143 (MQ=255)
                                                                                       aGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTTTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTGTCGTCAAAGAAGTTCGCTGGTTTcag                                                     >  1:795006/1‑143 (MQ=255)
                                                                                         tCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTcagca                                                   <  1:797184/143‑1 (MQ=255)
                                                                                         tCGGTGCCTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTcagca                                                   <  2:565640/143‑1 (MQ=255)
                                                                                            gtgtgtgtACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCagag                                                >  1:754250/1‑143 (MQ=255)
                                                                                             tgtgtgtACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAACGAAGTTCCATGGGTTACGCagagg                                               >  1:479265/1‑142 (MQ=255)
                                                                                                   tACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTag                                         >  1:275291/1‑142 (MQ=255)
                                                                                                   tACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGCATGTAc                                         >  1:468931/1‑143 (MQ=255)
                                                                                                      ggTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCa                                      >  2:22832/1‑143 (MQ=255)
                                                                                                      ggTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGACGTTTTTTTGTATGAAACAATCACTTATCTTTCTTCAGCGCACCTAGCGTTGGCGTTTCGACAAAGAAGTTGCAAGGTAACAGAAGAGTATGTaacaa                                      >  2:655896/1‑139 (MQ=255)
                                                                                                        tgttAATCTGGTTCTCTCATAGCTAATGAACGGCACCTGGTGCCGTTTTTTTGTCTGACACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCAtt                                    <  2:812704/141‑1 (MQ=255)
                                                                                                        tATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGACGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCACCTAGCGTTGACGTTTCGACAAAGACGATGCATGATTACAGCAGATTCTGTACCCAta                                    >  2:215280/1‑142 (MQ=255)
                                                                                                         aTTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATtc                                   <  1:22249/143‑1 (MQ=255)
                                                                                                         aTTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATtc                                   >  2:784034/1‑143 (MQ=255)
                                                                                                              tttGGTTCGCGCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATtcggtcg                              <  1:507776/141‑1 (MQ=255)
                                                                                                              tCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATtcggtcg                              >  1:226481/1‑143 (MQ=255)
                                                                                                              tCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTGTTTGTCTGAAAAAATCAGTGATGTTTCGTCAGCGCCCCTGGCGTTGGCGTTTCGTAAAAGAAGTTCGAGGGTTTGAGCAGAGGAGGTACCAGTgcggtcg                              >  2:116379/1‑143 (MQ=255)
                                                                                                              tCTGGTGCGCTCATAAGTAAAAAACAGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATtcggtcg                              <  1:231723/143‑1 (MQ=255)
                                                                                                                  gTTCGCTCATAAGTAAAAACCGGCACCTGGTTCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCTTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCAt                          <  2:542830/143‑1 (MQ=255)
                                                                                                                  gTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCAt                          <  2:60822/143‑1 (MQ=255)
                                                                                                                        tCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCgg                    >  2:188314/1‑143 (MQ=255)
                                                                                                                         cATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGc                   <  1:477539/143‑1 (MQ=255)
                                                                                                                          aTAAGTAAAAAACGGCACCTGGTGCCGTATTTTTGTCTGAAACAATCACTTAGCCATCTTCAGCGCCACTAGCGAGAGCGTTTCGTACAAGAAGTGCAATGATTACCGCAGAGAGTGTGCCCGATCGATGGGCATCATCgtac                  >  2:606852/1‑140 (MQ=255)
                                                                                                                           tAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTACATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATCATCGGCCa                 >  2:453201/1‑143 (MQ=255)
                                                                                                                              gTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGCTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCActc              >  1:527953/1‑143 (MQ=255)
                                                                                                                               tAAAAACCGGCCCCTCGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCActct             <  2:605877/143‑1 (MQ=255)
                                                                                                                               tAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCACAGAATTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGGCGGCATAATCGGCCActct             >  1:502400/1‑143 (MQ=255)
                                                                                                                               tAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGATTCAGCAGAGTATGTACCCATTCGGTCGGCTTAATCGGCCActct             >  2:607485/1‑143 (MQ=255)
                                                                                                                               aaaaaaaCGGCACCAGGTGCCGTTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCCTAATCGGGCActc              >  2:154917/2‑143 (MQ=255)
                                                                                                                                aacagaCGGCACCTGGTGGCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCGCCTAGCGTTGGCGGTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCtt            <  2:720941/138‑1 (MQ=255)
                                                                                                                                 aaaaaCGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTAACGTTTCGTCAAAGAAGTTCCCTGGTGTCAGCAGAGTATGTACCGATTCGGTCGGCCTAATCAGACACGCTTc           >  2:704523/1‑143 (MQ=255)
                                                                                                                                  aaaaCGGCACCAGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCTTCg          >  1:273103/1‑143 (MQ=255)
                                                                                                                                    aaCGGCACCAGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCCGCATAATCGGCCCCTCTTCgtc        >  1:593205/1‑141 (MQ=255)
                                                                                                                                     aCGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCTTCGgcg       >  1:632575/1‑143 (MQ=255)
                                                                                                                                      cGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTCTCTTTCTCCAGCGCCCCTAGGATTGGAGGGTTGTCAGAGCAGTCCCCCGGCTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCAACTCCTCGgcgc      >  2:820395/1‑143 (MQ=255)
                                                                                                                                      cGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCTGCCACTCCTCGgcgc      >  2:206200/1‑143 (MQ=255)
                                                                                                                                      cGGCACCTGGCGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCTTCGgcgc      <  1:184102/143‑1 (MQ=255)
                                                                                                                                        gCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCTTCGGCGCgg    <  1:787221/143‑1 (MQ=255)
                                                                                                                                          aCCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCTTCGGCGCGGGc  <  1:278374/143‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                            
ACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCTTCGGCGCGGGC  >  NC_000913/1448909‑1449189

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: