Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,416,862 G→A 9.2% R37Q (CGG→CAG)  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,416,8620GA 9.2% 184.8 / 6.5 65R37Q (CGG→CAG) PP_1239metallo‑beta‑lactamase family protein
Reads supporting (aligned to +/- strand):  ref base G (29/30);  new base A (4/2);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

TGCTGAGCGCTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATG  >  NC_002947/1416721‑1417010
                                                                                                                                              |                                                                                                                                                      
tGCTGAGCGCTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCG‑‑CCTg                                                                                                                                               <  5:339970/149‑1 (MQ=255)
            ttCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCAGCG‑‑CCAGGCGCCGCTCg                                                                                                                                     <  5:12788/147‑1 (MQ=255)
                ggACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGttt                                                                                                                               <  6:140843/149‑1 (MQ=255)
                 gACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTc                                                                                                                              >  4:40050/1‑149 (MQ=255)
                 gACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTc                                                                                                                              <  2:317173/149‑1 (MQ=255)
                  accaccTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCt                                                                                                                             >  3:105732/1‑149 (MQ=255)
                   ccaccTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTg                                                                                                                            <  5:106918/149‑1 (MQ=255)
                     accTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGca                                                                                                                          <  7:157685/149‑1 (MQ=255)
                        tGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGcagct                                                                                                                       >  7:165706/1‑149 (MQ=255)
                         gggCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGcagct                                                                                                                       >  6:14500/1‑148 (MQ=255)
                               aGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGCAgctcagctca                                                                                                                <  8:165706/149‑1 (MQ=255)
                                      ttCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCTCTGCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTcagtgcagt                                                                                                         >  4:388570/1‑149 (MQ=255)
                                      ttCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTcagtgcagt                                                                                                         >  7:101236/1‑149 (MQ=255)
                                          cgGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTg                                                                                                     <  4:68369/149‑1 (MQ=255)
                                            gTACTCGGAAGCGGCAGCCAAGGGATCGGCACTCTGATCGCCAGTGGTGACACGTTCATGCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGt                                                                                                   <  6:128010/149‑1 (MQ=255)
                                               cTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCa                                                                                                 <  6:158157/148‑1 (MQ=255)
                                               cTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCAc                                                                                                <  1:340083/149‑1 (MQ=255)
                                                tCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTg                                                                                                     <  7:51686/143‑1 (MQ=255)
                                                tCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTg                                                                                                     >  8:51686/1‑143 (MQ=255)
                                                      gCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTgg                                                                                                    <  3:292076/138‑1 (MQ=255)
                                                      gCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTgg                                                                                                    >  4:292076/1‑138 (MQ=255)
                                                       cGGCAGCCAAGGGAACGGCACGCTGATCGCCAATGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGaa                                                                                         >  1:250776/1‑148 (MQ=255)
                                                               aaGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCTCTGCGTGTAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGAc                                                                                <  3:388570/149‑1 (MQ=255)
                                                               aaGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGAc                                                                                <  5:14500/149‑1 (MQ=255)
                                                                  ggAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGc                                                                                    >  7:57621/1‑142 (MQ=255)
                                                                  ggAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGc                                                                                    <  8:57621/142‑1 (MQ=255)
                                                                     aCGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTg                                                                          >  5:119459/1‑149 (MQ=255)
                                                                       ggCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCa                                                                        >  7:95033/1‑149 (MQ=255)
                                                                          aCGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGc                                                                                                                           <  5:53131/95‑1 (MQ=255)
                                                                          aCGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGc                                                                                                                           >  6:53131/1‑95 (MQ=255)
                                                                              tGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTc                                                                 >  3:172696/1‑149 (MQ=255)
                                                                                aTCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCCCATGCACGTGGTCg                                                                <  8:101236/148‑1 (MQ=255)
                                                                                aTCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTc                                                                 <  6:118853/147‑1 (MQ=255)
                                                                                    ccAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCACTGCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTg                                                           >  6:216572/1‑149 (MQ=255)
                                                                                     cAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGAAGCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCtt                                                            >  1:322896/1‑149 (MQ=255)
                                                                                           tgACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCAc                                                    <  6:324119/149‑1 (MQ=255)
                                                                                                 gTTCATCCTGGTCGTTTGCGGCTTTTCCCTGCGTTAAACCTAGCGGCG‑‑CCTGGCGCTGGTGGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCt                                                         <  2:198545/138‑1 (MQ=255)
                                                                                                 gTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCt                                                         >  1:198545/1‑138 (MQ=255)
                                                                                                 gTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGgcgc                                               >  4:57283/1‑148 (MQ=255)
                                                                                                 gTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGgcgc                                               >  3:385927/1‑148 (MQ=255)
                                                                                                 gTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCt                                              <  8:95033/149‑1 (MQ=255)
                                                                                                     aTCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACaa                                          >  1:321423/1‑149 (MQ=255)
                                                                                                       ccTGGTCGATTGCGGCTTTTCCCTGCCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTc                                                      >  1:212551/1‑136 (MQ=255)
                                                                                                       ccTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTc                                                      <  2:212551/135‑1 (MQ=255)
                                                                                                         tGGTCGATTGCGGCTTTTCCCTTCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGc                                                          <  8:127930/129‑1 (MQ=255)
                                                                                                         tGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGc                                                          >  7:127930/1‑129 (MQ=255)
                                                                                                         tGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTa                                      >  1:388743/1‑149 (MQ=255)
                                                                                                           gTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTAcc                                    >  7:3516/1‑149 (MQ=255)
                                                                                                           gTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTAc                                     <  3:57283/148‑1 (MQ=255)
                                                                                                            tCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCg                                   <  5:216572/149‑1 (MQ=255)
                                                                                                              gATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCgg                                  >  4:93722/1‑148 (MQ=255)
                                                                                                                  gCGGCTTTTCCCTGCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTa                                             <  1:392112/133‑1 (MQ=255)
                                                                                                                  gCGGCTTTTCCCTGCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTa                                             >  2:392112/1‑133 (MQ=255)
                                                                                                                         ttCCCTGCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGc                       >  5:7863/1‑148 (MQ=255)
                                                                                                                            ccTGCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAg                   >  7:235869/1‑149 (MQ=255)
                                                                                                                             cTGCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCTGTCTACCTCAGCCAAgg                  <  4:172696/149‑1 (MQ=255)
                                                                                                                                  gTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCCCATGCACGTGGTCGGc                                                              <  6:319844/101‑1 (MQ=255)
                                                                                                                                  gTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCCCATGCACGTGGTCGGc                                                              >  5:319844/1‑101 (MQ=255)
                                                                                                                                   tGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACtt             <  8:3516/149‑1 (MQ=255)
                                                                                                                                     accccGAGCAGCG‑‑CCAGGCGCCGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATACCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCa                         >  6:129577/4‑135 (MQ=255)
                                                                                                                                     accccGAGCAGCG‑‑CCAGGCGCCGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATACCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCa                         <  5:129577/132‑1 (MQ=255)
                                                                                                                                       cccGAGCAGCG‑‑CCAGGCGCCGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTgcg         >  8:34529/2‑149 (MQ=255)
                                                                                                                                       cccGAGCAGCG‑‑CCAGGCGCCGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTgcg         >  6:190367/2‑149 (MQ=255)
                                                                                                                                        ccGAGCAGCG‑‑CCAGGCGCCGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTgcgc        >  1:43060/1‑149 (MQ=255)
                                                                                                                                         cGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGCAGCTCCGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGcgcg       <  3:93722/149‑1 (MQ=255)
                                                                                                                                            gcggcg‑‑CCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACttt            <  7:120237/141‑1 (MQ=255)
                                                                                                                                            gcggcg‑‑CCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACttt            >  8:120237/1‑141 (MQ=255)
                                                                                                                                              ggcg‑‑CCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATg  <  6:119459/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                      
TGCTGAGCGCTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCG‑‑CCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCATG  >  NC_002947/1416721‑1417010

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: