Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,451,660 A→G 6.1% E48G (GAG→GGG)  lnt → apolipoprotein N‑acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,451,6600AG 6.1% 244.7 / 4.3 84E48G (GAG→GGG) lntapolipoprotein N‑acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (43/35);  new base G (3/2);  total (46/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.40e-01

ACATGCGTTGGATCACCCGCCCCGGCTGGCCCGGTAACCTGCTGGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTACGGCGGCGCCTCGCCGCTGCTGGCAATCGTGTTGCTGCTGGCGTTCTTCGCCG  >  NC_002947/5451516‑5451809
                                                                                                                                                 |                                                                                                                                                     
aCATGCGTTGGATCACCCGCCCCGGCTGGCCCGGTAACCTGCTGGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGACTGCGCGAGCTCa                                                                                                                                                  <  1:1195691/150‑1 (MQ=255)
       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                                                                                                                               >  1:475693/1‑116 (MQ=37)
                                                 ggcggcCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGGCGCGGGCTGAGCTCGCGCAGGCCaaggt                                                                                                                               <  2:475693/120‑5 (MQ=37)
                                                 ggcggcCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTa                                                                                                 <  2:689513/150‑1 (MQ=255)
                                                        ggTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGcc                                                                                          >  1:330472/1‑150 (MQ=255)
                                                            cTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTggtttggg                                                                                                            <  2:169925/128‑1 (MQ=255)
                                                            cTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTggtttggg                                                                                                            >  1:169925/1‑128 (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                                                                                <  1:1039119/150‑1 (MQ=255)
                                                                  cccTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCtggt                                                                                <  2:1196637/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                                                                                                      <  1:472976/117‑1 (MQ=255)
                                                                             ctggcGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGc                                                                                                      >  2:472976/1‑117 (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                                                   <  2:599131/126‑1 (MQ=255)
                                                                                                                       gcgGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTAcggc                                                   >  1:599131/1‑126 (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)
                                                                                                                                                 |                                                                                                                                                     
ACATGCGTTGGATCACCCGCCCCGGCTGGCCCGGTAACCTGCTGGCCCTGGCGGCCGGTGCTTCCACCCTCCTGGCCCTGGCGCCGTTCGACATCTGGCCGTTGGCGGTGTTGTCCATCGCGGTGCTCTACCTTGGCCTGCGCGAGCTCAGCCCGCGCCAGGCCATGTGGCGTGGCTGGTGGTTTGGGTTTGGCCTGTATGGCGCCGGCACCTGGTGGATCTACGTCAGCATGAACACCTACGGCGGCGCCTCGCCGCTGCTGGCAATCGTGTTGCTGCTGGCGTTCTTCGCCG  >  NC_002947/5451516‑5451809

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