Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,416,868 T→A 8.6% L39Q (CTG→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,8680TA 8.6% 195.1 / 7.9 70L39Q (CTG→CAG) PP_1239metallo‑beta‑lactamase family protein
Reads supporting (aligned to +/- strand):  ref base T (30/34);  new base A (4/2);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 4.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.80e-01

TGAGCGCTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGC‑ATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCAT  >  NC_002947/1416724‑1417009
                                                                                                                                                |                                                                                                                                              
tGAGCGCTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGgcg                                                                                                                                            >  2:72380/1‑149 (MQ=255)
  aGCGCTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCgct                                                                                                                                          <  4:309297/149‑1 (MQ=255)
  aGCGCTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCAGCGCCAGgcgc                                                                                                                                           <  5:191697/148‑1 (MQ=255)
    cgcTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAGGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCgctgc                                                                                                                                        >  8:312021/1‑149 (MQ=255)
    cgcTATGCAGGACCACCTGGGCGAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCgctgc                                                                                                                                        <  8:163266/149‑1 (MQ=255)
       tATTCAGGACCCCCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCg                                                                                                                                     >  4:283950/1‑149 (MQ=255)
       tATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCAGCGCCAGGCGCCGCTc                                                                                                                                      >  8:324987/1‑148 (MQ=255)
         ttCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCAGCGCCAGGCGCCGCTc                                                                                                                                      <  7:324987/146‑1 (MQ=255)
         ttCAGGACCACCTGCGCTAAGTGCGCTTCACGGTACTCGGAAGCGGCAGCCAAGGGAACGCCACGATGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTCCTCggc                                                                                                                                   >  2:168534/1‑148 (MQ=255)
             ggACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGttt                                                                                                                               >  3:196908/1‑149 (MQ=255)
              gACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTc                                                                                                                              >  5:256116/1‑149 (MQ=255)
               accaccTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCt                                                                                                                             <  5:341179/149‑1 (MQ=255)
                ccaccTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTg                                                                                                                            <  3:75600/149‑1 (MQ=255)
                    cTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCATCCAAGGGAACGCCACGCTGATCGCAAGTGGTAACACGGTCATCCTCGTAGATTGCGGCTTTTCACTGAGTGAAACCTAGAGGCGCATGGCACTGCTCGGGATTTCTGcagc                                                                                                                        >  6:290430/1‑149 (MQ=255)
                    cTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGcagc                                                                                                                        >  4:121159/1‑149 (MQ=255)
                            aGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGcagctcagct                                                                                                                  <  6:294025/147‑1 (MQ=255)
                             gTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAAAGGCACGCTGATCGCCAGTGGTGACACGTTCATCATGGTCGAGTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGctcagctcag                                                                                                               >  6:123707/1‑149 (MQ=255)
                             gTGCGCTTCGCGGGACGCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCGGTGGGGACACGTTCATCCGGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAACGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGctcagctcag                                                                                                               <  7:179908/149‑1 (MQ=255)
                              tGCGCTTCGCGGTACTCGGAAGCGGCAGCCAGGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTcagt                                                                                                              <  7:312021/149‑1 (MQ=255)
                              tGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTcagt                                                                                                              <  2:279175/149‑1 (MQ=255)
                                  cTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTcagtgcag                                                                                                          >  7:198771/1‑149 (MQ=255)
                                       cgGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTg                                                                                                     <  2:359913/149‑1 (MQ=255)
                                            cTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCa                                                                                                 <  2:160375/148‑1 (MQ=255)
                                             tCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCAcc                                                                                               >  2:267695/1‑149 (MQ=255)
                                               ggAAGCGGCAGCCAAGGGAACGGCAAGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCAcc                                                                                               <  6:174297/147‑1 (MQ=255)
                                                gAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCAc                                                                                            >  3:298198/1‑149 (MQ=255)
                                                gAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCAc                                                                                            <  3:121159/149‑1 (MQ=255)
                                                         gCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGcc                                                                                   <  4:196908/149‑1 (MQ=255)
                                                           cAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAGCATGCCGa                                                                                 >  8:159553/1‑149 (MQ=255)
                                                               ggAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCAc                                                                             >  3:87787/1‑149 (MQ=255)
                                                                 aaCGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGGCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGt                                                                           <  2:21910/149‑1 (MQ=255)
                                                                       aCGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGC‑ATgg                                                                    <  6:256116/149‑1 (MQ=255)
                                                                            gATCGCCAGTGGTGACACGTTCATCCTGGTCGATAGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCCCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCCCATGC‑ACGTGGTCg                                                               >  1:50849/1‑149 (MQ=255)
                                                                             aTCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGC‑ATGGGGTCgg                                                              >  2:3104/1‑149 (MQ=255)
                                                                              tCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCCCATGC‑ACGTGGTCGGc                                                             <  8:38398/149‑1 (MQ=255)
                                                                                gCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTc                                                                                                  <  2:129650/111‑1 (MQ=255)
                                                                                gCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGATTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTCCAGTGCTGGTc                                                                                                  >  1:129650/1‑111 (MQ=255)
                                                                                 ccAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGC‑ATGGGGTCGGCTTg                                                          >  2:48684/1‑149 (MQ=255)
                                                                                   aGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGC‑ATGGGGTCGGCTTGCt                                                        <  4:298198/149‑1 (MQ=255)
                                                                                   aGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGC‑ATGGGGTCGGCTTGCt                                                        <  7:159553/149‑1 (MQ=255)
                                                                                       gtgACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGC‑ATGGGGTCGGCTTGCTGt                                                      <  8:338396/147‑1 (MQ=255)
                                                                                       gtgACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGC‑ATGGGGTCGGCTTGCTGt                                                      >  7:338396/1‑147 (MQ=255)
                                                                                       gtgACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGC‑ATGGGGTCGGCTTGCTGTc                                                     <  1:205204/148‑1 (MQ=255)
                                                                                       gtgACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGC‑ATGGGGTCGGCTTGCTGTCa                                                    >  3:102432/1‑149 (MQ=255)
                                                                                            acGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGC‑ATGGGGTCGGCTTGCTGTCa                                                    >  3:368428/1‑144 (MQ=255)
                                                                                            acGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGC‑ATGGGGTCGGCTTGCTGTCa                                                    <  4:368428/144‑1 (MQ=255)
                                                                                              gTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGC‑ATGGGGTCGGCTTGCt                                                        <  7:213940/138‑1 (MQ=255)
                                                                                              gTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGC‑ATGGGGTCGGCTTGCt                                                        >  8:213940/1‑138 (MQ=255)
                                                                                              gTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGC‑ATGGGGTCGGCTTGCTGTCACGgcgc                                              >  2:327807/1‑148 (MQ=255)
                                                                                                   tCCTGGTGGATTGCGGCTTTTCCCTGCGTGAAAGCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGC‑ATGGGGTCGGCTTGCTGTCACGGCGCTACAAt                                        <  1:168534/149‑1 (MQ=255)
                                                                                                   tCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGC‑ATGGGGTCGGCTTGCTGTCACGGCGCTACAAt                                        <  4:196754/149‑1 (MQ=255)
                                                                                                        gTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACTAACATGCCGACCACGTGC‑ATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTAcc                                   <  3:283950/149‑1 (MQ=255)
                                                                                                             ttGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGC‑ATGGGGTc                                                                >  7:151495/1‑115 (MQ=255)
                                                                                                             ttGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGC‑ATGGGGTc                                                                <  8:151495/115‑1 (MQ=255)
                                                                                                               gCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGC‑ATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTAc                            >  2:342736/1‑149 (MQ=255)
                                                                                                                cGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGCAATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTAc                            <  1:72380/149‑1 (MQ=255)
                                                                                                                cGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGC‑ATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTAcc                           >  3:330806/1‑149 (MQ=255)
                                                                                                                cGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGC‑ATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTAc                            <  1:3104/148‑1 (MQ=255)
                                                                                                                  gCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGC‑ATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCt                          >  2:287686/1‑148 (MQ=255)
                                                                                                                      ttCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGC‑ATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGc                      <  8:198771/148‑1 (MQ=255)
                                                                                                                        cccTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGC‑ATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGcc                     >  5:135608/1‑147 (MQ=255)
                                                                                                                         ccTGCGTGAAACCGAGCCGCGCCTGGCGCTGCTCGGTGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGAACACGTGC‑ATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAg                  <  5:290430/149‑1 (MQ=255)
                                                                                                                              gaaaccccGAGCCGCGCCAGGCGCCGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCCCCCACGCACATGCCGCCCACGTGC‑ATGGGGTCGGCTGGCTGTCACGGCGCTACAATGTACCGGTCTACCTCCGCCAAGGGACt             >  5:41686/7‑149 (MQ=255)
                                                                                                                                gAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGC‑ATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACttt           >  4:309481/1‑149 (MQ=255)
                                                                                                                                   cccGAGCAGCGCCAGGCGCCGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGC‑ATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTgcg        >  7:282419/2‑149 (MQ=255)
                                                                                                                                   cccGAGCAGCGCCAGGCGCCGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGC‑ATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTgcg        >  4:340783/2‑149 (MQ=255)
                                                                                                                                     cGAGCGGCGCCTGGCGCTGCGCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGC‑ATGGGGGCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGcgcg      <  5:123707/149‑1 (MQ=255)
                                                                                                                                         cggcgCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGC‑ATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCAt  <  1:48684/149‑1 (MQ=255)
                                                                                                                                           gcgcCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGC‑ATGGGGTCgg                                                              <  7:56791/87‑1 (MQ=255)
                                                                                                                                           gcgcCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGC‑ATGGGGTCgg                                                              >  8:56791/1‑87 (MQ=255)
                                                                                                                                                |                                                                                                                                              
TGAGCGCTATTCAGGACCACCTGGGCTAAGTGCGCTTCGCGGTACTCGGAAGCGGCAGCCAAGGGAACGGCACGCTGATCGCCAGTGGTGACACGTTCATCCTGGTCGATTGCGGCTTTTCCCTGCGTGAAACCGAGCGGCGCCTGGCGCTGCTCGGGGTTTCTGCAGCTCAGCTCAGTGCAGTGCTGGTCACCCACGAACATGCCGACCACGTGC‑ATGGGGTCGGCTTGCTGTCACGGCGCTACAATGTACCGGTCTACCTCAGCCAAGGGACTTTGCGCGGCAT  >  NC_002947/1416724‑1417009

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑