Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 402,065 A→C 5.7% I154S (ATC→AGC)  PP_0335 ← glutathione S‑transferase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947402,0650AC 5.7% 276.4 / 5.9 87I154S (ATC→AGC) PP_0335glutathione S‑transferase family protein
Reads supporting (aligned to +/- strand):  ref base A (39/43);  new base C (3/2);  total (42/45)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

GACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGATCTTGTCACCCTGCGCTTGCAGCCAGCCATCCCAGCGCTTGTCCTCTGGCCGCAGGAAGCTCTCATACCGGCTCGACACTGCCGCATCCATGATCGCATC  >  NC_002947/401923‑402213
                                                                                                                                              |                                                                                                                                                    
gACCATCGATGGCCGCTTGACCGCCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCt                                                                                                                                                >  6:329471/1‑149 (MQ=255)
  ccATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTcc                                                                                                                                              <  1:17683/149‑1 (MQ=255)
   cATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCg                                                                                                                                             >  4:268726/1‑149 (MQ=255)
   cATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCACTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCg                                                                                                                                             <  5:302587/149‑1 (MQ=255)
   cATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCACTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCg                                                                                                                                             >  6:302587/1‑149 (MQ=255)
       gATGGCCGCTTGACCGCCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCaa                                                                                                                                         <  5:329471/149‑1 (MQ=255)
             cGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGtgctc                                                                                                                                   <  4:190928/149‑1 (MQ=255)
             cGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCAAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGtgctc                                                                                                                                   >  2:214277/1‑149 (MQ=255)
               cTTGACCACCTCGGTATTCGACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGtgctcc                                                                                                                                  <  6:252039/148‑1 (MQ=255)
                  gACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGtgctcctg                                                                                                                                <  7:205586/147‑1 (MQ=255)
                  gACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGtgctcc                                                                                                                                  <  6:252036/145‑1 (MQ=255)
                  gACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGtgctcc                                                                                                                                  >  5:252036/1‑145 (MQ=255)
                  gACAACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCAAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGtgctcc                                                                                                                                  >  5:252039/1‑145 (MQ=255)
                    ccaccTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGtgctcctgctc                                                                                                                             >  4:151924/1‑148 (MQ=255)
                     caccTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGAGGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGtgctcctgctc                                                                                                                             >  8:113128/1‑147 (MQ=255)
                     caccTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGtgctcctgctc                                                                                                                             >  4:180862/1‑147 (MQ=255)
                      accTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAg                                                                                                                          <  4:158676/149‑1 (MQ=255)
                       ccTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAgg                                                                                                                         <  3:268726/149‑1 (MQ=255)
                            gCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTc                                                                                                                    <  3:180862/149‑1 (MQ=255)
                              aTACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTc                                                                                                                    >  8:256213/1‑147 (MQ=255)
                                   aCGCCGCCAACCCCGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTa                                                                                                             >  5:237383/1‑149 (MQ=255)
                                   aCGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGCTGATCTCCGGc                                                                                                                                           >  4:66257/1‑119 (MQ=255)
                                   aCGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGCTGATCTCCGGc                                                                                                                                           <  3:66257/119‑1 (MQ=255)
                                   aCGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAAGACAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTa                                                                                                             >  1:273127/1‑149 (MQ=255)
                                        gccAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTAcggcg                                                                                                        >  5:45608/1‑149 (MQ=255)
                                        gccAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTAcggcg                                                                                                        >  5:40227/1‑149 (MQ=255)
                                         ccAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTAcggcgg                                                                                                       <  8:286335/149‑1 (MQ=255)
                                            aCCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGATc                                                                                                    <  4:273984/149‑1 (MQ=255)
                                              ccAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCGAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTCTCCAGGCTACGGCGGATCtt                                                                                                  >  4:198089/1‑149 (MQ=255)
                                              ccAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGATCtt                                                                                                  >  7:273397/1‑149 (MQ=255)
                                               cAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGATCTTg                                                                                                 <  1:127612/149‑1 (MQ=255)
                                                      gcgcTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAg                                                                                                                 <  7:302420/126‑1 (MQ=255)
                                                      gcgcTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAg                                                                                                                 >  8:302420/1‑126 (MQ=255)
                                                        gcTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTACTCCAGGTTGTCCAGGCTACGGCGGATCTTGTCACCCTgc                                                                                        >  8:109086/1‑149 (MQ=255)
                                                            aCGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGATCTTGTCACCCTGCGCtt                                                                                    <  8:273397/149‑1 (MQ=255)
                                                              gccagccaAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGa                                                                                                      <  8:232342/129‑1 (MQ=255)
                                                              gccagccaAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGa                                                                                                      >  7:232342/1‑129 (MQ=255)
                                                              gccagccaAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGATCTTg                                                                                                 >  1:139367/1‑134 (MQ=255)
                                                              gccagccaAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGATCTTg                                                                                                 <  2:139367/134‑1 (MQ=255)
                                                               ccagccaAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGATCTTGTCAccc                                                                                           >  8:120393/1‑139 (MQ=255)
                                                               ccagccaAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGATCTTGTCAccc                                                                                           <  7:120393/139‑1 (MQ=255)
                                                               ccagccaAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGATCTTGTCACCCTGCGCTTGca                                                                                 <  6:45608/149‑1 (MQ=255)
                                                               ccagccaAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGATCTTGTCACCCTGCGCTTGca                                                                                 <  6:40227/149‑1 (MQ=255)
                                                                     aaaCCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAgg                                                                                                                         <  6:352237/103‑1 (MQ=255)
                                                                     aaaCCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAgg                                                                                                                         >  5:352237/1‑103 (MQ=255)
                                                                     aaaCCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGATCTTGTCACCCTGCGCTTGcagccagc                                                                           >  3:347769/1‑149 (MQ=255)
                                                                           gggCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGATCTTGTCACCCTGCGCTTGCAGCCAGCCATccc                                                                     <  4:347769/149‑1 (MQ=255)
                                                                            ggCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGATCTTGTCACCCTGCGCTTGCAGCCAGCCATCCCa                                                                    <  7:256213/149‑1 (MQ=255)
                                                                            ggCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGATCTTGTCACCCTGCGCTTGCAGCCAGCCATCCCa                                                                    <  3:198089/149‑1 (MQ=255)
                                                                                 aCGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGATCTTGt                                                                                                <  2:279676/116‑1 (MQ=255)
                                                                                 aCGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGATCTTGt                                                                                                >  1:279676/1‑116 (MQ=255)
                                                                                  cGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTACTCCAGGTTGTCCAGGCTACGGCGGATCTTGTCACCCTGCGCTTGCAGCCAGCCATCCCAGCGCtt                                                              <  7:109086/149‑1 (MQ=255)
                                                                                   gCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTcc                                                                                                                   <  6:278793/95‑1 (MQ=255)
                                                                                   gCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTcc                                                                                                                   >  5:278793/1‑95 (MQ=255)
                                                                                     aaaTCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGATCTTGTCACCCTGCGCTTGCAGCCAGCCATCCCAGCGCTTGTc                                                           >  7:236444/1‑149 (MQ=255)
                                                                                     aaaTCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGATCTTGTCACCCTGCGCTTGCAGCCAGCCATCCCAGCGCTTGTc                                                           >  1:18260/1‑149 (MQ=255)
                                                                                     aaaTCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGACGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGATCTTGTCACCCTGCGCTTGCAGCCAGCCATCCCAGCGCTTGTc                                                           <  6:315809/149‑1 (MQ=255)
                                                                                          cAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGATCTTGTCACCCTGCGCTTGCAGCCAGCCATCCCAGCGCTTGTCctct                                                       <  3:49849/148‑1 (MQ=255)
                                                                                          cAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGATCTTGTCACCCTGCGCTTGCAGCCAGCCATCCCAGCGCTTGTCCTCTg                                                      <  6:237383/149‑1 (MQ=255)
                                                                                                  ccAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGATCTTGTCACCCTGCGCTTGCAGCCAGCCATCCCAGCGCTTGTCCTCTGGCCGCAg                                               >  2:181036/1‑148 (MQ=255)
                                                                                                    aGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGATCTTGTCACCCTGCGCTTGCAGCCAGCCATCCCAGCGCTTGTCCTCTGGCCGCAgg                                              >  6:99794/1‑147 (MQ=255)
                                                                                                     gcgcACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGATCTTGTCACCCTGCGCTTGCAGCCAGCCATCCCAGCGCTTGTCCTCTGGCCGCAGGAAg                                           <  2:273127/149‑1 (MQ=255)
                                                                                                     gcgcACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGATCTTGTCACCCTGCGCTTGCAGCCAGCCATCCCAGCGCTTGTCCTCTGGCCGCAGGAAg                                           <  4:257269/149‑1 (MQ=255)
                                                                                                         acaGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGATCTTGTCACCCTGCGCTTGCAGCCAGCCATCCCAGCGCTTGTCCTCTGGCCGCAGGAAGctct                                       <  8:78354/149‑1 (MQ=255)
                                                                                                           aGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGATCTTGTCACCCTGCGCTTGCAGCCAGCCATCCCAGCGCTTGTc                                                           <  6:244578/127‑1 (MQ=255)
                                                                                                           aGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGATCTTGTCACCCTGCGCTTGCAGCCAGCCATCCCAGCGCTTGTc                                                           >  5:244578/1‑127 (MQ=255)
                                                                                                             gCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGATCTTGTCACCCTGCGCTTGCAGCCAGCCATCCCAGCGCTTGTCCTCTGGCCGCAGGAAGCTCTCATa                                   >  2:241015/1‑149 (MQ=255)
                                                                                                                aCGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGATCTTGTCACCCTGCGCTTGCAGCCAGCCATCCCAGCGCTTGTCCTCTGGCCGCAGGAAGCTCTCATACCg                                >  5:326144/1‑149 (MQ=255)
                                                                                                                  gCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGATCTTGTCACCCTGCGCTTGCAGCCAGCCATCCCAGCGCTTGTCCTCTGGCCGCAGGAAg                                           <  2:190571/136‑1 (MQ=255)
                                                                                                                  gCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGATCTTGTCACCCTGCGCTTGCAGCCAGCCATCCCAGCGCTTGTCCTCTGGCCGCAGGAAg                                           >  1:190571/1‑136 (MQ=255)
                                                                                                                  gCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGATCTTGTCACCCTGCGCTTGCAGCCAGCCATCCCAGCGCTTGTCCTCTGGCCGCAGGAAGCTCTCATACCGGc                              <  7:113128/149‑1 (MQ=255)
                                                                                                                   ccAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGATCTTGTCACCCTGCGCTTGCAGCCAGCCATCCCAGCGCTTGTCCTCTGGCCGCAGGAAGCTCTCATACCGGCt                             <  1:241015/149‑1 (MQ=255)
                                                                                                                     aaTCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGATCTTGTCACCCTGCGCTTGCAGCCAGCCATCCCAGCGCTTGTCCTCTGGc                                                    <  4:1769/124‑1 (MQ=255)
                                                                                                                     aaTCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGATCTTGTCACCCTGCGCTTGCAGCCAGCCATCCCAGCGCTTGTCCTCTGGc                                                    >  3:1769/1‑124 (MQ=255)
                                                                                                                     aaTCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGATCTTGTCACCCTGCGCTTGCAGCCAGCCATCCCAGCGCTTGTCCTCTGGCCGCAGGAAGCTCTCATACCGGCTc                            >  3:310648/1‑148 (MQ=255)
                                                                                                                               aGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGATCTTGTCACCCTGCGCTTGCAGCCAGCCATCCCAGCGCTTGTCCTCTGGCCGCAGGAAGCTCTCATACCGGCTCGACACTGCCGc                 <  1:181036/149‑1 (MQ=255)
                                                                                                                                ggTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGATCTTGTCACCCTGCGCTTGCAGCCAGCCATCCCAGCGCTTGTCCTCTGGCCGCAGGAAGCTCTCATACCGGCTCGACACTGCCGCa                >  6:299719/1‑149 (MQ=255)
                                                                                                                                  tCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCTGGCTACGGCGGATGTTGTCACCCTGCGCTTGCAGCCAGCCATGCCAGCGCTTGTCCTCTGGCCGGAGGAAGCTCTCATACCGGCTCGAGACTGCCGCATc              <  8:43874/149‑1 (MQ=255)
                                                                                                                                  tCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGATCTTGTCACCCTGCGCTTGCAGCCAGCCATCCCAGCGCTTGTCCTCTGGCCGCAGGAAGCTCTCATACCGGCTCGACACTGCCGCATc              <  8:38253/149‑1 (MQ=255)
                                                                                                                                  tCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGATCTTGTCACCCTGCGCTTGCAGCCAGCCATCCCAGCGCTTGTCCTCTGGCCGCAGGAAGCTCTCATACCGGCTCGACACTGCCGCATc              <  4:310648/149‑1 (MQ=255)
                                                                                                                                   cGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGATCTTGTCACCCTGCGCTTGCAGCCAGCCATCCCAGCGCTTGTCCTCTGGCCGCAGGAAGCTCTCATACCGGCTCGACACTGCCGCATc              >  3:163203/1‑148 (MQ=255)
                                                                                                                                    gAAGCCGGAGATCAGCTCCGGCAAGTGCTCGTGCTCCAGGTTGTCCAGGCTACGGCGGATCTTGTCACCCTGCGCTTGCAGCCAGCCATCCCAGCGCTTGTCCTCTGGTCGCAGGAAGCTCTCATACCGGCTCGACACTGCCGCATCCa            >  1:29168/1‑149 (MQ=255)
                                                                                                                                     ttgCCGGAGCTGATCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGATCTTGTCACCCTGCGCTTGCAGCCAGCCATCCCAGCGCTTGTCCTCTGGCCGCAGGAAGCTCTCATACCGGCTCGACACTGCCGCATCCAt           >  5:1505/3‑149 (MQ=255)
                                                                                                                                       gCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGATCTTGTCACCCTGCGCTTGCAGCCAGCCATCCCAGCGCTTGTCCTCTGGCCGCAGGAAGCTCTCATACCGGCTCGACACTGCCGCATCCATGa         >  2:184467/1‑149 (MQ=255)
                                                                                                                                         cGGAGCTGATCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGATCTTGTCACCCTGCGCTTGCAGCCAGCCATCCCAGCGCTTGTCCTCTGGCCGCAGGAAGCTCTCATACCGGCTc                            >  3:42283/1‑128 (MQ=255)
                                                                                                                                         cGGAGCTGATCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGATCTTGTCACCCTGCGCTTGCAGCCAGCCATCCCAGCGCTTGTCCTCTGGCCGCAGGAAGCTCTCATACCGGCTc                            <  4:42283/128‑1 (MQ=255)
                                                                                                                                         cGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGATCTTGTCACCCTGCGCTTGCAGCCAGCCATCCCAGCGCTTGTCCTCTGGCCGCAGGAAGCTCTCATACCGGCTCGACACTGCCGCATCCATGATc       <  5:299719/149‑1 (MQ=255)
                                                                                                                                           gagaTCAGCTCCGGCAAGTGCTCCTGCTCCATGTTGTCCAGGCTACGGCGGATCTTGTCACCCTGCGCTTGCAGCCAGCCATCCCAGCGCTTGTCCTCTGGCCGCAGGAAGCTCTCATACCGGCTCGACACTGCCGCATCCATGATCg      >  4:332332/1‑148 (MQ=255)
                                                                                                                                             gaTCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGATCTTGTCACCCTGCGCTTGCAGCCAGCCATCCCAGCGCTTGTCCTCTGGCCGCAGGAAGCTCTCATACCGGCTCGACACTGCCGCATCCATGATCGCAt   <  1:214277/149‑1 (MQ=255)
                                                                                                                                              aTCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGATCTTGTCACCCTGCGCTTGCAGCCAGCCATCCCAGCGCTTGTCCTCTGGCCGCAGGAAGCTCTCATACCGGCTCGACACTGCCGCATCCATGATCGCATc  <  2:18260/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
GACCATCGATGGCCGCTTGACCACCTCGGCATACCACGCCGCCAACCCAGGCTGGCGCTCACGCCAGCCAAACCCGGGCTGACGCAAATCCAGATAGCCCAGCGCACAGGCCACGCCAATCGCCGCCAGGTCGAAGCCGGAGATCAGCTCCGGCAAGTGCTCCTGCTCCAGGTTGTCCAGGCTACGGCGGATCTTGTCACCCTGCGCTTGCAGCCAGCCATCCCAGCGCTTGTCCTCTGGCCGCAGGAAGCTCTCATACCGGCTCGACACTGCCGCATCCATGATCGCATC  >  NC_002947/401923‑402213

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