Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,449,473 C→T 5.6% R626H (CGC→CAC)  tynA ← tyramine oxidase, copper‑requiring

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,449,4730CT 5.6% 226.8 / 2.5 71R626H (CGC→CAC) tynAtyramine oxidase, copper‑requiring
Reads supporting (aligned to +/- strand):  ref base C (29/38);  new base T (2/2);  total (31/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.53e-01

GGATGATAACGCGTTACCCAGAGCTGCTTGTCCATAAAGCTTAAACGATGATAGATCCACTCGTCCGGCGCGAACTGGGCACCTTTTGCTACCGGGTGAGTACCACCTGCATAAGGAATAATTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGTTACTCAACAGACGAATCGTGCCCGGATCAAATTTCTGTGCGGCATCCTGTTCATTGCCGATGTTGTACTGATTAACTTGCATGGTACTGGTGCGTGGGCCACCGGCAGTATTCGGTTTTACCACTG  >  NC_000913/1449326‑1449617
                                                                                                                                                   |                                                                                                                                                
ggATGATAACGCGTTACCCAGAGCTGCTTGTCCATAAAGCTTAAACGATGATAGATCCACTCGTCCGGCGCGAACTGGGCACCTTTTGCTACCGGGTGAGTACCACCTGCATAAGGAATAATTTGATAGGAAACCGGATTGCCCATGc                                                                                                                                                  >  8:327704/1‑148 (MQ=255)
 gatgatAACGCGTTACCCAGAGCTGCTTGTCCATAAAGCTTAAACGATGATAGATCCACTCGTCCGGCGCGAACTGGGCACCTTTTGCTACCGGGTGAGTACCACCTGCATAAGGAATAATTTGATAGGAAACCGGATTGCCCATGCg                                                                                                                                                 >  7:160078/1‑148 (MQ=255)
          gcgTTACCCAGAGCTGCTTGTCCATAAAGCTTAAACGATGATAGATCCACTCGTCCGGCGCGAACTGGGCACCTTTTGCTACCGGGTGAGTACCACCTGCATAAGGAATAATTTGATAGGAAACCGGATTGCCCATGCGGTTCTCttt                                                                                                                                        >  3:124907/1‑148 (MQ=255)
             ttACCCAGAGCTGCTTGTCCATAAAGCTTAAACGATGATAGATCCACTCGTCCGGCGCGAACTGGGCACCTTTTGCTACCGGGTGAGTACCACCTGCATAAGGAATAATTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGtt                                                                                                                                     <  3:331693/148‑1 (MQ=255)
                  cAGAGCTGCTTGTCCATAAAGCTTAAACGATGATAGATCCACTCGTCCGGCGCGAACTGGGCACCTTTTGCTACCGGGTGAGTACCACCTGCATAAGGAATAATTTGATAGGAAACCGGATTGCCCATGTGGTTCTCTTTGTTCGGGt                                                                                                                                >  8:316329/1‑148 (MQ=255)
                    gagCTGCTTGTCCATAAAGCTTAAACGATGATAGATCCACTCGTCCGGCGCGAACTGGGCACCTTTTGCTACCGGGTGAGTACCACCTGCATAAGGAATAATTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGTTa                                                                                                                              <  4:124907/148‑1 (MQ=255)
                          ctTGTCCATAAAGCTTAAACGATGATAGATCCACTCGTCCGGCGCGAACTGGGCACCTTTTGCTACCGGGTGAGTACCACCTGCATAAGGAATAATTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGTTACTCaa                                                                                                                         <  7:367863/147‑1 (MQ=255)
                           ttGTCCATAAAGCTTAAACGATGATAGATCCACTCGTCCGGCGCGAACTGGGCACCTTTTGCTACCGGGTGAGTACCACCTGCATAAGGAATAATTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGTTACTCaa                                                                                                                         >  2:410419/1‑146 (MQ=255)
                              tCCATAAAGCTTAAACGATGATAGATCCACTCGTCCGGCGCGAACTGGGCACCTTTTGCTACCGGGTGAGTACCACCTGCATAAGGAATAATTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGTTACTCAACAGAc                                                                                                                    >  2:183486/1‑148 (MQ=255)
                                cATAAAGCTTAAACGATGATAGATCCACTCGTCCGGCGCGAACTGGGCACCTTTTGCTACCGGGTGAGTACCACCTGCATAAGGAATAATTTGATAGGAAACCGGATTGCCCATGCGGt                                                                                                                                               >  7:207632/1‑119 (MQ=255)
                                cATAAAGCTTAAACGATGATAGATCCACTCGTCCGGCGCGAACTGGGCACCTTTTGCTACCGGGTGAGTACCACCTGCATAAGGAATAATTTGATAGGAAACCGGATTGCCCATGCGGt                                                                                                                                               <  8:207632/119‑1 (MQ=255)
                                  tAAAGCTTAAACGATGATAGATCCACTCGTCCGGCGCGAACTGGGCACCTTTTGCTACCGGGTGAGTACCACCTGCATAAGGAATAATTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGTTACTCAACAGACGAAt                                                                                                                <  6:425749/148‑1 (MQ=255)
                                   aaaGCTTAAACGATGATAGATCCACTCGTCCGGCGCGAACTGGGCACCTTTTGCTACCGGGTGAGTACCACCTGCATAAGGAATAATTTGATATGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGtt                                                                                                                               >  5:393509/1‑132 (MQ=255)
                                   aaaGCTTAAACGATGATAGATCCACTCGTCCGGCGCGAACTGGGCACCTTTTGCTACCGGGTGAGTACCACCTGCATAAGGAATAATTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGtt                                                                                                                               >  6:127517/1‑132 (MQ=255)
                                   aaaGCTTAAACGATGATAGATCCACTCGTCCGGCGCGAACTGGGCACCTTTTGCTACCGGGTGAGTACCACCTGCATAAGGAATAATTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGtt                                                                                                                               <  6:393509/132‑1 (MQ=255)
                                   aaaGCTTAAACGATGATAGATCCACTCGTCCGGCGCGAACTGGGCACCTTTTGCTACCGGGTGAGTACCACCTGCATAAGGAATAATTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGtt                                                                                                                               <  5:127517/132‑1 (MQ=255)
                                     aGCTTAAACGATGATAGATCCACTCGTCCGGCGCGAACTGGGCACCTTTTGATACCGGGTGAGTACCACCTGCATAAGGAATAATTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGTTACTCAACAGAAGAATCGt                                                                                                             >  2:172606/1‑148 (MQ=255)
                                           aaCGATGATAGATCCACTCGTCCGGCGCGAACTGGGCACCTTTTGCTACCGGGTGAGTACCACCTGCATAAGGAATAATTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGTTACTCAACAGACGAATCGTGCCCgg                                                                                                       <  6:177231/148‑1 (MQ=255)
                                           aaCGATGATAGATCCACTCGTCCGGCGCGAACTGGGCACCTTTTACTACCGGGTGAGTACCACCTGCATAAGGAATAATTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGTTACTCAACAGACGAATCGTGCCCg                                                                                                        >  2:365319/1‑147 (MQ=255)
                                               atgatAGATCCACTCGTCCGGCGCGAACTGGGCACCTTTTGCTACCGGGTGAGTACCACCTGCATAAGGAATAATTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGTTACTCAACAGACGAATCGTGCCCGGATCa                                                                                                   <  3:219450/148‑1 (MQ=255)
                                                 gatAGATCCACTCGTCCGGCGCGAACTGGGCACCTTTTGCTACCGGGTGAGTACCACCTGCATAAGGAATAATTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGTTACTCAACAGACGAATCGTGCCCGGATCaaa                                                                                                 >  3:69160/1‑148 (MQ=255)
                                                     gATCCACTCGTCCGGCGCGAACTGGGCACCTTTTGCTACCGGGTGAGTACCACCTGCATAAGGAATAATTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGTTACTCAACAGACGAATCGTGCCCGGATCAAAttt                                                                                              <  8:160078/147‑1 (MQ=255)
                                                     gATCCACTCGTCCGGCGCGAACTGGGCACCTTTTGCTACCGGGTGAGTACCACCTGCATAAGGAATAATTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGTTACTCAACAGACGAATCGTGCCCGGATCAAATATc                                                                                             >  8:390027/1‑148 (MQ=255)
                                                       tCCACTCGTCCGGCGCGAACTGGGCACCTTTTGCTACCGGGTGAGTACCACCTGCATAAGGAATAATTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGTTACTCAACAGACGAATCGTGCCCGGATCAAATTTCtg                                                                                           <  6:156212/148‑1 (MQ=255)
                                                           cTCGTCCGGCGCGAACTGGGCACCTTTTGCTACCGGGTGAGTACCACCTGCATAAGGAATAATTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGTTACTCAACAGACGAATCGTGCCCGGATCAAATTTCTGTGCg                                                                                       <  4:410897/148‑1 (MQ=255)
                                                                  ggCGCGAACTGGGCACCTTTTGCTACCGGGTGAGTACCACCTGCATAAGGAATAATTTGATAGGAAACCGGATTGCCCATGTGGTTCTCTTTGTTCGGGTTACTCAACAGACGAATCGTGCCCGGATCAAATTTCTGTGCGGCATCCt                                                                                <  7:316329/148‑1 (MQ=255)
                                                                   gcgcGAACTGGGCACCTTTTGCTACCGGGTGAGTACCACCTGCATAAGGAATAATTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGTTACTCAACAGACGAATCGTGCCCGGATCAAATTTCTGTGCGGCATCCTg                                                                               >  5:390151/1‑148 (MQ=255)
                                                                     gcgAACTGGGCACCTTTTGCTACCGGGTGAGTACCACCTGCATAAGGAATAATTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGTTACTCAACAGACGAATCGTGCCCGGATCAAATTTCTGTGCGGCATCCTGtt                                                                             <  6:390151/148‑1 (MQ=255)
                                                                      cgAACTGGGCACCTTTTGCTACCGGGTGAGTCCCACCTGCATAAGGAATAATTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTGCGGGTTACTCAACAGACGAATCGTGCCCGGATCAAATTTCTGTGCGGCATCCTGTTc                                                                            <  4:361132/148‑1 (MQ=255)
                                                                      cgAACTGGGCACCTTTTGCTACCGGGTGAGTACCACCTGCATAAGTAATAATTTGATAGGAAACCGTATTGCCCATGCGGTTCTCTTTGTTCGGGTTACTCAACAGCCGAATCGTGCCCGGATCAAATTTCTGTGCGGCATCCTGTTc                                                                            <  6:310377/148‑1 (MQ=255)
                                                                      cgAACTGGGCACCTTTTGCTACCGGGTGAGTACCACCTGCATAAGGAATAATTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGTTACTCAACAGACGACTCGTGCCCGGATCAAATTTCTGTGCGGCATCCTGTTc                                                                            <  6:316121/148‑1 (MQ=255)
                                                                      cgAACTGGGCACCTTTTGCTACCGGGTGAGTACCACCTGCATAAGGAATAAATAGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGTTCCTCAACAGACGAATCGTGCCCGGATCAAATTTCTGTGCGGCATCCTGtt                                                                             >  7:186552/1‑147 (MQ=255)
                                                                        aaCTGGGCTCCTTTTGCTAGCGGGTGAGTAGCACCGGCATAAGGAATAATTTGATAGGAAACAGGATTGCCCATGCGGTTCTCTTTGTTCGCGGTACTCAACAGACGAATCGTGCCCGGATCAAATTTCTGTGCGGCATCCTGTTCAt                                                                          <  5:75537/148‑1 (MQ=255)
                                                                         aCTGGGCACCTTTTGCTACCGGGTGAGTACCACCTGCATAAGGAATAATTTGATAGGAAACCGGATTGCCCATGTGGTTCTCTTTGTTCGGGTTACTCAACAGACGAATCGTGCCCGGATCAAATTTCTGTGCGGCATCCTGTTCAtt                                                                         >  3:179669/1‑148 (MQ=255)
                                                                             ggCACCTTTTGCTACGGGGTGAGTACCACCTGCATAAGGAATAATTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGTTACTCAACAGACGAATCGTGCCCGGATCAAATTTCTGTGCGGCATCCTGTTCATTGcc                                                                      <  7:390027/147‑1 (MQ=255)
                                                                              gccccTTTTGCTACCGGGAGAGTACCACCTGCATAACGAATAATTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGTTACTCAACAGACGAATCGTGCCCGGATCAAATTTCTGTGCGGCATCCTGTTCATTGCCGa                                                                    <  8:476448/145‑1 (MQ=255)
                                                                              gCACCTTTTGCTACCGGGTGAGTACCACCTGCATAAGGAATAATTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGTTACTCAACAGACGAATCGTGCCCGGATCAAATTTCTGTGCGGCATCCTGtt                                                                             >  1:408205/1‑139 (MQ=255)
                                                                              gCACCTTTTGCTACCGGGTGAGTACCACCTGCATAAGGAATAATTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGTTACTCAACAGACGAATCGTGCCCGGATCAAATTTCTGTGCGGCATCCTGtt                                                                             <  2:408205/139‑1 (MQ=255)
                                                                                aCCTTTTGCTACCGGGTGAGTACCACCTGCATAAGGAATAATTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGTTACTCAACAGACGAATCGTGCCCGGATCAAATTTCTGTGCGGCATCCTGTTCATTGCCGAtg                                                                  >  1:403777/1‑148 (MQ=255)
                                                                                 ccTTTTGCTACCGGGTGAGTACCACCTGCATAAGGAATAATTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGTTACTCAACAGACGAATCGTGCCCGGATCAAATTTCTGTGCGGCATCCTCTTCATTGCCGAtgt                                                                 >  7:370089/1‑148 (MQ=255)
                                                                                  cTTTTGCTACCGGGTGAGTACCACCTGCATAAGGAATAATTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGTTACTCAACAGACGAATCGTGCCCGGATCAAATTTCTGTGCGGCATCCTGTTCATTGCCGAtgt                                                                 >  5:236632/1‑147 (MQ=255)
                                                                                  cTTTTGCTACCGGGTGAGTACCACCTGCATAAGGAATAATTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGTTACTCAACAGACGAATCGTGCCCGGATCAAATTTCTGTGCGGCATCCTGTTCATTGCCGAtgt                                                                 <  6:236632/147‑1 (MQ=255)
                                                                                    tttGCTACCGGGTGAGTACCACCTGCATAAGGAATAATTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGTTACTCAACAGACGAATCGTGCCCGGATCAAATTTCTGTGCGGCATCCTGTTCATTGCCGAtgt                                                                 >  8:289988/1‑145 (MQ=255)
                                                                                     ttGCTACCGGGTGAGTACCACCTGCATAAGGAATAATTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGTTACTCAACAGACGAATCGTGCCCGGATCAAATTTCTGTGCGGCATCCTGTTCATTGCCGAtgttgt                                                              <  3:27755/147‑1 (MQ=255)
                                                                                     ttGCTACCGGGTGAGTACCACCTGCATAAGGAATAATTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGTTACTCAACAGACGAATCGTGCCCGGATCAAATTTCTGTGCGGCATCCTGTTCATTGCCGAtgtt                                                                <  7:289988/145‑1 (MQ=255)
                                                                                               gTGAGTACCACCTGCATAAGGAATAATTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGTTACTCAACAGACGAATCGTGCCCGGATCAAATTTCTGTGCGGCATCCTGTTCATTGCCGATGTTGTACTGATTAACt                                                   <  4:69160/148‑1 (MQ=255)
                                                                                                  aGTACCACCTTCATATGTAATAATTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCTTGTTACTCAAAAGACGAATCGTGCCCGGATCAAATTTCTGTGCGGCATCATGTTCATTGCCGATGTTGTACTGATTAACTTGc                                                <  7:135432/148‑1 (MQ=255)
                                                                                                  aGTACCACCTGCATAAGGAATAATTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGTTACTCAACAGACGAATCGTGCCCGGATCAAATTTCTGTGCGGCATCCTGTTCATTGCCGATGTTGTACTGATTAACTTGc                                                <  1:183486/148‑1 (MQ=255)
                                                                                                   gTACCACCTGCATAAGGAATAATTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGTTACTCAACAGACGAATCGTGCCCGGATCAAATTTCTGTGCGGCATCCTGTTCATTGCCGATGTTGTACTGATTAACTTGCa                                               <  1:410419/148‑1 (MQ=255)
                                                                                                    tACCACCTGCATTAGGAATAATTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGTTACTCAACAGACGAATCGTGCCCGGATCAAATTTCTGTGCGGCATCCTGTTCATTGCCGATGTTGTACTGATTAACTTGCAt                                              >  8:365029/1‑148 (MQ=255)
                                                                                                       caccTGCATAAGGAATAATTTGATAGGAAACCGGATTGCCCATGTGGTTCTCTTTGTTCGGGTTACTCAACAGACGAATCGTGCCCGGATCAAATTTCTGTGCGGCATCCTGTTCATTGCCGATGTTGTACTGATTAACTTGCATGGt                                           <  4:179669/148‑1 (MQ=255)
                                                                                                        accTGCATAAGGAATAATTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGTTACTCAACAGACGAATCGTGCCCGGATCAAATTTCTGTGCGGCATCCTGTTCATTGCCGATGTTGTACTGATTAACTTGCATGGTa                                          <  1:365319/148‑1 (MQ=255)
                                                                                                            gCATAAGGAATAATTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGTTACTCAACAGACGAATCGTGCCCGGATCAAATTTCTGTGCGGCATCCTGTTCATTGCCGATGTTGTACTGATTaa                                                     >  7:96221/1‑133 (MQ=255)
                                                                                                            gCATAAGGAATAATTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGTTACTCAACAGACGAATCGTGCCCGGATCAAATTTCTGTGCGGCATCCTGTTCATTGCCGATGTTGTACTGATTaa                                                     <  8:96221/133‑1 (MQ=255)
                                                                                                                aaGGAATAATTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGTTACTCAACAGACGAATCGTGCCCGGATCAAATTTCTGTGCGGCATcc                                                                                 <  6:457713/101‑1 (MQ=255)
                                                                                                                aaGGAATAATTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGTTACTCAACAGACGAATCGTGCCCGGATCAAATTTCTGTGCGGCATcc                                                                                 >  5:457713/1‑101 (MQ=255)
                                                                                                                aaGGAATAATTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGTTACTCAACAGACGAATCGTGCCCGGATCAAATTTCTGTGCGGCATCCTGTTCATTGCCGATGTTGTACTGATTAACTTGc                                                <  4:254513/134‑1 (MQ=255)
                                                                                                                aaGGAATAATTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGTTACTCAACAGACGAATCGTGCCCGGATCAAATTTCTGTGCGGCATCCTGTTCATTGCCGATGTTGTACTGATTAACTTGc                                                >  3:254513/1‑134 (MQ=255)
                                                                                                                 aGGAATAATTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGTTACTCAACAGACGAATCGTGCCCGGATCAAATTTCTGTGCGGCATCCTGTTCATTGCCGATGTTGTACTGATTAACTTGCATGGTACTGGTGCg                                  <  7:327704/147‑1 (MQ=255)
                                                                                                                  ggAATAATTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGTTACTCAACAGACGAATCGTGCCCGGATCAAATTTCTGTGCGGCATCCTGTTCATTGCCGATGTTGTACTGATTAACTTGCATGGTACTGGTGCGTg                                >  1:448668/1‑148 (MQ=255)
                                                                                                                  ggAATAATTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGTTACTCAACAGACGAATCGTGCCCGGATCAAATTTCTGTGCGGCATCCTCTTCATTGCCGATGTTGTACTGATTAACTTGCATGGTACTGGTGCGTg                                <  8:370089/148‑1 (MQ=255)
                                                                                                                       aatTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGTTACTCAACAGACGAATCGTGCCCGGATCAAATTTCTGTGCGGCATCCTGTTc                                                                            >  3:7986/1‑99 (MQ=255)
                                                                                                                       aatTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGTTACTCAACAGACGAATCGTGCCCGGATCAAATTTCTGTGCGGCATCCTGTTc                                                                            <  4:7986/99‑1 (MQ=255)
                                                                                                                             aTAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGTTACTCAACAGACGAATCGTGCCCGGATCAAATTTCTGTGCGGCATCCTGTTCATTGCCGATGTTGTACTGATTAACTTGCATGGTACTGGTGCGTGGGCCACCGGc                      <  3:8751/147‑1 (MQ=255)
                                                                                                                               aGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGTTACTCAACAGACGAATCGTGCCCGGATCAAATTTCTGTGCGGCATCCTGTTCATTGCCGATGTTGTACTGATTAACTTGCATGGt                                           <  4:173978/124‑1 (MQ=255)
                                                                                                                               aGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGTTACTCAACAGACGAATCGTGCCCGGATCAAATTTCTGTGCGGCATCCTGTTCATTGCCGATGTTGTACTGATTAACTTGCATGGt                                           >  3:173978/1‑124 (MQ=255)
                                                                                                                                       ggATTGCCCATGCGGTTCTCTTTGTTCGGGTTACTCAACAGACGAATCGTGCCCGGATCAAATTTCTGTGCGGCATCCTGTTCATTGCCGATGTTGTACTGATTAACTTGCATGGTACTGGTGCGTGGGCCACCGGCAGTATTCGGtt           >  5:377194/1‑148 (MQ=255)
                                                                                                                                         aTTGCCCATGCGGTTCTCTTTGTTCGGGTTACTCAACAGACGAATCGTGCCCGGATCAAATTTCTGTGCGGCATCCTGTTCATTGCAGATGTTGTACTGATTAACTTGCATGGTACTGGTGCGTGGGCCACCGGCAGTATTCGGttt          >  2:332692/1‑147 (MQ=255)
                                                                                                                                             cccATGCGGTTCTCTTTGTTCGGGTTACTCAACAGACGAATCGTGCCCGGATCAAATTTCTGTGCGGCATCCTGTTCATTGCCGATGTTGTACTGATTACCTTGCATGGTACTGGTGCGTGGGCCACCGGCAGTATTCGGTTTTAcc      <  1:130545/147‑1 (MQ=255)
                                                                                                                                              ccATGCGGTTCTCTTTGTTCGGGTTACTCAACAGACGAATCGTGCCCGGATCAAATTTCTGTGCGGCATCCTGTTCATTGCCGATGTTGTACTGATTAACTTGCATGGTACTGGTGCGTGGGCCACCGGCAGTATTCGGTTTTACCAc    >  8:388529/1‑148 (MQ=255)
                                                                                                                                               cATGCGGTTCTCTTTGTTCGGGTTACTCAACAGACGAATCGTGCCCGGATCAAATTTCTGTGCGGCATCCTGTTCATTGCCGATGTTGTACTGATTAACTTGCATGGTACTGGTGCGTGGGCCACCGGCAGTATTCGGTTTTACCACt   >  6:449108/1‑148 (MQ=255)
                                                                                                                                                aTGCGGTTCTCTTTGTTCGGGTTACTCAACAGACGAATCGTGCCCGGATCAAATTTCTGTGCGGCATCCTGTTCATTGCCGATGTTGTACTGATTAACTTGCATGGTACTGGTGCGTGGGCCACCGGCAGTATTCGGTTTTACCACTg  <  2:403777/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GGATGATAACGCGTTACCCAGAGCTGCTTGTCCATAAAGCTTAAACGATGATAGATCCACTCGTCCGGCGCGAACTGGGCACCTTTTGCTACCGGGTGAGTACCACCTGCATAAGGAATAATTTGATAGGAAACCGGATTGCCCATGCGGTTCTCTTTGTTCGGGTTACTCAACAGACGAATCGTGCCCGGATCAAATTTCTGTGCGGCATCCTGTTCATTGCCGATGTTGTACTGATTAACTTGCATGGTACTGGTGCGTGGGCCACCGGCAGTATTCGGTTTTACCACTG  >  NC_000913/1449326‑1449617

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