Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,726,416 C→A 7.8% Q357K (CAA→AAA)  PP_3291 → metallo‑beta‑lactamase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,726,4160CA 7.8% 260.0 / 8.3 85Q357K (CAA→AAA) PP_3291metallo‑beta‑lactamase family protein
Reads supporting (aligned to +/- strand):  ref base C (41/37);  new base A (3/4);  total (44/41)
Fisher's exact test for biased strand distribution p-value = 7.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCGATGAGTCCAGGCACCTGGACCAGCTGCGCGAGCCCGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGATGGAGAATCTTTCTGCGCACGCCGATGCTGATGAAATCATGGAATGGTT  >  NC_002947/3726272‑3726561
                                                                                                                                                |                                                                                                                                                 
gTCGATGAGTCCAGGCACCTGGACCAGCTGCGCGAGCCCGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGc                                                                                                                                               >  4:144075/1‑149 (MQ=255)
 tCGATGAGTCCAGGCACCTGGACCAGCTGCGCGAGCCCGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGcc                                                                                                                                              >  3:43482/1‑149 (MQ=255)
   gATGAGTCCAGGCACCTGGACCAGCTGCGCGAGCCCGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCgg                                                                                                                                            <  4:43482/149‑1 (MQ=255)
   gATGAGTCCAGGCACCTGGACCAGCTGCGCGAGCCCGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTGAAAGCCg                                                                                                                                             <  8:35448/148‑1 (MQ=255)
    aTGAGTCCAGGCACCTGGACCAGCTGCGCGAGCCCGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCggg                                                                                                                                           >  7:66087/1‑149 (MQ=255)
    aTGAGTCCAGGCACCTGGACCAGCTGCGCGAGCCCGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCggg                                                                                                                                           <  2:41408/149‑1 (MQ=255)
     tGAGTCCAGGCACCTGGACCAGCTGCGCGAGCCCGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCgggg                                                                                                                                          <  7:257498/149‑1 (MQ=255)
     ctggTCCAGGTGCCTGGACCAGCTGCGCGAGCCCGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCgggg                                                                                                                                          >  7:33692/4‑149 (MQ=255)
      tggTCCAGGTGCCTGGACCAGCTGCGCGAGCCCGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCggggg                                                                                                                                         >  1:39416/3‑149 (MQ=255)
       aGTCCAGGCACCTGGACCAGCTGCGCGAGCCCGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCggggg                                                                                                                                         <  5:245766/148‑1 (MQ=255)
         tCCAGGCACCTGGACCAGCTGCGCGAGCCCGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGAc                                                                                                                                      >  3:156911/1‑149 (MQ=255)
         tCCAGGCACCTGGACCAGCTGCGCGAGCCCGCAGTGATCATCACCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGAc                                                                                                                                      >  3:372086/1‑149 (MQ=255)
           cAGGCACCTGGACCAGCTGCGCGAGCCCGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGAcgc                                                                                                                                    >  7:166826/1‑149 (MQ=255)
           cAGGCACCTGGACCAGCTGCGCGAGCCCGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCAGGCTTTCAAGCAGGGGGGCcgc                                                                                                                                    >  5:205960/1‑149 (MQ=255)
            aGGCACCTGGACCATCTGCGCGATCGCGCAGTGATCATCGCCGGCAGCGGGATGGCGACCGGCGGACGGGTGCTTCACGACCTGAATGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGGCGGGGGGAcgcg                                                                                                                                   <  8:107394/149‑1 (MQ=255)
            aGGCACCTGGACCAGCTGCGCGAGCCCGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGAcgcg                                                                                                                                   >  3:299144/1‑149 (MQ=255)
            aGGCACCTGGACCAGCTGCGCGAGCCCGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCCCCGGCTTGAAAGCCggagaa                                                                                                                                        <  7:233815/144‑5 (MQ=255)
            aGGCACCTGGACCAGCTGCGCGAGCCCGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCCCCGGCTTGAAAGCCggagaa                                                                                                                                        >  8:233815/1‑140 (MQ=255)
             ggCACCTGGACCAGCTGCGCGAGCCCGCATTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCgg                                                                                                                                  <  6:300391/149‑1 (MQ=255)
               cACCTGGACCAGCTGCGCGAGCCCGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCgggg                                                                                                                                <  4:299144/149‑1 (MQ=255)
                  cTGGACCAGCTGCGCGAGCCCGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGc                                                                                                                             <  1:339099/149‑1 (MQ=255)
                      aCCAGCTGCGCGAGCCCGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGTGTGCCGa                                                                                                                          >  8:352532/1‑148 (MQ=255)
                      aCCAGCTGCGCGAGCCCGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGa                                                                                                                          >  3:175629/1‑148 (MQ=255)
                      aCCAGCTGCGCGAGCCCGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGa                                                                                                                          >  3:33949/1‑148 (MQ=255)
                              cgcgAGCCCGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGACCCAAACCCGCGCAACAGCCTTCTGTTCTCCGGCTTACAAGCCGGGGGGACGCGGGGTGCCGATATCGTcgc                                                                                                                 >  7:35738/1‑149 (MQ=255)
                               gcgAGCCCGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTcgcg                                                                                                                >  6:83182/1‑149 (MQ=255)
                                    tccGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGgcgc                                                                                                           >  8:342965/2‑149 (MQ=255)
                                    cccGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGgcgc                                                                                                           <  4:213445/149‑1 (MQ=255)
                                     ccGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTGAAAGCCggag                                                                                                                                          >  1:246856/1‑115 (MQ=255)
                                     ccGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTGAAAGCCggag                                                                                                                                          <  2:246856/117‑3 (MQ=255)
                                       gCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCca                                                                                                        <  4:156911/149‑1 (MQ=255)
                                       gCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCca                                                                                                        >  7:337117/1‑149 (MQ=255)
                                              tcatcGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGtgtg                                                                                                 >  8:133602/1‑149 (MQ=255)
                                                atcGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGc                                                                                                <  6:44789/148‑1 (MQ=255)
                                                atcGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCg                                                                                               <  7:352532/149‑1 (MQ=255)
                                                 tcGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCgggg                                                                                                                                          >  6:157579/1‑105 (MQ=255)
                                                 tcGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCgggg                                                                                                                                          <  5:157579/105‑1 (MQ=255)
                                                 tcGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCgg                                                                                              <  7:342965/149‑1 (MQ=255)
                                                 tcGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCgg                                                                                              >  2:216873/1‑149 (MQ=255)
                                                     cgccAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGc                                                                                          <  8:33692/149‑1 (MQ=255)
                                                          gCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGcc                                                                                                          >  6:36511/1‑128 (MQ=255)
                                                          gCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGcc                                                                                                          <  5:36511/128‑1 (MQ=255)
                                                               aTGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGa                                                                                <  3:144075/149‑1 (MQ=255)
                                                                  gCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGAtgt                                                                             <  4:33949/149‑1 (MQ=255)
                                                                  gCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGAtgt                                                                             <  3:307613/149‑1 (MQ=255)
                                                                  gCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGAtgt                                                                             <  4:175629/149‑1 (MQ=255)
                                                                    gACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGc                                                                           <  5:83182/149‑1 (MQ=255)
                                                                      gcggcggACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGc                                                                           <  5:86918/146‑1 (MQ=255)
                                                                             aCGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGt                                                                   >  3:229066/1‑148 (MQ=255)
                                                                                    cTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGt                                                           >  3:156048/1‑149 (MQ=255)
                                                                                       caccacCTGAAGGCACTGGCCCCAAACCCGCGCAACGGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGCCGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCa                                                        >  7:2641/1‑149 (MQ=255)
                                                                                       caccacCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCa                                                        >  3:101964/1‑149 (MQ=255)
                                                                                        accaccTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCAt                                                       >  5:393344/1‑149 (MQ=255)
                                                                                         ccaccTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATg                                                      <  4:312362/149‑1 (MQ=255)
                                                                                          caccTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGc                                                     <  4:101964/149‑1 (MQ=255)
                                                                                             cTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGAt                                                  >  2:334953/1‑149 (MQ=255)
                                                                                                aaGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGATgg                                                <  1:216873/148‑1 (MQ=255)
                                                                                                 aGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGATGga                                               <  4:229066/148‑1 (MQ=255)
                                                                                                     aCTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCca                                                                                                        >  1:385971/1‑87 (MQ=255)
                                                                                                     aCTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCca                                                                                                        <  2:385971/87‑1 (MQ=255)
                                                                                                      cTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGATGGAGAATCt                                         >  5:283370/1‑149 (MQ=255)
                                                                                                      cTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGATGGAGAATCt                                         >  8:249388/1‑149 (MQ=255)
                                                                                                           cccAAACCCGCGCAACCGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGATGGAGAATCTTTCTg                                    <  7:209187/149‑1 (MQ=255)
                                                                                                           cccAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGATGGAGAATCTTTCTg                                    >  8:191014/1‑149 (MQ=255)
                                                                                                                        aaCAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTg                                                                  <  6:317346/106‑1 (MQ=255)
                                                                                                                        aaCAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTg                                                                  >  5:317346/1‑106 (MQ=255)
                                                                                                                          cAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGATGGAGAATCTTTCt                                     >  2:194196/1‑133 (MQ=255)
                                                                                                                          cAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGATGGAGAATCTTTCt                                     >  4:313214/1‑133 (MQ=255)
                                                                                                                                  tGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGAGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGATGGAGAATCTTTCTGCGCACGCCGATGCTGATGAAAt              >  3:245449/1‑148 (MQ=255)
                                                                                                                                    tggtcCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGTCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGATGGAGAATCTTTCTGCGCACGCCGATGCTGATGAAAt              <  8:337117/143‑1 (MQ=255)
                                                                                                                                    ccccccGGCTTGAAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGATGGAGAATCTTTCTGCGCACGCCGATGCTGATGAAATCAt           >  1:272918/5‑149 (MQ=255)
                                                                                                                                     tctcCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGATGGAGAATCTTTCTGCGCACGCCGa                          <  6:296704/133‑1 (MQ=255)
                                                                                                                                     tctcCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGATGGAGAATCTTTCTGCGCACGCCGa                          >  5:296704/1‑133 (MQ=255)
                                                                                                                                      ctcCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGATGGAGAATCTTTCTGCGCACGCCGATGCTGATGAAATCATg          <  2:162275/148‑1 (MQ=255)
                                                                                                                                      ctcCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGATGGAGAATCTTTCTGCGCACGCCGATGCTGATGAAATCATg          <  5:242620/148‑1 (MQ=255)
                                                                                                                                         cGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGATGGAGAATCTTTCTGCGCACGCCGATGCTGATGAAATCATgg         <  7:191014/146‑1 (MQ=255)
                                                                                                                                         cGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGATGGAGAATCTTTCTGCGCACGCCGATGCTGATGAAATCATGGaa       <  8:166826/148‑1 (MQ=255)
                                                                                                                                         cGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGATGGAGAATCTTTCTGCGCACGCCGATGCTGATGAAATCATGGaa       <  7:133602/148‑1 (MQ=255)
                                                                                                                                         cGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGATGGAGAATCTTTCTGCGCACGCCGATGCTGATGAAATCATGGAAt      >  6:250965/1‑149 (MQ=255)
                                                                                                                                         cGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGATGGAGAATCTTTCTGCGCACGCCGATGCTGATGAAATCATGGAAt      >  8:115933/1‑149 (MQ=255)
                                                                                                                                         cGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGATGGAGAATCTTTCTGCGCACGCCGATGCTGATGAAATCATGGAAt      >  1:336730/1‑149 (MQ=255)
                                                                                                                                          ggCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGATGGAGAATCTTTCTGCGCACGCCGATGCTGATGAAATCATGGAATg     >  2:146671/1‑149 (MQ=255)
                                                                                                                                          ggCTTGAAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGATGGAGAATCTTTCTGCGCACGCCGATGCTGATGAAATCATGGAATg     <  2:272918/149‑1 (MQ=255)
                                                                                                                                           gCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGATGGAGAATCTTTCTGCGCACGCCGATGCTGATGAAATCATGGAATgg    <  2:39416/149‑1 (MQ=255)
                                                                                                                                             tttCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACGGTGTGCGGTTGCACGGCGTGGATGTGCCTATCCGTGCAGAGGTGCATGCGATGGAGAATCTTTCTGCGCCCGCCGATGCTGATGAAATCATGGAATGGtt  >  4:313526/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
GTCGATGAGTCCAGGCACCTGGACCAGCTGCGCGAGCCCGCAGTGATCATCGCCGCCAGCGGCATGGCGACCGGCGGACGGGTGCTTCACCACCTGAAGGCACTGGCCCCAAACCCGCGCAACAGCATTCTGTTCTCCGGCTTTCAAGCCGGGGGGACGCGGGGTGCCGATATCGTCGCGGGCGCCCACAGTGTGCGGTTGCACGGCGAGGATGTGCCTATCCGTGCACAGGTGCATGCGATGGAGAATCTTTCTGCGCACGCCGATGCTGATGAAATCATGGAATGGTT  >  NC_002947/3726272‑3726561

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