Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,142,771 T→C 7.5% F302L (TTC→CTC)  PP_3646 → aldehyde dehydrogenase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,142,7710TC 7.5% 181.0 / 5.3 68F302L (TTC→CTC) PP_3646aldehyde dehydrogenase family protein
Reads supporting (aligned to +/- strand):  ref base T (33/29);  new base C (3/2);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.57e-01

ACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGCCCGATGGCCACCGCGCAACAGCTGGCAGTCGTCGAAGGCCTGGTGGCGGCGGCCAAGGCCGAAGGCGCCAAAC  >  NC_002947/4142623‑4142918
                                                                                                                                                    |                                                                                                                                                     
aCTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGt                                                                                                                                                       <  4:73315/149‑1 (MQ=255)
aCTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGt                                                                                                                                                       <  4:68168/149‑1 (MQ=255)
aCTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGt                                                                                                                                                       <  1:383177/149‑1 (MQ=255)
    ggcggcAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTtcgt                                                                                                                                                   >  6:264149/1‑149 (MQ=255)
      cggcAAGTCGCCGAACATCTTCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTtcgtc                                                                                                                                                  >  5:96555/1‑148 (MQ=255)
       ggcAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTtcgtc                                                                                                                                                  >  1:167275/1‑147 (MQ=255)
        gcAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAg                                                                                                                                               <  6:303646/149‑1 (MQ=255)
        gcAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAACTCGTCGaa                                                                                                                                               <  2:160373/149‑2 (MQ=255)
         cAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAgc                                                                                                                                              <  1:248024/149‑1 (MQ=255)
            gTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGcc                                                                                                                                           <  2:127039/149‑1 (MQ=255)
            gTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGcc                                                                                                                                           >  7:354385/1‑149 (MQ=255)
            gTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGcc                                                                                                                                           <  2:167275/149‑1 (MQ=255)
               gCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGa                                                                                                                                        >  7:329699/1‑149 (MQ=255)
                 cGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGAtt                                                                                                                                     <  2:197209/149‑1 (MQ=255)
                     catcatCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGccc                                                                                                                                 <  2:229801/149‑1 (MQ=255)
                       tcatcTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCcgc                                                                                                                               >  3:172227/1‑149 (MQ=255)
                         atcTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCgcgc                                                                                                                             >  7:32635/1‑149 (MQ=255)
                              cgccgacgccgaCCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCaa                                                                                                                          <  7:233611/147‑1 (MQ=255)
                               gccgacgccgacCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAgc                                                                                                                        <  5:207427/148‑1 (MQ=255)
                                     gccgacCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGcgcatccg                                                                                                                 >  3:160834/1‑149 (MQ=255)
                                      ccgacCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGcgcatccgc                                                                                                                <  4:86797/149‑1 (MQ=255)
                                         acCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGcgcatccgcatc                                                                                                             <  2:165869/149‑1 (MQ=255)
                                         acCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGcgcatccgcatc                                                                                                             >  4:333403/1‑149 (MQ=255)
                                                    gcCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCa                                                                                                  <  3:333403/149‑1 (MQ=255)
                                                     ccATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAg                                                                                                 <  5:218153/149‑1 (MQ=255)
                                                      cATCAATGGCGCCGTGGCCGGCATATATGCCGCGTCCGGCCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGAGCAGGATGAAATCTTCGACGAGTTCGACGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAg                                                                                                 >  2:219569/1‑148 (MQ=255)
                                                       aTCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGTCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGga                                                                                               <  1:276161/149‑1 (MQ=255)
                                                          aaTGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGgacga                                                                                           >  5:178189/1‑149 (MQ=255)
                                                          aaTGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGgacga                                                                                           <  6:122589/149‑1 (MQ=255)
                                                                   gTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCg                                                                                   >  3:97078/1‑148 (MQ=255)
                                                                        cGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATgg                                                                             >  6:214806/1‑149 (MQ=255)
                                                                         ggCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATggg                                                                            <  8:32635/149‑1 (MQ=255)
                                                                          gCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGc                                                                           <  8:354385/149‑1 (MQ=255)
                                                                               tATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGCCCGAt                                                                      >  4:262043/1‑149 (MQ=255)
                                                                                aTGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGCCCGATg                                                                     >  5:83165/1‑149 (MQ=255)
                                                                                 tGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGCCCGATg                                                                     <  5:31910/148‑1 (MQ=255)
                                                                                   ccGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGCCCGATGGc                                                                   >  7:130869/1‑148 (MQ=255)
                                                                                    cgcgTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGCCCGATGGCCa                                                                 >  8:109396/1‑149 (MQ=255)
                                                                                       gTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGCCCGATGGCCACcg                                                              >  1:1228/1‑149 (MQ=255)
                                                                                               aaaGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCATCGAAATGGGCCCGATGGCCACCGCGCAACAg                                                      >  6:161220/1‑149 (MQ=255)
                                                                                               aaaGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGCCCGATGGCCACCGCGCAACAg                                                      >  4:173605/1‑149 (MQ=255)
                                                                                                     gCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGCGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGCCCGATGGCCACCGCGCAACAGCTgg                                                  >  3:39135/1‑147 (MQ=255)
                                                                                                      cGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGAAGACGCCAGCGAAATGGGCCCGATGGCCACCGCGCAACAGCTGGCa                                                <  1:276574/148‑1 (MQ=255)
                                                                                                        tGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGCCCGATGGCCACCGCGCAACAGCTGGCAgtc                                             >  2:198884/1‑149 (MQ=255)
                                                                                                           ccGGCTCGCGCCTGGTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGTACGCCAGCGAAATGGGCCCGATGGCCACCGCGCAACAGCTGGCAgtcgt                                           <  1:91829/149‑1 (MQ=255)
                                                                                                            cGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCACAGGACGACGCCAGCGAAATGGGCCCGATGGCCACCGCGCAACAGCTGGCAgtcgtc                                          >  8:303185/1‑148 (MQ=255)
                                                                                                              gCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGCCCGATGGCCACCGCGCAACAGCTGGCAGTCGTCGaa                                       <  3:262043/149‑1 (MQ=255)
                                                                                                                  gcgcCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGCCCGATGGCCACCGCGCAACAGCTGGCAGTCGTCGAAgg                                     <  6:178189/147‑1 (MQ=255)
                                                                                                                  gcgcCTGCTGGTGCAGGATGAAATCTTCGACGAGTCCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGCCCGATGGCCACCGCGCAACAGCTGGCAGTCGTCGAAGGc                                    <  8:49800/148‑1 (MQ=255)
                                                                                                                     ccTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCCCAGGACGACGCCAGCGAAATGGGCCCGATGGCCACCGCGCAACAGCTGGCAGTCGTCGAAGGCCtgg                                >  4:113058/1‑149 (MQ=255)
                                                                                                                        gctgGTGCAGGATGAAATCTTCGACGAGTTCGTCAAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGCCCGATGGCCACCGCGCAACAGCTGGCAGTCGTCGAAGGCCtggt                               <  1:146087/147‑1 (MQ=255)
                                                                                                                        gctgGTGCAGGATGAAATCTTCGACGAGTTCGTCAAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGCCCGATGGCCACCGCGCAACAGCTGGCAGTCGTCGAAGGCCtggt                               >  2:146087/1‑147 (MQ=255)
                                                                                                                                 ggATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAATGGACCCGATGGCCACCGCGCAACAGCTGGCAGTCGTCGAAGGCCTGGTGGCGGCGGcc                     >  2:353547/1‑147 (MQ=255)
                                                                                                                                     gAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGCCCGATGGCCACCGCGCAACAGCTGGCAGTCGTCGAAGGCCTGGTGGCGGCGGCCAAGGc                >  6:54506/1‑149 (MQ=255)
                                                                                                                                     gAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGCCCGATGGCCACCGCGCAACAGCTGGCAGTCGTCGAAGGCCTGGTGGCGGCGGCCAAGGc                >  5:22056/1‑149 (MQ=255)
                                                                                                                                       aaTCTTCTACGAGTTCGTCGAGCGCCTGATTGCCCTCGCCAAGCGCATCCGCATGTGCAACCGGCAGGACGACGCCAGCGAAATGGGCCCGATGGCCACCGCGCAAGAGCTGGCAGTCGTTGAAGGCCTGGTGGCGGCGGCCAAGGCCg              <  1:219569/149‑1 (MQ=255)
                                                                                                                                        aTCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGCCCGATGGCCACCGCGCAACAGCTGGCAGTCGTCGAAGGCCTGGTGGCGGCGGCCAAGGCCGa             <  6:96555/149‑1 (MQ=255)
                                                                                                                                          gctCGACGAACTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGCCCGATGGCCACCGCGCAACAGCTGGCAGTCGTCGAAGGCCTGGTggcggcg                        >  5:137204/3‑136 (MQ=255)
                                                                                                                                          gctCGACGAACTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGCCCGATGGCCACCGCGCAACAGCTGGCAGTCGTCGAAGGCCTGGTggcggcg                        <  6:137204/134‑1 (MQ=255)
                                                                                                                                           ttCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGCCCGATGGCCACCGCGCAACAGCTGGCAGTTGTCGAAGGCCTGGTGGCGGCGGCCAAGGCCGAAgg          <  8:130869/149‑1 (MQ=255)
                                                                                                                                           ttCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGCCCGATGGCCACCGCGCAACAGCTGGCAGTCGTCGAAGGCCTGGTGGCGGCGGCCAAGGCCGAAgg          >  6:220190/1‑149 (MQ=255)
                                                                                                                                           ctCGACGAACTCGTCGAGCGCCTGATGTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGCCCGATGGCCACCGCGCAACAGCTGGCAGTCGTCGAAGGCCTGGTGGCGGCGGCCAAGGCCGAAg           >  1:5158/2‑149 (MQ=255)
                                                                                                                                           ctCGACGAACTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGCCCGATGGCCACCGCGCAACAGCTGGCAGTCGTCGAAGGCCTGGTGGCGGCGGCCAAGGCCGAAgg          >  7:92619/2‑149 (MQ=255)
                                                                                                                                             cgacgaGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGCCCGATGGCCACCGCGCAACAGCTGGCAGTCGTCGAAGGCCTGGTGGCGGCGGCCAAGGCCGAAGgcg        >  2:122337/1‑149 (MQ=255)
                                                                                                                                             cgacgaGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGCCCGATGGCCACCGCGCAACAGCTGGCAGTCGTCGAAGGCCTGGTGGCGGCGGCCAAGGCCGAAGgc         <  6:83165/148‑1 (MQ=255)
                                                                                                                                             cgacgaGATCGTCGAGCGGCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGCCCGATGGCCACCGCGCAACAGCTGGCAGTCGTCGAAGGCCTGGTGGCGGCGGCCAAGGCCGAAGgc         <  8:329699/148‑1 (MQ=255)
                                                                                                                                              gacgaGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGCCCGATGGCCACCGCGCAACAGCTGGCAGTCGTCGAAGGCCTGGTGGCGGCGGCCAAGGCCGAAGgcgc       >  4:67155/1‑149 (MQ=255)
                                                                                                                                               acgaGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGCCCGATGGCCACCGCGCAACAGCTGGCAGTCGTCGAAGGCCTGGTGGCGGCGGCCAAGGCCGAAGgcgc       >  7:293283/1‑148 (MQ=255)
                                                                                                                                                   gTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGCCCGATGGCCACCGCGCAACAGCTGGCAGTCGTCGAAGGCCTGGTGGCGGCGGCCAAGGCCGAAGGCGCCAAAc  <  4:172227/149‑1 (MQ=255)
                                                                                                                                                   gTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGCCCGATGGCCACCGCGCAACAGCTGGCAGTCGTCGAAGGCCTGGTGGCGGCGGCCAAGGCCGAAGGCGCCAAAc  <  1:122337/149‑1 (MQ=255)
                                                                                                                                                   gTTCGTCGAGCGCCTGATTGCCCGAGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGCCCGATGGCCACCGCGCAACAGCTGGCAGTCGTCGAAGGCCTGGTGGAGGCGGCCAAGGCCGAAGGCGCCaacc  >  3:340554/1‑147 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ACTGGGCGGCAAGTCGCCGAACATCATCTTCGCCGACGCCGACCTGGACAGCGCCATCAATGGCGCCGTGGCCGGCATCTATGCCGCGTCCGGGCAAAGCTGCGTGGCCGGCTCGCGCCTGCTGGTGCAGGATGAAATCTTCGACGAGTTCGTCGAGCGCCTGATTGCCCGCGCCAAGCGCATCCGCATCGGCAACCCGCAGGACGACGCCAGCGAAATGGGCCCGATGGCCACCGCGCAACAGCTGGCAGTCGTCGAAGGCCTGGTGGCGGCGGCCAAGGCCGAAGGCGCCAAAC  >  NC_002947/4142623‑4142918

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: