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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,449,0500TA17.4% 337.5 / 55.6 139intergenic (+32/+26)feaB/tynAphenylacetaldehyde dehydrogenase/tyramine oxidase, copper‑requiring
Reads supporting (aligned to +/- strand):  ref base T (52/61);  new base A (9/15);  total (62/77)
Fisher's exact test for biased strand distribution p-value = 5.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 143 total aligned reads displayed.

GCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGT  >  NC_000913/1448951‑1449148
                                                                                                   |                                                                                                   
gCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCt                                                                                                     >  2:2739529/1‑100 (MQ=255)
gCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCt                                                                                                     <  1:2126929/100‑1 (MQ=255)
gCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCt                                                                                                     <  1:3802981/100‑1 (MQ=255)
 ccGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTg                                                                                                    >  1:1872742/1‑100 (MQ=255)
     gATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGc                                                                                                <  1:2451341/100‑1 (MQ=255)
      aTTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCTGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGcc                                                                                               <  2:1838063/100‑1 (MQ=255)
       ttttGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCg                                                                                              <  2:1206933/100‑1 (MQ=255)
       ttttGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCg                                                                                              >  1:162312/1‑100 (MQ=255)
        tttGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑t                                                                                            >  2:2567587/1‑100 (MQ=255)
           ggCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCAGGTGCCG‑tttt                                                                                         <  1:2659476/100‑1 (MQ=255)
            gccggggctgggtggAGAGTGGGTGTGAAGATAAGTGGGTGTGTGTAGGGTATTAATCTGATTCGCTCATAAGTAAAAAACGGCACTAGGTGCCG‑tttt                                                                                         <  2:2278763/92‑1 (MQ=255)
            gCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑ttttt                                                                                        <  1:576219/100‑1 (MQ=255)
            gCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑ttttt                                                                                        <  1:2997039/100‑1 (MQ=255)
            gCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCAGGTGCCG‑ttttt                                                                                        <  1:3237633/100‑1 (MQ=255)
            gCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCAGGTGCCG‑ttttt                                                                                        <  1:221630/100‑1 (MQ=255)
             ccccGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑tttttt                                                                                       <  1:1889251/100‑1 (MQ=255)
              cccGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑ttttttt                                                                                      <  1:1679995/100‑1 (MQ=255)
               ccGACTGGCTGGACGGTTGGTGTGAAACTAAGGCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTg                                                                                     <  2:398965/100‑1 (MQ=255)
                cggctggctggACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑tttttt                                                                                       <  2:3754555/94‑1 (MQ=255)
                cGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCAGGTGCCG‑ttttttact                                                                                    <  1:1090236/100‑4 (MQ=255)
                 gACTGGCTGGACGGTTGGTGTGAATCTAAGTCTGTGTCTGTACTGTATTAATCTGGTTCGCTCATAAGTAAACCCCGTCACCAGGTGCCG‑ttttttactt                                                                                   <  2:2581783/100‑5 (MQ=255)
                 gACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTc                                                                                   <  1:2763293/100‑1 (MQ=255)
                 gACGGGCTGTACGGGTTGTGTGAGACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCAGGTGCCG‑ttttttact                                                                                    <  1:3571586/99‑4 (MQ=37)
                   ctggctggACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCAGGTGCCG‑ttttttactta                                                                                  <  2:1333576/99‑6 (MQ=255)
                       ctgtACGTTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTTATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAAc                                                                             <  2:916981/100‑1 (MQ=255)
                       ctggACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAAc                                                                             >  2:2108867/1‑100 (MQ=255)
                       ctggACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAAc                                                                             <  1:2897393/100‑1 (MQ=255)
                        tggACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACa                                                                            <  2:1310220/100‑1 (MQ=255)
                        tggACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACa                                                                            <  2:3760510/100‑1 (MQ=255)
                         ggACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACaa                                                                           >  2:1986763/1‑100 (MQ=255)
                           aCGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACAATc                                                                         <  1:2739549/100‑1 (MQ=255)
                           aCGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACAATc                                                                         >  1:1902803/1‑100 (MQ=255)
                             ggTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACAATCAc                                                                       >  1:523256/1‑100 (MQ=255)
                              gTTGGTGTGACACTAAGTCGGTGTGTGTACGGTCTTAATCTGGTTAGCTCATAAGTAAAACACCGCACACGCTGCAG‑TCTTTTTTCCTGAAACaacccc                                                                       >  2:3270991/1‑95 (MQ=255)
                              gTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACAATCACt                                                                      <  1:3295141/100‑1 (MQ=255)
                              gTTGGTGTGAAACTAAGGCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACAATCACt                                                                      >  1:2867632/1‑100 (MQ=255)
                                tGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTa                                                                    <  1:856532/100‑1 (MQ=255)
                                  gtgtGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTAtc                                                                  >  1:1169131/1‑100 (MQ=255)
                                    gtgAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTAtctt                                                                <  1:1986774/100‑1 (MQ=255)
                                    gtgAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTAtctt                                                                >  2:3477783/1‑100 (MQ=255)
                                    gtgAAACTAAGTCGGTGTGTGTACGGGATTAATCTGGTTCGCTCATAAGTAACAAACGGGACGAGGTGCCG‑TTTTGTTGTCTGAAACAAGCACTTAtctt                                                                >  2:3576443/1‑100 (MQ=255)
                                     tgAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTAtcttt                                                               >  1:3034867/1‑100 (MQ=255)
                                      gAAACTAAGTCGGTGTGTGTACGGTGTTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTAtcttt                                                               >  2:1797141/1‑99 (MQ=255)
                                      gAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTAtctttc                                                              >  2:3399303/1‑100 (MQ=255)
                                       aaaCTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTAtctttct                                                             <  2:3252446/100‑1 (MQ=255)
                                        aaCTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTAtctttctt                                                            >  1:799876/1‑100 (MQ=255)
                                         aCTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTATCTttcttc                                                           <  2:3034843/100‑1 (MQ=255)
                                          cTAAGTCGGTGTGTGTACGGTATTAATCTGGTTGGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTATGTCTGAAACAATAACTTATCTTTCTTCa                                                          >  2:2290563/1‑100 (MQ=255)
                                          cTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTATCTTTCTTCa                                                          <  2:165928/100‑1 (MQ=255)
                                             atgcGGTGTGTGGACGGTTTTAGGCTGGTTCGTTCAGGGGTAAGAGAGGGCTTCTGGAGCAG‑TTTTTATGTGTGAAACAATCACTTATCTTTCTTCAgcg                                                       >  2:2675627/4‑100 (MQ=255)
                                             aGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAgcg                                                       <  1:2370895/100‑1 (MQ=255)
                                              gTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAgcgc                                                      >  2:1367834/1‑100 (MQ=255)
                                                 ggTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGcccc                                                   <  2:2157124/100‑1 (MQ=255)
                                                      gtgtACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCg                                              <  1:436379/100‑1 (MQ=255)
                                                      gtgtACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCg                                              <  1:436390/100‑1 (MQ=255)
                                                        gtACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGtt                                            >  1:1054820/1‑100 (MQ=255)
                                                        gtACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGtt                                            <  1:3167863/100‑1 (MQ=255)
                                                          aCGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTgg                                          >  2:1931014/1‑100 (MQ=255)
                                                           cGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGc                                         <  1:2187587/100‑1 (MQ=255)
                                                           cGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGc                                         <  2:2430285/100‑1 (MQ=255)
                                                                ttAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTc                                    <  1:1464977/100‑1 (MQ=255)
                                                                ttAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTc                                    <  1:2334712/100‑1 (MQ=255)
                                                                 tAATCTGGTTCGCTCATAAGTAAAAACCGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCg                                   <  1:3280266/100‑1 (MQ=255)
                                                                 tAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCg                                   >  1:472065/1‑100 (MQ=255)
                                                                  aaTCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGt                                  >  1:3498352/1‑100 (MQ=255)
                                                                   aTCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTc                                 <  1:3280626/100‑1 (MQ=255)
                                                                    tCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCa                                <  2:1780192/100‑1 (MQ=255)
                                                                    tCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCa                                <  1:2567604/100‑1 (MQ=255)
                                                                     cTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCaa                               <  2:2887089/100‑1 (MQ=255)
                                                                      tGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCaaa                              <  2:116348/100‑1 (MQ=255)
                                                                      tGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCaaa                              <  2:691387/100‑1 (MQ=255)
                                                                      tGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCaaa                              >  2:204985/1‑100 (MQ=255)
                                                                       ggTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCaaa                              >  2:867244/1‑99 (MQ=255)
                                                                          tCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAaagaag                          >  2:2565601/1‑100 (MQ=255)
                                                                           cGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGt                         <  1:2113544/100‑1 (MQ=255)
                                                                           cGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGt                         >  1:792231/1‑100 (MQ=255)
                                                                            gCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGtt                        >  2:535333/1‑100 (MQ=255)
                                                                            gCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGtt                        >  1:1537793/1‑100 (MQ=255)
                                                                             cTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTc                       >  2:1540839/1‑100 (MQ=255)
                                                                               cATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCa                     >  1:1956298/1‑100 (MQ=255)
                                                                                  gacaaaaaaaCGGCACCAGGTGCCGTTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATg                   >  2:2815775/5‑100 (MQ=38)
                                                                                  gacaaaaaaaCGGCACCAGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATgg                  >  1:3590651/5‑100 (MQ=255)
                                                                                   acaaaaaaaCGGCACCAGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGt                 >  1:873951/4‑100 (MQ=255)
                                                                                   aGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGt                 >  2:2397028/1‑100 (MQ=255)
                                                                                    caaaaaaaCGGCACCAGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGtt                >  2:1674373/3‑100 (MQ=255)
                                                                                    caaaaaaaCGGCACCAGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGtt                >  1:1558831/3‑100 (MQ=255)
                                                                                    caaaaaaaCGGCACCAGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGtt                >  1:1049921/3‑100 (MQ=255)
                                                                                      aaaactcGGCTCCTGGTGCCG‑TTTTTTTGTCTGAAACATTCACTTATCTTTCTTCGGCGCCCATGGAGTTGGCGTTTCGTCAAGGTAGTTCCATGGTTTc              <  2:3188144/94‑1 (MQ=255)
                                                                                      aaaaaaCGGCACCAGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTc              >  2:3121628/1‑100 (MQ=255)
                                                                                         aaaCGGCACCAGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTcagc           >  2:1193516/1‑100 (MQ=255)
                                                                                          aaCGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTcagca          >  2:694157/1‑100 (MQ=255)
                                                                                            cGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTcagcag         <  1:3345402/99‑1 (MQ=255)
                                                                                             ggCACCTGGTGCCG‑TTTTTTTGTCTGATACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGAGTTTCGTCAAAGAAGTTCGATGGTTTCAGCagag       >  2:2009951/1‑100 (MQ=255)
                                                                                               cACCAGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTa     <  1:1674383/100‑1 (MQ=255)
                                                                                               cACCAGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTa     <  1:1193525/100‑1 (MQ=255)
                                                                                                aCCTGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTAt    <  1:2893396/100‑1 (MQ=255)
                                                                                                 ccTGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATg   <  1:3349363/100‑1 (MQ=255)
                                                                                                  cTGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGt  >  1:1286187/1‑100 (MQ=255)
                                                                                                  cTGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGGTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGt  >  1:3510694/1‑100 (MQ=255)
                                                                                                   tGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATg   <  2:1928643/98‑1 (MQ=255)
                                                                                                   |                                                                                                   
GCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCG‑TTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGT  >  NC_000913/1448951‑1449148

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