Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,451,620 C→T 6.2% P51S (CCG→TCG)  lnt → apolipoprotein N‑acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,451,6200CT 6.2% 244.9 / 3.6 82P51S (CCG→TCG) lntapolipoprotein N‑acyltransferase
Reads supporting (aligned to +/- strand):  ref base C (37/39);  new base T (3/2);  total (40/42)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.92e-01

GCGTTGGATCACCCGCCCCGGCTGGCCCGGTAACCTGCTGGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTACGGCGGCGCCTCGCCGCTGCTGGCAATCGTGTTGCTGCTGGCGTTCTTCGCCGCACT  >  NC_002947/5451472‑5451765
                                                                                                                                                     |                                                                                                                                                 
gCGTTGGATCACCCGCCCCGGCTGGCCCGGTAACCTGCTGGCCCTGGCGGCCGGTCCTTCCACCCTCCTGGCCCAGGCGCCGTTCCACAACTGGACGTTGGCGGTGTTCTCCATCGAGGTGCACTAACTTGGACTGCGCGAGCACAccc                                                                                                                                                   >  2:227230/1‑149 (MQ=255)
      gATCACCCGCCCCGGCTGGCCCGGTAACCTGCTGGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCcgcgag                                                                                                                                             >  3:72481/1‑147 (MQ=255)
      gATCACCCGCCCCGGCTGGCCCGGTAACCTGCTGGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCcgc                                                                                                                                                <  4:72481/146‑1 (MQ=255)
        tCACCCGCCCCGGCTGGCCCGGTAACCTGCTGGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAg                                                                                                                                           <  2:312318/149‑1 (MQ=255)
          aCCCGCCCCGGCTGGCCCGGTAACCTGCTGGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGGCGCGGGCTGAGCTCgcgcagg                                                                                                                                           <  2:61358/148‑4 (MQ=255)
           cccGCCCCGGCTGGCCCGGTAACCTGCTGGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGc                                                                                                                                         >  2:102298/1‑148 (MQ=255)
           cccGCCCCGGCTGGCCCGGTAACCTGCTGGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGGTCGACATCTGGCCGTTGGCGGTGGTGTCCCTCGCGGTGAGCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGcc                                                                                                                                        >  4:263113/1‑149 (MQ=255)
            ccGCCCCGGCTGGCCCGGTAACCTGCTGGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACGTCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCa                                                                                                                                       >  8:72401/1‑149 (MQ=255)
            ccGCCCCGGCTGGCCCGGTAACCTGCTGGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCa                                                                                                                                       >  6:56325/1‑149 (MQ=255)
             cGCCCCGGCTGGCCCGGTAACCTGCTGGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGGCGCGGGCTGAGCTCgcgcagg                                                                                                                                           >  2:150349/1‑142 (MQ=255)
             cGCCCCGGCTGGCCCGGTAACCTGCTGGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGGCGCGGGCTGAGCTCgcgcagg                                                                                                                                           <  1:150349/145‑4 (MQ=255)
               ccccGGCTGGCCCGGTAACCTGCTGGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCAtgt                                                                                                                                    <  7:34456/149‑1 (MQ=255)
                    gCTGGCCCGGTAACCTGCTGGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATgtggcgt                                                                                                                               >  5:298441/1‑149 (MQ=255)
                          ccGGTAACCTGCTGGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCtgg                                                                                                                         >  4:66338/1‑149 (MQ=255)
                          ccGGTAACCTGCTGGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCtgg                                                                                                                         <  5:111701/149‑1 (MQ=255)
                          ccGGTAACCTGCTGGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCtgg                                                                                                                         <  2:78672/149‑1 (MQ=255)
                           cGGTAACCTGCTGGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCtggt                                                                                                                        >  7:354482/1‑149 (MQ=255)
                           cGGTAACCTGCTGGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCtggt                                                                                                                        >  5:281491/1‑149 (MQ=255)
                            ggTAACCTGCTGGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCttgtg                                                                                                                       <  1:387361/149‑1 (MQ=255)
                                 ccTGCTGGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTggtttg                                                                                                                  <  8:325533/149‑1 (MQ=255)
                                 ccTGCTGGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAACCCGCGCCAGGCCATGTGGCGTGGCTGGTggtttg                                                                                                                  >  4:74850/1‑149 (MQ=255)
                                  ctgctgGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTggtttgg                                                                                                                 >  2:135010/1‑149 (MQ=255)
                                  ctgctgGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTggtttgg                                                                                                                 <  3:320805/149‑1 (MQ=255)
                                     ctgGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTggtttgggtt                                                                                                              <  1:102298/149‑1 (MQ=255)
                                     ctgGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTggtttggg                                                                                                                <  5:313587/147‑1 (MQ=255)
                                        gCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACGTCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTggtttgggttt                                                                                                             <  7:72401/147‑1 (MQ=255)
                                         cccTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGc                                                                                                          <  6:281491/149‑1 (MQ=255)
                                              gcggcCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGt                                                                                                      <  1:227230/148‑1 (MQ=255)
                                                ggcCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGCGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATg                                                                                                   <  1:106614/149‑1 (MQ=255)
                                                    ggTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGgcgc                                                                                               >  1:371241/1‑149 (MQ=255)
                                                    ggTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTCGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGgcgc                                                                                               <  6:276231/149‑1 (MQ=255)
                                                     gTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGcc                                                                                              <  8:354482/149‑1 (MQ=255)
                                                     gTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGcc                                                                                              >  5:226902/1‑149 (MQ=255)
                                                             aCCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCtg                                                                                      >  6:300994/1‑149 (MQ=255)
                                                              cccTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCtgg                                                                                     >  7:170449/1‑149 (MQ=255)
                                                               cctcctGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGGAGCtggt                                                                                    >  2:319511/1‑149 (MQ=255)
                                                               cctcctGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCtggt                                                                                    <  5:56325/149‑1 (MQ=255)
                                                                 tcctGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCtggtg                                                                                   >  7:333808/1‑148 (MQ=255)
                                                                 tcctGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCtggtg                                                                                   >  7:29482/1‑148 (MQ=255)
                                                                    tggccctggcGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATc                                                                               >  5:58190/1‑149 (MQ=255)
                                                                      gccctggcGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTa                                                                             >  4:253187/1‑149 (MQ=255)
                                                                       ccctggcGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTAc                                                                            <  4:43602/149‑1 (MQ=255)
                                                                       ccctggcGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTAc                                                                            <  6:226902/149‑1 (MQ=255)
                                                                           ggcggggTTCGAGTTCTGGCTGTTGGCGGTGTTGTGCTTCGCGGTGGTCTACGTTGGCCTGAGCGAGGTGAGCGCGGTGCAGGCCATGTGGCGTGGATGGTGGTTTGGGTTTGTGCTGTATTTCGCCGGCACCTGGTGGATCTACGTCa                                                                        <  8:333808/143‑1 (MQ=255)
                                                                            gcgcCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAg                                                                       <  6:298441/149‑1 (MQ=255)
                                                                              gcCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCa                                                                     <  8:170449/149‑1 (MQ=255)
                                                                              gcCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAACCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGACGCCGGCACCTGGTGGATCTACGTCAGCa                                                                     <  3:74850/149‑1 (MQ=255)
                                                                                 gTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTTGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATg                                                                   >  6:135148/1‑148 (MQ=255)
                                                                                    cGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAac                                                                <  3:66338/148‑1 (MQ=255)
                                                                                          cTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCAc                                                                                         <  7:164894/117‑1 (MQ=255)
                                                                                          cTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCAc                                                                                         >  8:164894/1‑117 (MQ=255)
                                                                                          cTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTAcg                                                         >  8:204415/1‑149 (MQ=255)
                                                                                               cGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTACggcggc                                                    <  5:338413/149‑1 (MQ=255)
                                                                                                            gtCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTACGGCGGCGCCTCGCCgctgc                                       >  6:98469/1‑149 (MQ=255)
                                                                                                             ctcATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTACGGCGGCGCCTCGCCgctgct                                      <  7:331876/147‑1 (MQ=255)
                                                                                                               cATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTACGGCGGCGCCTCGCCGCTGCTgg                                    <  7:204415/149‑1 (MQ=255)
                                                                                                               cATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTACGGCGGCGCCTCGCCGCTGCTgg                                    <  2:371241/149‑1 (MQ=255)
                                                                                                                aTCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTACGGCGGCGCCTCGCCgctgct                                      >  8:331876/1‑146 (MQ=255)
                                                                                                                aTCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCAGGCACCTGGTGGATCTACGTCAGCATGAACACCTACGGCGGCGCCTCGCCGCTGCTgg                                    >  8:134540/1‑148 (MQ=255)
                                                                                                                 tCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTACGGCGGCGCCTCGCCGCTGCTGGc                                   >  3:58075/1‑148 (MQ=255)
                                                                                                                    cgGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTACGGCGGCGCCTCGCCGCTGCTGGCAATc                               >  5:380430/1‑149 (MQ=255)
                                                                                                                    cgGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTACGGCGGCGCCTCGCCGCTGCTGGCAATc                               >  3:12669/1‑149 (MQ=255)
                                                                                                                    cgGAGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTACGGCGGCGCCTCGCCGCTGCTGGCAATc                               <  4:110090/149‑1 (MQ=255)
                                                                                                                       tGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTACGGCGGCGCCTCGCCGCTGCTGGCAATCgtg                            >  8:267467/1‑149 (MQ=255)
                                                                                                                         ctctACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTACGGCGGCGCCTCGCCGCTGCTGGCAATCGTGtt                          >  3:59203/1‑149 (MQ=255)
                                                                                                                             cacatGGCCTGGCGCGGGCTGAGCTCGCGCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTACGGCGGCGCCTCGCCGCTGCTGGCAATCGTGTtgctg                      >  6:18002/5‑149 (MQ=38)
                                                                                                                              ccTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCa                                                                        <  6:329987/98‑1 (MQ=255)
                                                                                                                              ccTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCa                                                                        >  5:329987/1‑98 (MQ=255)
                                                                                                                              ccTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTACGGCGGCGCCTCGCCGCTGCTGGCAATCGTGTtgctgc                     >  4:326383/1‑149 (MQ=255)
                                                                                                                               cTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTACGGCGGCGCCTCGCCGCTGCTGGCAATCGTGTTgctgct                    <  6:58190/149‑1 (MQ=255)
                                                                                                                                 tGGCCTGCGCGAGCTCAGCCCGCGCCGGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTACGGCGGCGCCTCGCCGCTGCTGGCAATCGTGTTGCTGCTgg                  >  1:208463/1‑149 (MQ=255)
                                                                                                                                   gCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTACGGCGGCGCCTCGCCGCTGCTGGCAATCGTGTTGCTGCTGGc                 <  6:380430/148‑1 (MQ=255)
                                                                                                                                   gCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTACGGCGGCGCCTCGCCGCTGCTGGCAATCGTGTTGCTGCTGGCg                <  5:135148/149‑1 (MQ=255)
                                                                                                                                    cctggcgcGGGCTGAGCTCGCGCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTACGGCGGCGCCTCGCCGCTGCTGGCAATCGTGTTGCTGCTGGCGt               >  5:96653/5‑149 (MQ=255)
                                                                                                                                    ccTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTACGGCGGCGCCTCGCCGCTGCTGGCAATCGTGTTGCTGCTGGCGt               <  3:252506/149‑1 (MQ=255)
                                                                                                                                     cTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTACGGCGGCGCCTCGCCGCTGCTGGCAATCGTGTTGCTGCTGGCGtt              >  8:160225/1‑149 (MQ=255)
                                                                                                                                         cgcgAGCTCTGTCCGCTCGAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTTGTGGATCTACGTCAGCATGAACACCTACGGCGGGGCCTCGCCGCTGCTGGCAATCGTGTTGCTGCTGGCGttc             <  7:282324/147‑1 (MQ=255)
                                                                                                                                          gcgAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTACGGCGGCGCCTCGCCGCTGCTGGCAATCGTGTTGCTGCTGGCGttcttc          <  4:12669/149‑1 (MQ=255)
                                                                                                                                               cTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTACGGCGGCGCCTCGCCGCTGCTGGCAATCGTGtt                          >  6:83912/1‑128 (MQ=255)
                                                                                                                                               cTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTACGGCGGCGCCTCGCCGCTGCTGGCAATCGTGtt                          <  5:83912/128‑1 (MQ=255)
                                                                                                                                                tCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTACGGCGGCGCCTCGCCGCTGCTGGCAATCGTGTTGCTGCTGGCGTTCTTcgcc       <  3:300583/146‑1 (MQ=255)
                                                                                                                                                tCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTACGGCGGCGCCTCGCCGCTGCTGGCAATCGTGTTGCTGCTGGCGTTCTTcgcc       >  4:300583/1‑146 (MQ=255)
                                                                                                                                                  aGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTACGGCGGCGCCTCGCCGCTGCTGGCAATCGTGTTGCTGCTGGCGTTCTTCGCCGCACt  <  5:98469/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GCGTTGGATCACCCGCCCCGGCTGGCCCGGTAACCTGCTGGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTACGGCGGCGCCTCGCCGCTGCTGGCAATCGTGTTGCTGCTGGCGTTCTTCGCCGCACT  >  NC_002947/5451472‑5451765

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: