Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,417,328 Δ1 bp 8.1% coding (576/759 nt) PP_1239 → metallo‑beta‑lactamase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,417,3280T. 8.1% 222.2 / 14.4 74coding (576/759 nt)PP_1239metallo‑beta‑lactamase family protein
Reads supporting (aligned to +/- strand):  ref base T (32/35);  new base . (2/4);  total (35/39)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.02e-01

GCGCTGTTGCTGGAGCGTTACCAGGGCCTGGATGCACTGCTGATCGAGGCTAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGTGGCGAAT  >  NC_002947/1417182‑1417475
                                                                                                                                                   |                                                                                                                                                    
gcgcTGTTGCTGGAGCGTTACCAGGGCCTGGATGCACTGCTGATCGAGGCTAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTgc                                                                                                                                                    >  2:134692/1‑149 (MQ=255)
  gcTGTTGCTGGAGCGTTACCAGGGCCTGGATGCACTGCTGATCGAGGCTAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTgcgc                                                                                                                                                  <  6:263254/149‑1 (MQ=255)
       tGCTGGAGCGTTACCAGGGCCTGGATGCACTGCTGATGGAGGCTAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTg                                                                                                                                            >  7:193818/1‑149 (MQ=255)
        gCTGGAGCGTTACCAGGGCCTGGATGCACTGCTGATCGAGGCTAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTg                                                                                                                                            <  1:134692/148‑1 (MQ=255)
        gCTGGAGCGGTACCAGGGCCTGGATGCACTGCTGATCGAGGCTAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTg                                                                                                                                            <  4:135983/148‑1 (MQ=255)
          tGGAGCGTTACCAGGGCCTGGATGCACTGCTGATCGAGGCTAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTg                                                                                                                                         >  7:379653/1‑149 (MQ=255)
          tGGAGCGTTACCAGGGCCTGGATGCACTGCTGATCGAGGCTAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTg                                                                                                                                         >  1:95310/1‑149 (MQ=255)
            gAGCGTTACCTGGGCCTGGATGCACTGCTGATCGAGGCTAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCGGAGGCAGCGGGTGGGTGGCATGCATGGGCATTTCAACAACCACCAGGCGGGGCGGCTGGTGCa                                                                                                                                       >  8:296691/1‑149 (MQ=255)
            gAGCGTTACCAGGGCCTGGATGCACTGCTGATCGAGGCTAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTg                                                                                                                                         <  1:11408/147‑1 (MQ=255)
            gAGCGTTACCAGGGCCTGGATGCACTGCTGATCGAGGCTAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCGCGCAGCCTGGTg                                                                                                                                         <  2:334497/147‑1 (MQ=255)
            gAGCGTTACCAGGGCCTGGATGCACTGCTGATCGAGGCTAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCGCGCAGCCTGGTg                                                                                                                                         <  7:51466/147‑1 (MQ=255)
             aGCGTTACCAGGGCCTGGATGCACTGCTGATCGAGGCTAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCAc                                                                                                                                      <  2:295246/149‑1 (MQ=255)
             aGCGTTACCAGGGCCTGGATGCACTGCTGATCGAGGCTAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCGCGCAGCCTGGTggtt                                                                                                                                      <  5:110975/149‑4 (MQ=255)
                 ttACCAGGGCCTGGATGCACTGCTGATCGAGGCTAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGt                                                                                                                                  >  5:420098/1‑149 (MQ=255)
                 ttACCAGGGCCTGGATGCACTGCTGATCGAGGCTAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGt                                                                                                                                  <  8:116337/149‑1 (MQ=255)
                 ttACCAGGGCCTGGATGCACTGCTGATCGAGGCTAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGt                                                                                                                                  <  8:122480/149‑1 (MQ=255)
                    ccAGGGCCTGGATGCACTGCTGATCGAGGCTAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGa                                                                                                                                    >  4:424351/1‑144 (MQ=255)
                    ccAGGGCCTGGATGCACTGCTGATCGAGGCTAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGa                                                                                                                                    <  3:424351/144‑1 (MQ=255)
                     cAGTGCCTGGATGCACTGCTGATCGAGGCTAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACTACCACCAGGCTGCGCGCCTGGTGCACGAGTTggg                                                                                                                              <  5:292915/149‑1 (MQ=255)
                       gggCCTGGATGCACTGCTGATCGAGGCTAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCt                                                                                                                            <  4:221834/149‑1 (MQ=255)
                        ggCCTGGATGCACTGCTGATCGAGGCTAACCACTGTCGCGACCTGGTGGCAGGTGGTCACTACGCGGCGTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACTACGACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTg                                                                                                                           <  7:168610/149‑1 (MQ=255)
                           cTGGATGCACTGCTGATCGAGGCTAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGaa                                                                                                                        >  2:317656/1‑149 (MQ=255)
                              gATGCACTGCTGATCGAGGCTAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAACGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGaa                                                                                                                        <  1:114797/147‑1 (MQ=255)
                              gATGCACTGCTGATCGAGGCTAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGaacaa                                                                                                                     <  2:95310/149‑1 (MQ=255)
                               aTGCACTGCTGATCGAGGCTAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGaacaac                                                                                                                    <  5:245521/149‑1 (MQ=255)
                                tGCACTGCTGATCGAGGCTAATCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAAcc                                                                                                                   >  7:76066/1‑149 (MQ=255)
                                  cACTGCTGATCGAGGCTAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTg                                                                                                                 <  3:120682/149‑1 (MQ=255)
                                   aCTGCTGATCGAGGCTAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTg                                                                                                                 <  8:345635/148‑1 (MQ=255)
                                   aCTGCTGATCGAGGCTAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTg                                                                                                                 >  7:345635/1‑148 (MQ=255)
                                         gATCGAGGCTAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACAcc                                                                                                          >  2:327506/1‑149 (MQ=255)
                                          aTCGAGGCTAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCt                                                                                                         >  5:292695/1‑149 (MQ=255)
                                           tCGAGGCTAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTg                                                                                                        >  2:319417/1‑149 (MQ=255)
                                            cGAGGCTAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTg                                                                                                        <  8:379653/148‑1 (MQ=255)
                                                    aCCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGcc                                                                                               >  8:39125/1‑149 (MQ=255)
                                                    aCCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCGCGCAGCCTGGTggt                                                                                                                                       <  3:422666/109‑3 (MQ=255)
                                                    aCCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCGCGCAGCCTGGTggt                                                                                                                                       >  4:422666/1‑107 (MQ=255)
                                                      cACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCa                                                                                             >  2:55110/1‑149 (MQ=255)
                                                       aCTGTCGCGACCTGCTGGCAAGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGTGTGGCATGCAATGGCATTTTAACAACCACCAGGCTGCGCGCCTGGTGAACGAGTTTGGCTGTAACAACCTGCAACACCTGGTGCTGGCCCAc                                                                                            <  6:439906/149‑1 (MQ=255)
                                                           tCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCaa                                                                                                              <  7:75715/127‑1 (MQ=255)
                                                           tCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCaa                                                                                                              >  8:75715/1‑127 (MQ=255)
                                                             gcgACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTcag                                                                                       >  3:155354/1‑148 (MQ=255)
                                                              cgACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCt                                                                                                   <  5:49725/135‑1 (MQ=255)
                                                              cgACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCt                                                                                                   >  6:49725/1‑135 (MQ=255)
                                                                 ccTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAgcagca                                                                                  <  6:420098/149‑1 (MQ=255)
                                                                  ctgctgGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCATCAGGCTGCGCGCCTGGt                                                                                                                                          >  4:389032/1‑92 (MQ=255)
                                                                  ctgctgGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCATCAGGCTGCGCGCCTGGt                                                                                                                                          <  3:389032/92‑1 (MQ=255)
                                                                     ctgGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAg                                                                                <  3:50552/147‑1 (MQ=255)
                                                                     ctgGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGaa                                                                              <  5:135211/149‑1 (MQ=255)
                                                                        gCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGaacaa                                                                           >  1:262342/1‑149 (MQ=255)
                                                                        gCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGaac                                                                             <  1:319417/147‑1 (MQ=255)
                                                                               gTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGTCCACCTCAGCAGCAAGAACAACCTGCca                                                                    >  3:187015/1‑149 (MQ=255)
                                                                                 cACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCcaca                                                                  >  1:5310/1‑149 (MQ=255)
                                                                                     aCCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTg                                                              >  1:410622/1‑149 (MQ=255)
                                                                                       ccGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGc                                                            >  1:87643/1‑149 (MQ=255)
                                                                                          gCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCg                                                         >  6:384493/1‑149 (MQ=255)
                                                                                           cTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGc                                                        <  6:292695/149‑1 (MQ=255)
                                                                                                    agcagcGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCttt                                               <  1:55110/149‑1 (MQ=255)
                                                                                                       agcGTGTGGGTGGCATGCAAGGGCATTTGAACAACCAGCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCGCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTc                                            <  8:81918/149‑1 (MQ=255)
                                                                                                          gggtgggtGTCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGATCAACCTGCTACACCTGGTGCTGGCCCACCTCAGCAGCTAGAACAACCTGCCAGACTTGGCGCGCCAGTGCTTTGTCGac                                         <  8:193818/149‑1 (MQ=255)
                                                                                                                gtggCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTc                                   <  2:5310/149‑1 (MQ=255)
                                                                                                                   gCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCgg                                 <  3:131730/148‑1 (MQ=255)
                                                                                                                         aaGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGAc                          >  6:388493/1‑149 (MQ=255)
                                                                                                                         aaGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGAc                          <  1:401887/149‑1 (MQ=255)
                                                                                                                                aTTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGAc                    >  1:403168/1‑149 (MQ=255)
                                                                                                                                 tttGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGTACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGACCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACt                   >  1:95858/1‑149 (MQ=255)
                                                                                                                                 tttGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACt                   >  7:326565/1‑149 (MQ=255)
                                                                                                                                      acaacCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCt                                                                        >  5:422259/1‑91 (MQ=255)
                                                                                                                                      acaacCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCt                                                                        <  6:422259/91‑1 (MQ=255)
                                                                                                                                       caacCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTc              >  3:104923/1‑148 (MQ=255)
                                                                                                                                         tgcaccaGGCGCGCAGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAg           >  1:213/3‑149 (MQ=255)
                                                                                                                                          ccaccaGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAgg          >  6:114788/1‑149 (MQ=255)
                                                                                                                                           caccaGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGt         >  5:189421/1‑149 (MQ=255)
                                                                                                                                              caGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGTGGc      <  4:155354/149‑1 (MQ=255)
                                                                                                                                                ggCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGTGGCGa    <  7:64178/149‑1 (MQ=255)
                                                                                                                                                   tGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGTGGCGAAt  >  7:319708/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCGCTGTTGCTGGAGCGTTACCAGGGCCTGGATGCACTGCTGATCGAGGCTAACCACTGTCGCGACCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGTGGCGAAT  >  NC_002947/1417182‑1417475

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: