Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,451,668 C→T 5.7% 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,6680CT 5.7% 258.9 / 4.4 87P51S (CCG→TCG) lntapolipoprotein N‑acyltransferase
Reads supporting (aligned to +/- strand):  ref base C (43/39);  new base T (3/2);  total (46/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

TTGGATCACCCGCCCCGGCTGGCCCGGTAACCTGCTGGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTACGGCGGCGCCTCGCCGCTGCTGGCAATCGTGTTGCTGCTGGCGTTCTTCGCCGCACTG  >  NC_002947/5451523‑5451814
                                                                                                                                                  |                                                                                                                                                  
ttGGATCACCCGCCCCGGCTGGCCCGGTAACCTGCTGGCCCTGGCGGCCGGTGCTTCCACCATCCTGGCCCTGGGTGCGTTCGACATCTGGGCGTTGGGTGTTTGGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCcgcg                                                                                                                                                >  1:571042/1‑150 (MQ=255)
 tGGATCACCCGCCCCGGCTGGCCCGGTAACCTGCTGGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCgcgc                                                                                                                                               <  2:478365/150‑1 (MQ=255)
   gATCACCCGCCCCGGCTGGCCCGGTAACCTGCTGGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTCCCTTGGCCTGCGCGAGCTCAGCCCgcgcag                                                                                                                                             <  1:630545/150‑3 (MQ=255)
   gATCACCCGCCCCGGCTGGCCCGGTAACCTGCTGGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGGGCTGAGCTCgcgcag                                                                                                                                             <  2:93255/150‑3 (MQ=255)
   gATCACCCGCCCCGGCTGGCCCGGTAACCTGCTGGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCa                                                                                                                                             >  2:301662/1‑150 (MQ=255)
    aTCACCCGCCCCGGCTGGCCCGGTAACCTGCTGGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAg                                                                                                                                            >  2:846580/1‑150 (MQ=255)
     tCACCCGCCCCGGCTGGCCCGGTAACCTGCTGGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAgg                                                                                                                                           >  2:290767/1‑150 (MQ=255)
     tCACCCGCCCCGGCTGGCCCGGTAACCTGCTGGCCCTGGCGGCCGGCGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAgg                                                                                                                                           >  2:421051/1‑150 (MQ=255)
        cccGCCCCGGCTGGCCCGGTAACCTGCTGGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCa                                                                                                                                        >  1:787391/1‑150 (MQ=255)
        cccGCCCCGGCTGGCCCGGTAACCTGCTGGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCa                                                                                                                                        >  2:934682/1‑150 (MQ=255)
         ccGCCCCGGCTGGCCCGGTAACCTGCTGGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCAt                                                                                                                                       >  1:1179797/1‑150 (MQ=255)
         ccGCCCCGGCTGGCCCGGTAACCTGCTGGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCAt                                                                                                                                       >  1:392166/1‑150 (MQ=255)
         ccGCCCCGGCTGGCCCGGTAACCTGCTGGCCCTGGCGGCCGGTGCGTCCACCCTCCTGGCCCTGCCGCCGTTCGCCAGCTGGCCGTCGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTCCGCGAGCTCCGCCCGCGCCAGGCCac                                                                                                                                       >  2:490580/1‑149 (MQ=255)
            ccccGGCTGGCCCGGTAACCTGCTGGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCAtgtg                                                                                                                                    <  2:490194/150‑1 (MQ=255)
              ccGGCTGGCCCGGTAACCTGCTGGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATgtggc                                                                                                                                  >  2:33920/1‑150 (MQ=255)
                  cTGGCCCGGTAACCTGCTGGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATgtggcgtgg                                                                                                                              <  2:594241/150‑1 (MQ=255)
                       ccGGTAACCTGCTGGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCtggt                                                                                                                         >  1:689513/1‑150 (MQ=255)
                        cGGTAACCTGCTGGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCtggtg                                                                                                                        <  1:1003320/150‑1 (MQ=255)
                        cGGTAACCTGCTGGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCtggtg                                                                                                                        >  1:997958/1‑150 (MQ=255)
                                tgctgGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCtggtg                                                                                                                        >  2:1171415/1‑142 (MQ=255)
                                tgctgGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCtggtg                                                                                                                        <  1:1171415/142‑1 (MQ=255)
                                tgctgGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTggtttgggt                                                                                                                >  2:1153976/1‑150 (MQ=255)
                                 gctgGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTggtttgggtt                                                                                                               >  2:1071075/1‑150 (MQ=255)
                                  ctgGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTggtttgggttt                                                                                                              >  2:883044/1‑150 (MQ=255)
                                   tgGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTggtttgggtttg                                                                                                             <  1:934682/150‑1 (MQ=255)
                                      cccTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGcc                                                                                                          >  2:9374/1‑150 (MQ=255)
                                       ccTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCt                                                                                                         <  1:1071075/150‑1 (MQ=255)
                                          ggcggcCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGGCGCGGGCTGAGCTCGCGCAGGCCaaggt                                                                                                                                    <  2:475693/120‑5 (MQ=37)
                                          ggcggcCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGGCGCGGGCTGAGCTCGCGCAGGCCaaggt                                                                                                                                    >  1:475693/1‑116 (MQ=37)
                                          ggcggcCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTa                                                                                                      <  2:689513/150‑1 (MQ=255)
                                                 ggTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGcc                                                                                               >  1:330472/1‑150 (MQ=255)
                                                     cTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTggtttggg                                                                                                                 >  1:169925/1‑128 (MQ=255)
                                                     cTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTggtttggg                                                                                                                 <  2:169925/128‑1 (MQ=255)
                                                     cTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCa                                                                                           >  2:426759/1‑150 (MQ=255)
                                                     cTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCa                                                                                           <  1:9374/150‑1 (MQ=255)
                                                      ttCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCAc                                                                                          >  1:242351/1‑150 (MQ=255)
                                                         cACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCtg                                                                                       <  1:1180918/150‑1 (MQ=255)
                                                          aCCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGAGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGCCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCtgg                                                                                      >  1:1147433/1‑150 (MQ=255)
                                                           cccTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCtggt                                                                                     <  2:1196637/150‑1 (MQ=255)
                                                           cccTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCtggt                                                                                     <  1:1039119/150‑1 (MQ=255)
                                                            cctcctGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCtggtg                                                                                    <  1:301662/150‑1 (MQ=255)
                                                            cctcctGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCtggtg                                                                                    >  1:464170/1‑150 (MQ=255)
                                                              tcctGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGt                                                                                                       >  1:649793/1‑129 (MQ=255)
                                                              tcctGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGt                                                                                                       <  2:649793/129‑1 (MQ=255)
                                                              tcctGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGa                                                                                  >  2:820961/1‑150 (MQ=255)
                                                                ctggccctggcGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATc                                                                                <  1:421051/150‑1 (MQ=255)
                                                                ctggccctggcGCCGTTCGAAATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCtggtg                                                                                    <  2:110106/146‑1 (MQ=255)
                                                                ctggccctggcGCCGTTCGAAATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCtggtg                                                                                    >  1:110106/1‑146 (MQ=255)
                                                                 tggccctggcGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCtggtg                                                                                                                        <  2:110040/109‑1 (MQ=255)
                                                                 tggccctggcGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCtggtg                                                                                                                        >  1:110040/1‑109 (MQ=255)
                                                                   gccctggcGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTAc                                                                             <  1:303506/150‑1 (MQ=255)
                                                                      ctggcGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGc                                                                                                           >  2:472976/1‑117 (MQ=255)
                                                                      ctggcGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGc                                                                                                           <  1:472976/117‑1 (MQ=255)
                                                                       tggcGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCa                                                                         >  2:76651/1‑150 (MQ=255)
                                                                            ccGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATg                                                                    <  2:997958/150‑1 (MQ=255)
                                                                             cGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGa                                                                   >  1:177894/1‑150 (MQ=255)
                                                                                    cATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCt                                                                                                         <  2:320570/105‑1 (MQ=255)
                                                                                    cATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCt                                                                                                         >  1:320570/1‑105 (MQ=255)
                                                                                     aTCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGcc                                                                                               >  2:786025/1‑114 (MQ=255)
                                                                                     aTCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGcc                                                                                               <  1:786025/114‑1 (MQ=255)
                                                                                        tGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTAcggc                                                        >  2:762847/1‑150 (MQ=255)
                                                                                          gCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTAcggcgg                                                      <  2:571042/150‑1 (MQ=255)
                                                                                           ccGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGTGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTACggcggc                                                     <  2:31679/150‑1 (MQ=255)
                                                                                             gTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCTCCTGGTGGATCTACGTCAGCATGAACACCTACGGCGgcgc                                                   <  1:1153976/150‑1 (MQ=255)
                                                                                             gTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGATGGATTGGGGTTGCCCTGTGTGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTGCGGCGgcgc                                                   >  1:574150/1‑150 (MQ=255)
                                                                                                ggCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTACGGCGGCGCCTc                                                <  1:33920/150‑1 (MQ=255)
                                                                                                       ttgtCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGCTCTACTTCAGCATGAACACCTACGGCGGCGCCTCGCCgctg                                         >  2:1189106/1‑150 (MQ=255)
                                                                                                          tCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTggtttgg                                                                                                                  >  1:553782/1‑74 (MQ=255)
                                                                                                          tCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTggtttgg                                                                                                                  <  2:553782/74‑1 (MQ=255)
                                                                                                            cATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTACGGCGGCGCCTCGCCGCTGCTGGc                                    >  2:328197/1‑150 (MQ=255)
                                                                                                                gcgGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTAcggc                                                        >  1:599131/1‑126 (MQ=255)
                                                                                                                gcgGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTAcggc                                                        <  2:599131/126‑1 (MQ=255)
                                                                                                                   gTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTCTAGCGCCGGCACCTAGTGGAACTACTTGAGCAGCAACACCTACGGCGGCGCCTCGCCGCTGCTGGCAATCgtg                             <  1:328197/150‑1 (MQ=255)
                                                                                                                      ctctACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTACGGCGGCGCCTCGCCGCTGCTGGCAATCGTGTtg                          >  2:612446/1‑150 (MQ=255)
                                                                                                                         tACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTACGGCGGCGCCTCGCCGCTGCTGGCAATCGTGTtgctg                       >  2:548220/1‑150 (MQ=255)
                                                                                                                         tACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTACGGCGGCGCCTCGCCGCTGCTGGCAATCGTGTtgctg                       >  2:485286/1‑150 (MQ=255)
                                                                                                                            cTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTACGGCGGCGCCTCGCCGCTGCTGGCAATCGTGTTGctgctg                    >  1:1018452/1‑150 (MQ=255)
                                                                                                                              tGGCTTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTACGGCGGCGCCTTGCCGCTGCTGGCAATCGTGTTGCTGCTGGc                  <  2:491391/150‑1 (MQ=255)
                                                                                                                                    ggcgcGGGCTGAGCTCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTACGGCGGCGCCTCGCCGCTGCTGGCAATCGTGTTGCTGCTGGCGttct             >  2:129684/2‑150 (MQ=255)
                                                                                                                                     gcgcGGGCTGAGCTCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTACGGCGGCGCCTCGCCGCTGCTGGCAATCGTGTTGCTGCTGGCGttctt            >  2:977728/1‑150 (MQ=255)
                                                                                                                                      cgcgAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTACGGCGGCGCCTCGCCGCTGCTGGCAATCGTGTTGCTGCTGGCGttcttc           <  1:426759/150‑1 (MQ=255)
                                                                                                                                        cgAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTACGGCGGCGCCTCGCCGCTGCTGGCAATCGTGTTGCTGCTGGCGTTCTTcgc         <  2:242351/150‑1 (MQ=255)
                                                                                                                                         gAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTACGGCGGCGCCTCGCCGCTGCTGGCAATCGTGTTGCTGCTGGCGTTCTTcgcc        <  2:1018452/150‑1 (MQ=255)
                                                                                                                                          aGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTACGGCGGCGCCTCGCCGCTGCTGGCAATCGTGTTGCTGCTGGCGTTCTTcgccg       <  1:883044/150‑1 (MQ=255)
                                                                                                                                               aGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTACGGCGGCGCCTCGCCGCTGCTGGCAATCGTGTTGCTGCTGGCGTTCTTCGCCGCACTg  <  1:846580/150‑1 (MQ=255)
                                                                                                                                               aGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTACGGCGGCGCCTCGCCGCTGCTGGCAATCGTGTTGCTGCTGGCGTTCTTCGCCGCACTg  <  2:464170/150‑1 (MQ=255)
                                                                                                                                               aGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGACAGCACCTGGTGGATCTACGTCAGCATGAACACCTACGGCGGCGCCTCGCCGCTGCTGGCAATCGTGTTGCTGCTGGCGTTCTTCGCCGCACTg  <  1:317816/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TTGGATCACCCGCCCCGGCTGGCCCGGTAACCTGCTGGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTACGGCGGCGCCTCGCCGCTGCTGGCAATCGTGTTGCTGCTGGCGTTCTTCGCCGCACTG  >  NC_002947/5451523‑5451814

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