Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,065,730 T→A 6.2% I39F (ATC→TTC)  PP_0923 ← acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,065,7300TA 6.2% 242.5 / 4.8 81I39F (ATC→TTC) PP_0923acyltransferase
Reads supporting (aligned to +/- strand):  ref base T (36/40);  new base A (3/2);  total (39/42)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-01

CCAAGCAGGGGGATGTCGGTCCAGGAAACATGGTTGCTGACCCACAACATCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGCATCAGACCGCCGCCTTGAAGTGGCGGGCGTAGCG  >  NC_002947/1065582‑1065876
                                                                                                                                                    |                                                                                                                                                  
ccAAGCAGGGGGATGTCGGTCCAGGAAACATGGTTGCTGACCCACAACATCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGAt                                                                                                                                                    <  8:183084/149‑1 (MQ=255)
 cAAGCAGGGGGATGTCGGTCCAGGAAACATGGTTGCTGACCCACAACATCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACTGGCAGGTCCAGCGCTGGCGCCGTTCGATc                                                                                                                                                   <  1:28017/149‑1 (MQ=255)
 cAAGCAGGGGGATGTCGGTCCAGGAAACATGGTTGCTGACCCACAACATCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATc                                                                                                                                                   <  8:115769/149‑1 (MQ=255)
 cAAGCAGGGGGATGTCGGTCCAGGAAACATGGTTGCTGACCCACAACATCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATAAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATc                                                                                                                                                   >  7:178497/1‑149 (MQ=255)
  aaGCAGGGGGATGTCGGTCCAGGAAACATGGTTGCTGACCCACAACATCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCg                                                                                                                                                  >  5:130531/1‑149 (MQ=255)
       gggggATGTCGGTCCAGGAAACATGGTTGCTGACCCACAACATCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGcc                                                                                                                                             <  7:147141/149‑1 (MQ=255)
        ggggATGTCGGTCCAGGAAACATGGTTGCTGACCCACAACATCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCt                                                                                                                                            <  3:339582/149‑1 (MQ=255)
         gggATGTCGGTCCAGGAAACATGGTTGCTGACCCACAACATCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCtt                                                                                                                                           <  6:130531/149‑1 (MQ=255)
         gggATGTCGGTCCAGGAAACATGGTTGCTGACCCACAACATCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCtt                                                                                                                                           <  7:110098/149‑1 (MQ=255)
         gggATGTCGGTCCAGGAAACATGGTTGCTGACCCACAACATCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCtt                                                                                                                                           <  7:9598/149‑1 (MQ=255)
         gggATGTCGGTCCAGGAAACATGGTTGCTGACCCACAACATCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCtt                                                                                                                                           <  8:232097/149‑1 (MQ=255)
                  gTCCAGGAAACATGGTTACTGACCCACAACATCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGCTCGGCGCCTTGAAGccc                                                                                                                                    >  2:227894/1‑147 (MQ=255)
                   tCCAGGAAACATGGTTGCTGACCCACAACATCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAg                                                                                                                                 >  3:297003/1‑149 (MQ=255)
                        gAAACATGGTTGCTGACCCACAACATCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTc                                                                                                                            >  4:208646/1‑149 (MQ=255)
                        gAAACATGGTTGCTGACCCACAACATCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTc                                                                                                                            >  3:364449/1‑149 (MQ=255)
                          aaCATGGTTGCTGACCCACAACATCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGc                                                                                                                          <  2:137849/149‑1 (MQ=255)
                               ggTTGCTGACCCACAACATCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAg                                                                                                                     >  3:335126/1‑149 (MQ=255)
                                gTTGCTGACCCACAACATCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAg                                                                                                                     <  2:321207/148‑1 (MQ=255)
                                 ttGCTGACCCACAACATCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGc                                                                                                                   >  3:113651/1‑149 (MQ=255)
                                      gACCCACAACATCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATaa                                                                                                              >  4:95178/1‑149 (MQ=255)
                                       aCCCACAACATCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATaa                                                                                                              <  8:290380/148‑1 (MQ=255)
                                       aCCCACAACATCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAAt                                                                                                             >  4:204659/1‑149 (MQ=255)
                                       aCCCACAACATCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAAt                                                                                                             >  7:266276/1‑149 (MQ=255)
                                         ccACAACATCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGc                                                                                                           >  6:200288/1‑149 (MQ=255)
                                           acaacaTCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTg                                                                                                         >  3:93106/1‑149 (MQ=255)
                                           acaacaTCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTg                                                                                                         <  4:364449/149‑1 (MQ=255)
                                            caacaTCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTgg                                                                                                        <  2:181288/149‑1 (MQ=255)
                                                   ggCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCa                                                                                                  <  3:164972/148‑1 (MQ=255)
                                                   ggCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAg                                                                                                 <  6:51762/149‑1 (MQ=255)
                                                    gCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGc                                                                                                <  5:200288/149‑1 (MQ=255)
                                                           cGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGccc                                                                                          <  1:253586/148‑1 (MQ=255)
                                                           cGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCt                                                                                         >  6:361461/1‑149 (MQ=255)
                                                            ggAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTg                                                                                                         <  2:282450/132‑1 (MQ=255)
                                                            ggAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTg                                                                                                         >  1:282450/1‑132 (MQ=255)
                                                             gAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTgc                                                                                       <  1:159540/149‑1 (MQ=255)
                                                                  tCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGccag                                                                                  >  1:290550/1‑149 (MQ=255)
                                                                   cGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGccag                                                                                  <  6:227534/148‑1 (MQ=255)
                                                                    gCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGccagca                                                                                >  6:198252/1‑149 (MQ=255)
                                                                     ccGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGccagcac                                                                               <  2:145624/149‑1 (MQ=255)
                                                                        aCTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGccagcaccag                                                                            <  3:262611/149‑1 (MQ=255)
                                                                            ccTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGccagcaccagca                                                                          >  6:339899/1‑147 (MQ=255)
                                                                            ccTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGccagcaccagca                                                                          <  5:339899/147‑1 (MQ=255)
                                                                            ccTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACcagcag                                                                         >  5:165204/1‑148 (MQ=255)
                                                                                  cGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTc                                                                    >  7:213889/1‑147 (MQ=255)
                                                                                  cGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTc                                                                    <  8:213889/147‑1 (MQ=255)
                                                                                   gTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGa                                                                 >  8:279188/1‑149 (MQ=255)
                                                                                      gAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGAt                                                                >  5:349424/1‑147 (MQ=255)
                                                                                           gcagggcagCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCgg                                                         >  1:160187/1‑149 (MQ=255)
                                                                                           gcagggcagCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCgg                                                         <  3:204659/149‑1 (MQ=255)
                                                                                            cagggcagCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCggg                                                        <  4:335126/149‑1 (MQ=255)
                                                                                               ggcagCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCa                                                     <  5:361461/149‑1 (MQ=255)
                                                                                               ggcagCGACCAGGCGCTTCATAAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCa                                                     <  8:178497/149‑1 (MQ=255)
                                                                                                gcagCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATc                                                                                                   >  6:176974/1‑102 (MQ=255)
                                                                                                gcagCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATc                                                                                                   <  5:176974/102‑1 (MQ=255)
                                                                                                   gCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTAc                                                 <  7:279188/149‑1 (MQ=255)
                                                                                                       ccAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGc                                             <  4:93106/149‑1 (MQ=255)
                                                                                                             gcTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGGTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTc                                       <  2:290550/149‑1 (MQ=255)
                                                                                                              cTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCg                                      >  7:1822/1‑149 (MQ=255)
                                                                                                                 cATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCgg                                     <  6:165204/147‑1 (MQ=255)
                                                                                                                   tGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGCAt                                  >  5:74528/1‑148 (MQ=255)
                                                                                                                   tGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGCAt                                  >  5:188508/1‑148 (MQ=255)
                                                                                                                   tGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGCATc                                 <  6:188508/149‑1 (MQ=255)
                                                                                                                        aggcaggTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGCATCAGAcc                            <  6:74528/149‑1 (MQ=255)
                                                                                                                          gcaggTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGCAt                                  <  3:10214/141‑1 (MQ=255)
                                                                                                                          gcaggTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGCAt                                  >  4:10214/1‑141 (MQ=255)
                                                                                                                          gcaggTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGCATCAGAccgc                          >  1:164923/1‑149 (MQ=255)
                                                                                                                          gcaggTCCAGCGCTGGCGCCGTTCGATCGGAGCCTTGCAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCAAGTACCCGCAACTTCGGCATCAGAccgc                          >  5:349900/1‑149 (MQ=255)
                                                                                                                               tCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCCATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGACGCGGGCCAGTACCCGCACCTTCGGCATCAGACCGCCGCCt                     >  4:12114/1‑149 (MQ=255)
                                                                                                                                 cAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGCATCAGACCGCCGCCtt                    >  7:221786/1‑148 (MQ=255)
                                                                                                                                  cttcaaggCGCCGATCGAACGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCAGCAACTTCGGCATCAGACcgccgc                       <  8:325588/138‑1 (MQ=255)
                                                                                                                                  cttcaaggCGCCGATCGAACGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCAGCAACTTCGGCATCAGACcgccgc                       >  7:325588/7‑144 (MQ=255)
                                                                                                                                  aGCGCTGGCGCCGTTCGATCGGCTCCTTTAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCTGGATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGGAACTTCGGCATCAGACCGCCGCCTTGa                  <  6:349900/149‑1 (MQ=255)
                                                                                                                                   gcgcTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGCATCAGACCGCCGCCTTg                   >  5:15538/1‑147 (MQ=255)
                                                                                                                                     caaggCGCCGATCGAACGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGCATCAGACCGCCGCCTTGAAGt               >  7:211247/4‑149 (MQ=255)
                                                                                                                                      aaggCGCCGATCGAACGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGCATCAGACCGCCGCCTTGAAGTg              >  6:105715/3‑149 (MQ=255)
                                                                                                                                         gcgcCGTTCTATCGTCGCCTTGAAGCCCAAGCGTTAGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGCATCAGACCGCCGCCTTGAAGTggc            <  8:266276/148‑1 (MQ=255)
                                                                                                                                         gcgcCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGCATCAGACCGCCGCCTTGAAGTggcg           <  3:208646/149‑1 (MQ=255)
                                                                                                                                          cgcCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGCATCAGACCGCCGCCTTGAAGTggcgg          <  1:227894/149‑1 (MQ=255)
                                                                                                                                          cgcCGATCGAACGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGCATCAGACCGCCGCCTTGAAGTggcgg          <  5:105715/149‑1 (MQ=255)
                                                                                                                                             cGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGCATCAGACCGCCGCCTTGAAGTggcgggcg       >  5:245132/1‑149 (MQ=255)
                                                                                                                                               atCGATCGGCGCCTTGAAGCAAAAGCGTTAGCCCAAAGCGATAATGCTGGCCATCAGCATGCCCTGCGCCATAACCCGCATCACTCGGATGAGGCGCGCCAGTACCCGAAACTACAGCATCAGACCGCCCCAATGAAGTGGAAGGCGt      >  8:156481/2‑148 (MQ=255)
                                                                                                                                                  gATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGCATCAGACCGCCGCCTTGAAGTGGCGGGCGTAGCg  >  2:292408/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CCAAGCAGGGGGATGTCGGTCCAGGAAACATGGTTGCTGACCCACAACATCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGCATCAGACCGCCGCCTTGAAGTGGCGGGCGTAGCG  >  NC_002947/1065582‑1065876

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: ‑