Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 990,598 G→T 5.1% G160V (GGC→GTC)  PP_0856 → lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947990,5980GT 5.1% 258.1 / 2.7 79G160V (GGC→GTC) PP_0856lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (40/35);  new base T (2/2);  total (42/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

CCTCAAGGGTGAAGTCATTGCCCTGGACTCCAGCACCGGTGAGGAGCGTTGGCGCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGT‑‑GCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCGGCAAGGTAGTGGCGGTGGACATCAA  >  NC_002947/990458‑990745
                                                                                                                                             |                                                                                                                                                     
ccTCAAGGGTGAAGTCATTGCCCTGGACTCCAGCACCGGTGAGGAGCGTTGGCGCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATg                                                                                                                                               >  3:763/1‑149 (MQ=255)
 cTCAAGGGTGAAGTCATTGCCCTGGACTCCAGCACCGGTGAGGAGCGTTGGCGCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTAGTGGTGCAGACCCAGGCCGATCAGGCGATCGTCCTgg                                                                                                                                                  >  5:342702/1‑145 (MQ=255)
 cTCAAGGGTGAAGTCATTGCCCTGGACTCCAGCACCGGTGAGGAGCGTTGGCGCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTAGTGGTGCAGACCCAGGCCGATCAGGCGATCGTCCTgg                                                                                                                                                  <  6:342702/145‑1 (MQ=255)
    aaGGGTGAAGTCATTGCCCTGGACTCCAGCACCGGTGAGGAGCGTTGGCGCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATgcc                                                                                                                                             <  5:7403/147‑1 (MQ=255)
      gggTGAAGTCATTGCCCTGGACTCCAGCACCGGTGAGGAGCGTTGGCGCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCa                                                                                                                                         <  5:16167/149‑1 (MQ=255)
       ggTGAAGTCATTGCCCTGGACTCCAGCACCGGTGAGGAGCGTTGGCGCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCt                                                                                                                                                    >  8:165528/1‑137 (MQ=255)
       ggTGAAGTCATTGCCCTGGACTCCAGCACCGGTGAGGAGCGTTGGCGCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCt                                                                                                                                                    <  7:165528/137‑1 (MQ=255)
         tGAAGTCATTGCCCTGGACTCCAGCACCGGTGAGGAGCGTTGGCGCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTg                                                                                                                                      >  2:335483/1‑149 (MQ=255)
                       tGGACTCCAGCACCGGTGAGGAGCGTTGGCGGCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTg                                                                                                                         >  1:226146/1‑149 (MQ=255)
                       tGGACTCCAGCACCGGTGAGGAGCGTTGGCGCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCAACCGCCGTTg                                                                                                                         >  1:344132/1‑148 (MQ=255)
                         gACTCCAGCACCGGTGAGGAGCGTTGGCGCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTgg                                                                                                                        <  1:358031/147‑1 (MQ=255)
                         gACTCCAGCACCGGTGAGGAGCGTTGGCGCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGa                                                                                                                       <  2:226146/148‑1 (MQ=255)
                         gACTCCAGCACCGGTGAGGAGCGTTGGCGCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCAACCGCCGTTg                                                                                                                         <  2:344132/146‑1 (MQ=255)
                           cTCCAGCACCGGTGAGGAGCGTTGGCGCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCAGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATc                                                                                                                     >  5:113975/1‑148 (MQ=255)
                            tCCAGCACCGGTGAGGAGCGTTGGCGCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTa                                                                                                                   >  2:266113/1‑149 (MQ=255)
                              cAGCACCGGTGAGGAGCGTTGGCGCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGACAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTAc                                                                                                                  <  2:127200/148‑1 (MQ=255)
                                   ccGGTGAGGAGCGTTGGCGCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAAc                                                                                                            <  2:79132/149‑1 (MQ=255)
                                    cGGTGAGGAGCGTTGGCGCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACt                                                                                                           >  2:191800/1‑149 (MQ=255)
                                           gagCGTTGGCGCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCgccggc                                                                                                    >  5:356830/1‑149 (MQ=255)
                                           gagCGTTGGCGCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCgccggc                                                                                                    >  7:262190/1‑149 (MQ=255)
                                            agCGTTGGCGCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCgccggcc                                                                                                   >  5:151394/1‑149 (MQ=255)
                                              cGTTGGCGCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGt                                                                                                 >  2:153026/1‑149 (MQ=255)
                                                   gcgcTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCgc                                                                                                        >  1:352788/1‑137 (MQ=255)
                                                   gcgcTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCgc                                                                                                        <  2:352788/137‑1 (MQ=255)
                                                      gcTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACg                                                                                          >  2:273696/1‑149 (MQ=255)
                                                      gcTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACg                                                                                          <  1:153026/149‑1 (MQ=255)
                                                        tCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGTGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCt                                                                                        >  6:84210/1‑149 (MQ=255)
                                                            gcgTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGt                                                                                    <  1:335483/149‑1 (MQ=255)
                                                            gcgTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGt                                                                                    >  3:213684/1‑149 (MQ=255)
                                                                 aCCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGCCGATCAGGCGATCGTCCTgg                                                                                                                                                  <  5:169186/82‑1 (MQ=37)
                                                                 aCCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGCCGATCAGGCGATCGTCCTgg                                                                                                                                                  >  6:169186/1‑82 (MQ=37)
                                                                  ccAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCgccgg                                                                                                     >  2:205275/1‑126 (MQ=255)
                                                                  ccAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCgccgg                                                                                                     <  1:205275/126‑1 (MQ=255)
                                                                   cAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCg                                                                             <  4:213684/149‑1 (MQ=255)
                                                                    aGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCgg                                                                            >  6:105018/1‑149 (MQ=255)
                                                                         aGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGa                                                                                            >  1:217069/1‑128 (MQ=255)
                                                                         aGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGa                                                                                            <  2:217069/128‑1 (MQ=255)
                                                                                gCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGT‑‑GCGCCg                                                                   <  5:217834/144‑1 (MQ=255)
                                                                                gCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGT‑‑GCGCCg                                                                   >  6:217834/1‑144 (MQ=255)
                                                                                gCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGT‑‑GCGCCGATCGc                                                              <  5:105018/149‑1 (MQ=255)
                                                                                 ccccACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCCCCGGT‑‑GCGCCGATCGc                                                              >  3:70492/1‑148 (MQ=255)
                                                                                   ccACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGT‑‑GCGCCGATCGccac                                                           >  6:294618/1‑149 (MQ=255)
                                                                                      ccTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGT‑‑GCGCCGATCGCCACCaa                                                        <  4:763/149‑1 (MQ=255)
                                                                                            aacaacGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGT‑‑GCGCCGATCGCCACCAACCGCCt                                                  >  2:369519/1‑149 (MQ=255)
                                                                                             acaacGGTGACGTGGTGGAGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGT‑‑GCGCCGATCGCCACCAACCGCCTg                                                 >  7:20277/1‑149 (MQ=255)
                                                                                                acGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGT‑‑GCGCCGATCGCCACCAACCGCCTGGCg                                              <  6:356830/149‑1 (MQ=255)
                                                                                                  ggTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGT‑‑GCGCCGATCGCCACCAACCGCCTGGCGGt                                            >  2:108952/1‑149 (MQ=255)
                                                                                                    tGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGT‑‑GCGCCGATCGCCACCAACCGCCTGGCGGTTg                                          >  3:23820/1‑149 (MQ=255)
                                                                                                      aCGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGT‑‑GCGCCGATCGCCACCAACCGCCTGGCGGTTGc                                         <  8:262190/148‑1 (MQ=255)
                                                                                                           gtggtggtgCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCTCCGTTGGATCTACGAATACTCGCCGGCCGTGCTGACGGTGCGTGGCACCGGT‑‑GCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCt                                   <  6:113975/149‑1 (MQ=255)
                                                                                                           gtggtggtgCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGT‑‑GCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCt                                   >  8:85392/1‑149 (MQ=255)
                                                                                                           gtggtggtgCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGCCACCGGT‑‑GCACCGATCGCCACCAACCGCCTGGCGGTTACCGGCCt                                   <  1:266113/149‑1 (MQ=255)
                                                                                                               tggtgCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGT‑‑GCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTc                                >  3:266481/1‑148 (MQ=255)
                                                                                                                 gtgCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGT‑‑GCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTAc                             >  8:293897/1‑149 (MQ=255)
                                                                                                                   gCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGT‑‑GCGCCGATCGCCACCAACCGCCAGGCGGTTggtggc                                     <  3:134369/139‑6 (MQ=255)
                                                                                                                   gCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGT‑‑GCGCCGATCGCCACCAACCGCCAGGCGGTTggtggc                                     >  4:134369/1‑134 (MQ=255)
                                                                                                                   gCAGAAACAGGACGATCGCCTGATCGGCCTGGATGCAGCCACTGGCGAACGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGT‑‑GCGCCGATCGCAACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCg                           >  8:304815/1‑149 (MQ=255)
                                                                                                                       aCCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGCCCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGTGTGGCACCGGT‑‑gc                                                                       >  4:29163/1‑101 (MQ=255)
                                                                                                                       aCCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGTGTGGCACCGGT‑‑gc                                                                       <  3:29163/101‑1 (MQ=255)
                                                                                                                         ccAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGT‑‑GCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCGGCAAgg                     >  8:155465/1‑149 (MQ=255)
                                                                                                                         ccAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGT‑‑GCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCGGCAAgg                     >  3:44274/1‑149 (MQ=255)
                                                                                                                           aGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGT‑‑GCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCgg                          <  3:82941/142‑1 (MQ=255)
                                                                                                                           aGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGT‑‑GCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCgg                          >  4:82941/1‑142 (MQ=255)
                                                                                                                           aGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGT‑‑GCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCGGCAAGGTa                   <  1:369519/149‑1 (MQ=255)
                                                                                                                             gACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGT‑‑GCGCCGATCGCCACCAACCGCCTGGCGGTTGCCgg                                      <  1:296965/128‑1 (MQ=255)
                                                                                                                             gACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGT‑‑GCGCCGATCGCCACCAACCGCCTGGCGGTTGCCgg                                      >  2:296965/1‑128 (MQ=255)
                                                                                                                             gACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGT‑‑GCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCGGCAAGGTAGt                 <  1:229805/149‑1 (MQ=255)
                                                                                                                               cGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGT‑‑GCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCGGCAAGGTAGTgg               <  1:191800/149‑1 (MQ=255)
                                                                                                                                  tCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGT‑‑GCGCCGATCGccac                                                           <  8:190108/102‑1 (MQ=255)
                                                                                                                                  tCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGT‑‑GCGCCGATCGccac                                                           >  7:190108/1‑102 (MQ=255)
                                                                                                                                  tCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGT‑‑GCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCGGCAAGGTAGTGGCgg            >  6:263638/1‑149 (MQ=255)
                                                                                                                                   cGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGGGCTGACGCTGCGTGGCACCGGT‑‑GCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCGGCAAGGTAGTGGCGGt           <  4:229580/149‑1 (MQ=255)
                                                                                                                                    gCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTCCGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCGGCAAGGTAGTGGCgg            <  1:108952/149‑1 (MQ=255)
                                                                                                                                     ccTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGT‑‑GCGCCGATCGCCACCAACCGCCTGGCGGTTGCCg                                       >  1:294953/1‑119 (MQ=255)
                                                                                                                                     ccTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGT‑‑GCGCCGATCGCCACCAACCGCCTGGCGGTTGCCg                                       <  2:294953/119‑1 (MQ=255)
                                                                                                                                       tGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGT‑‑GCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGt                                 <  6:385014/123‑1 (MQ=255)
                                                                                                                                       tGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGT‑‑GCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGt                                 >  5:385014/1‑123 (MQ=255)
                                                                                                                                         aTCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGT‑‑GCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCGGCAAGGTAGTGGCGGTGGAc       <  7:293897/147‑1 (MQ=255)
                                                                                                                                          tCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGT‑‑GCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCGGCAAGGTAGTGGCGGTGGACATc    <  7:85392/149‑1 (MQ=255)
                                                                                                                                            ggCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGT‑‑GCGCCGATCGCCACCAACCGCCTGGCGGTTGCCg                                       <  3:298885/112‑1 (MQ=255)
                                                                                                                                            ggCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGT‑‑GCGCCGATCGCCACCAACCGCCTGGCGGTTGCCg                                       >  4:298885/1‑112 (MQ=255)
                                                                                                                                             gCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGT‑‑GCGCCGATCGcc                                                             >  8:247627/1‑89 (MQ=255)
                                                                                                                                             gCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGT‑‑GCGCCGATCGcc                                                             <  7:247627/89‑1 (MQ=255)
                                                                                                                                             gCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGT‑‑GCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCGGCAAGGTAGTGGCGGTGGACATcaa  >  4:256880/1‑148 (MQ=255)
                                                                                                                                             gCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGT‑‑GCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCGGCAAGGTAGTGGCGGTGGACATcaa  >  4:56794/1‑148 (MQ=255)
                                                                                                                                             |                                                                                                                                                     
CCTCAAGGGTGAAGTCATTGCCCTGGACTCCAGCACCGGTGAGGAGCGTTGGCGCTCGCGCGTGACCAGCGAAGTGCTGGCCCCACCTGCCAACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGT‑‑GCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTGTCTACCGGCAAGGTAGTGGCGGTGGACATCAA  >  NC_002947/990458‑990745

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: