Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,211,318 A→T 6.4% F211Y (TTC→TAC)  PP_4592 ← metal dependent hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,211,3180AT 6.4% 182.8 / 2.9 62F211Y (TTC→TAC) PP_4592metal dependent hydrolase
Reads supporting (aligned to +/- strand):  ref base A (29/29);  new base T (2/2);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

GAGATTGGGCGTGGTATCGCGGTCATAGCCCAACACGCCGAAATGGTCTGCCCGGGCACTTTCCCCGACCACCAGCACGGTCAGCGACTTGCGCTCATGTTTTAGCCAG‑GCCGCATCACGCTTGGCATCTTCGCCATAATGCTGGAACGGGCGCGAGGCCGTACCGACACGTTCGTCGACATAACCAATGGAGGCCCCGACGATATTGCTGGGGGTGAGCATCAGACGCAGTTCATGGTGATTGCGAAACAGCGACGACAGCCCTTGATAATTGACCAACGCGACCGAGCCGA  >  NC_002947/5211173‑5211465
                                                                                                                                                  |                                                                                                                                                   
gagaTTGGGCGTGGTATCGCGGTCATAGCCCAACACGCCGAAATGGTCTGCCCGGGCACTTTCCCCGACCACCAGCACGGTCAGCGACTTGCGCTCATGTTTTAGCCAG‑GCCGCATCACGCTTGGCATCTTCGCCATAATGCTGGAACg                                                                                                                                                  >  7:197330/1‑149 (MQ=255)
gagaTTGGGCGTGGTATCGCGGTCATAGCCCAACACGCCGAAATGGTCTGCCCGGGCACTTTCCCCGACCACCAGCACGGTCAGCGACTTGCGCTCATGTTTTAGCCAG‑GCCGCATCACGCTTGGCATCTTCGCCATAATGCTGGAACg                                                                                                                                                  <  2:11672/149‑1 (MQ=255)
 agaTTGGGCGTGGTATCGCGGTCATAGCCCAACACGCCGAAATGGTCTGCCCGGGCACTTTCCCCGACCACCAGCACGGTCAGCGACTTGCGCTCATGTTTTAGCCAG‑GCCGCATCACGCTTGGCATCTTCGCCATAATGCTGGAACgg                                                                                                                                                 <  2:123700/149‑1 (MQ=255)
   aTTGGGCGTGGTATCGCGGTCATAGCCCAACACGCCGAAATGGTCTGCCCGGGCACTTTCCCCGACCACCAGCACGGTCAGCGACTTGCGCTCATGTTTTAGCCAG‑GCCGCATCACGCTTGGCATCTTCGCCATAATGCTGGAACGGgc                                                                                                                                               >  2:237597/1‑149 (MQ=255)
        gCGTGGTATCGCGGTCATAGCCCAACACGCCGAAATGGTCTGCCCGGGCACTTTCCCCGACCACCAGCACGGTCAGCGACTTGCGCTCATGTTTTAGCCAG‑GCCGCATCACGCTTGGCATCTTCGCCATAATGCTGGAACGGGCGCGAg                                                                                                                                          <  8:197330/149‑1 (MQ=255)
        gCGTGGTATCGCGGTCATAGCCCAACACGCCGAAATGGTCTGCCCGGGCACTATCCCCGACCACCAGCACGGTCAGCGACTTGCGCTCATGTTTTAGCCAG‑GCCGCATCACGCTTGGCATCTTCGCCATAATGCTGGAACGGGCGCGAg                                                                                                                                          >  2:280957/1‑149 (MQ=255)
                tCGCGGTCATAGCCCAACACGCCGAAATGGTCTGCCCGGGCACTTTCCCCGACCACCAGCACGGTCAGCGACTTGCGCTCATGTTTTAGCCAG‑GCCGCATCACGCTTGGCATCTTCGCCATAATGCTGGTACGGCCTCGCGCCCGttc                                                                                                                                   >  6:75799/1‑146 (MQ=255)
                  gcgGTCATAGCCCAACACGCCGAAATGGTCTGCCCGGGCACTTTCCCCGACCACCAGCACGGTCAGCGACTTGCGCTCATGTTTTAGCCAG‑GCCGCATCACGCTTGGCATCTTCGCCATAATGCTGGTACGGCCTCGCGCCCGTTcca                                                                                                                                 <  5:75799/148‑2 (MQ=255)
                  gcgGTCATAGCCCAACACGCCGAAATGGTCTGCCCGGGCACTTTCCCCGACCACCAGCACGGTCAGCGACTTGCGCTCATGTTTTAGCCAG‑GCCGCATCACGCTTGGCATCTTCGCCATAATGCTGGAACGGGCGCGAGGCCGTACCGa                                                                                                                                >  2:227761/1‑149 (MQ=255)
                  gcgGTCATAGCCCAACACGCCGAAATGGTCTGCCAGGGCACTTTCCCCGACCACCAGCACGGTCAGCGACTTGCGCTCATGTTTTAGCCAG‑GCCGCATCACGCTAGGCATCTTCGCCAGAATCCTGGAACGGGCGCGAGGCCGTACCGa                                                                                                                                >  8:322780/1‑149 (MQ=255)
                   cgGTCATAGCCCAACACGCCGAAATGGTCTGCCCGGGCACTTTCCCCGACCACCAGCACGGTCAGCGACTTGCGCTCATGTTTTAGCCAG‑GCCGCATCACGCTTGGCATCTTCGCCATAATGCTGGAACGGGCGCGAGGCCGTACCGac                                                                                                                               <  1:280957/149‑1 (MQ=255)
                   cgGTCATAGCCCAACACGCCGAAATGGTCTGCCCGGGCACTTTCCCCGACCACCAGCACGGTCAGCGACTTGCGCTCATGTTTTAGCCAG‑GCCGCATCACGCTTGGCATCTTCGCCATAATGCTGGAACGGGCGCGAGGCCGTACCGa                                                                                                                                >  2:75401/1‑148 (MQ=255)
                        aTAGCCCAACACGCCGAAATGGTCTGCCCGGGCACTTTCCCCGACCACCAGCACGGTCAGCGACTTGCGCTCATGTTTTAGCCAG‑GCCGCATCACGCTTGGCATCTTCGCCATAATGCTGGAACGGGCGCGAGGCCGTAcc                                                                                                                                  >  6:121638/1‑141 (MQ=255)
                        aTAGCCCAACACGCCGAAATGGTCTGCCCGGGCACTTTCCCCGACCACCAGCACGGTCAGCGACTTGCGCTCATGTTTTAGCCAG‑GCCGCATCACGCTTGGCATCTTCGCCATAATGCTGGAACGGGCGCGAGGCCGTAcc                                                                                                                                  <  5:121638/141‑1 (MQ=255)
                        aTAGCCCAACACGCCGAAATGGTCTGCCCGGGCACTTTCCCCGACCACCAGCACGGTCAGCGACTTGCGCTCATGTTTTAGCCAG‑GCCGCATCACGCTTGGCATCTTCGCCATAATGCTGGAACGGGCGCGAGGCCGTACCGACACGtt                                                                                                                          >  4:25979/1‑149 (MQ=255)
                         tAGCCCAACACGCCGAAATGGTCTGCCCGGGCACTTTCCCCGACCACCAGCACGGTCAGCGACTTGCGCTCATGTTTTAGCCAG‑GCCGCATCACGCTTGGCATCTTCGCCATAATGCTGGAACGGGCGCGAGGCCGTACCGACACGTtc                                                                                                                         <  8:31135/149‑1 (MQ=255)
                           gCCCAACAGGCCGAACTGGGCTGCCCGGGCACTTTCCCCGACCACCAGCTCGGTCAGCGACTTGCGCGGATGTTTTAGCCAG‑GCCGCATCACGCTTGGCATCTTCGCCATAATGCTGGAACGTGCGCGTGGCCGTACCGACACGTtcgt                                                                                                                       <  6:346626/149‑1 (MQ=255)
                             ccAACACGCCGAAATGGTCTGCCCGGGCACTTTCCCCGACCACCAGCACGGTCAGCGACTTGCGCTCATGTTTTAGCCAG‑GCCGCATCACGCTTGGCATCTTCGCCATAATGCTGGAACGGGCGCGAGGCCGTACCGACACGTtcgtcg                                                                                                                     <  3:25979/149‑1 (MQ=255)
                              cAACACGCCGAAATGGTCTGCCCGGGCACTTTCCCCGACCACCAGCACGGTCAGCGACTTGCGCTCATGTTTTAGCCAG‑GCCGCATCACGCTTGGCATCTTCGCCATAATGCTGGAACGGGCGCGAGGCCGTACCGACACGTTCGTCGa                                                                                                                    <  1:237597/149‑1 (MQ=255)
                               aaCACGCCGAAATGGTCTGCCCGGGCACTTTCCCCGACCACCAGCACGGTCAGCGACTTGCGCTCATGTTTTAGCCAG‑GCCGCATCACGCTTGGCATCTTCGCCATAATGCTGGAACGGGCGCGAGGCCGTACCGACACGTTCGTCGAc                                                                                                                   <  1:75401/149‑1 (MQ=255)
                                    gCCGAAATGGTCTGCCCGGGCACTTTCCCCGACCACCAGCACGGTCAGCGACTTGCGCTCATGTTTTAGCCAG‑GCCGCATCACGCTTGGCATCTTCGCCATAATGCTGGAACGGGCGCGAGGCCGTACCGACACGTTCGTCGACATAAc                                                                                                              >  5:205689/1‑149 (MQ=255)
                                      cGAAATGGTCTGCCCGGGCACTTTCCCCGACCACCAGCACGGTCAGCGACTTGCGCTCATGTTTTAGCCAG‑GCCGCATCACGCTTGGCATCTTCGCCATAATGCTGGAACGGGCGCGAGGCCGTACCGACACGTTCGTCGa                                                                                                                    <  1:61830/141‑1 (MQ=255)
                                      cGAAATGGTCTGCCCGGGCACTTTCCCCGACCACCAGCACGGTCAGCGACTTGCGCTCATGTTTTAGCCAG‑GCCGCATCACGCTTGGCATCTTCGCCATAATGCTGGAACGGGCGCGAGGCCGTACCGACACGTTCGTCGa                                                                                                                    >  2:61830/1‑141 (MQ=255)
                                        aaaTGGTCTGCCCGGGCACTTTCCCCGACCACCAGCACGGTCAGCGACTTGCGCTCATGTTTTAGCCAG‑GCCGCATCACGCTTGGCATCTTCGCCATAATGCTGGAACGGGCGCGAGGCCGTACCGACACGTTCGTCGACATAACCAAt                                                                                                          >  1:214580/1‑149 (MQ=255)
                                                   ccGGGCACTTTCCCCGACCACCAGCACGGTCAGCGACTTGCGCTCATGTTTTAGCCAG‑GCCGCATCACGCTTGGCATCTTCGCCATAATGCTGGAACGGGCGCGAGGCCGTACCGACACGTTCGTCGACATAACCAATGGAGGCCCcga                                                                                               >  6:71119/1‑149 (MQ=255)
                                                         aCTTTCCCCGACAACCAGAACGGTCAGCTACTTGGGCTCATTTTTTTGCCAGCGCAAGATCTTATATAATAACTTCGCCATAATGCTGGAACGGGCGCGAGGCCGTACCGACACGTTCGTCGACATAACCAATGGAGGCCCCGACGata                                                                                          <  1:111945/149‑1 (MQ=255)
                                                               cccGACCACCAGCACGGTCAGCGACTTGCGCTCATGTTTTAGCCAG‑GCCGCATCACGCTTGGCATCTTCGCCATAATGCTGGAACGGGCGCGAGGCCGTACCGACACGTTCGTCGACATAACCAATGGAGGCCCCGACGATATTGCTgg                                                                                   >  4:337798/1‑149 (MQ=255)
                                                                 cGACCACCAGCACGGTCAGCGACTTGCGCGCATGTTTTAGCCAG‑GCCGCATCACGCTTGGCATCTTCGCCATAATGCTGGAACGGGCGCGAGGCCGTACCGACACGTTCGTCGACATAACCAATGGAGGCCCCGACGATATTGCTgggg                                                                                 <  7:322780/149‑1 (MQ=255)
                                                                       ccaGCACGGTCAGCGACTTGCGCTCATGTTTTAGCCAG‑GCCGCATCACGCTTGGCATCTTCGCCATAATGCTGGAACGGGCGCGAGGCCGTACCGACACGTTCGTCGACATAACCAATGGAGGCCCCGACGATATTGCTGGGGGTGAg                                                                            >  7:305389/1‑148 (MQ=255)
                                                                       ccaGCACGGTCAGCGACTTGCGCTCATGTTTTAGCCAG‑GCCGCATCACGCTTGGCATCTTCGCCATAATGCTGGAACGGGCGCGAGGCCGTACCGACACGTTCGTCGACATAACCAATGGAGGCCCCGACGATATTGCTGGGGGTGAg                                                                            >  1:69353/1‑148 (MQ=255)
                                                                       ccaGCACGGTCAGCGACTTGCGCTCATGTTTTAGCCAG‑GCCGCATCACGCTTGGCATCTTCGCCATAATGCTGGAACGGGCGCGAGGCCGTACCGACACGTTCGTCGACATAACCAATGGAGGCCCCGACGATATTGCTGGGGGTGAg                                                                            <  2:69353/148‑1 (MQ=255)
                                                                           cACGGTCAGCGACTTGCGCTCATGTTTTAGCCAG‑GCCGCATCACGCTTGGCATCTTCGCCATAATGCTGGAACGGGCGCGAGGCCGTACCGACACGTTCGTCGACATAACCAATGGAGGCCCCGACGATATTGCTGGGGGTGAGCATCa                                                                       <  5:71119/149‑1 (MQ=255)
                                                                               gTCAGCGACTTGCGCTCATGTTTTAGCCAG‑GCCGCATCACGCTTGGCATCTTCGCCATAATGCTGGAACGGGCGCGAGGCCGTACCGACACGTTCGTCGCCATAACCAATGGAGGCCCCGACGCTATGGCGGGGGGTGATCATCAGACg                                                                   >  6:29526/1‑149 (MQ=255)
                                                                                  aGCGACTTGCGCTCATGTTTTAGCCAG‑GCCGCATCACGCTTGGCATCTTCGCCATAATGCTGGAACGGGCGCGAGGCCGTACCGACACGTTCGTCGACATAACCAATGGAGGCCCCGACGATATTGCTGGGGGTGAGCATCAGACGCAg                                                                >  4:153449/1‑149 (MQ=255)
                                                                                  aGCGACTTGCGCTCATGTTTTAGCCAG‑GCCGCATCACGCTTGGCATCTTCGCCATAATGCTGGAACGGGCGCGAGGCCGTACCGACACGTTCGTCGACATAACCAATGGAGGCCCCGACGATATTGCTGGGGGTGAGCATCAGACGCAg                                                                >  6:133801/1‑149 (MQ=255)
                                                                                          gcgcTCATGTTTTAGCCAG‑GCCGCATCACGCTTGGCATCTTCGCCATAATGCTGGAACGGGCGCGAGGCCGTACCGACACGTTCGTCGACATAACCAATGGAGGCCCCGACGATATTGCTGGGGGTGAGCATCAGACGCAGTTCATg                                                          <  6:205689/147‑1 (MQ=255)
                                                                                           cgcTCATGTTTTAGCCAG‑GCCGCATCACGCTTGGCATCTTCGCCATAATGCTGGAACGGGCGCGAGGCCGTACCGACACGTTCGTCGACATAACCAATGGAGGCCCCGACGATATTGCTGGGGGTGAGCATCAGACGCAGTTCATGGTg                                                       <  3:337109/149‑1 (MQ=255)
                                                                                             cTCATGTTTTAGCCAG‑GCCGCATCACGCTTGGCATCTTCGCCATAATGCTGGAACGGGCGCGAGGCCGTACCGACACGTTCGTCGACATAACCAATGGAGGCCCCGACGATATTGCTGGGGGTGAGCATCAGACGCAGTTCATGGTGAt                                                     >  5:238573/1‑149 (MQ=255)
                                                                                              tCATGTTTTAGCCAG‑GCCGCATCACGCTTGGCATCTTCGCCATAATGCTGGAACGGGCGCGAGGCCGTACCGACACGTTCGTCGACATAACCAATGGAGGCCCCGACGATATTGCTGGGGGTGAGCATCAGACGCAGTTCATGGTGAtt                                                    >  4:292256/1‑149 (MQ=255)
                                                                                                aTGTTTTAGCCAG‑GCCGCATCACGCTTGGCATCTTCGCCATAATGCTGGAACGGGCGCGAGGCCGTACCGACACGTTCGTCGACATAACCAATGGAGGCCCCGACGATATTGCTGGGTGTGAGCATCAGACGCAGTTCATGGTGAtt                                                    <  1:227761/147‑1 (MQ=255)
                                                                                                aTGTTTTAGCCAG‑GCCGCATCACGCTTGGCATCTTCGCCATAATGCTGGAACGGGCGCGAGGCCGTACCGACACGTTCGTCGACATAACCAATGGAGGCCCCGACGATATTGCTGGGGGTGAGCATCAGACGCAGTTCATGGTGATTGc                                                  >  4:284443/1‑149 (MQ=255)
                                                                                                aTGTTTTAGCCAG‑GCCGCATCACGCTTGGCATCTTCGCCATAATGCTGGAACGGGAGCGAGGCCGTACCTACCCGTTCGTCGACATAACCAATGGAGGCCCCGACGATATTGCTTGGGGTGAGCCTCAGACGCAGTTAATGGTGAtt                                                    >  5:289947/1‑147 (MQ=255)
                                                                                                   ttttAGCCAG‑GCCGCATCACGCTTGGCATCTTCGCCATAATGCTGGAACGGGCGCGAGGCCGTACCGACACGTTCGTCGACATAACCAATGGAGGCCCCGACGATATTGCTGGGGGTGAGCATCAGACGCAGTTCATGGTGATTGCGa                                                >  2:284839/1‑148 (MQ=255)
                                                                                                     ttAGCCAG‑GCCGCATCACGCTTGGCATCTTCGCCATAATGCTGGAACGGGCGCGAGGCCGTACCGACACGTTCGTCGACATAACCAATGGAGGCCCCGACGATATTGCTGGGGGTGAGCATCAGACGCAGTTCATGGTGATTGCGaaa                                              <  1:284839/148‑1 (MQ=255)
                                                                                                        gCCAG‑GCCGCATCACGCTTGGCATCTTCGCCATAATGCTGGAACGGGCGCGCGGCCGTACCGACACGTTCGTCGACATAACCAATGGAGGCCCCGACGATATTGCTGGGGGTGAGCATCAGACGCAGTTCATGGTGATTGCGAAACAGc                                          >  5:241994/1‑149 (MQ=255)
                                                                                                          cAG‑GCCGCATCACGCTTGGCATCTTCGCCATAATGCTGGTACGGGCGCGAGGCCGTACCGACACGTTCGTCGACATAACCAATGGAGGCCCCGACGATATTGCTGGGGGTGAGCATCAGACGCAGTTCATGGTGATTGCGAAACAGcga                                        >  4:111688/1‑149 (MQ=255)
                                                                                                          cAG‑GCCGCATCACGCTTGGCATCTTCGCCATAATGCTGGAACGGGCGCGAGGCCGTACCGACACGTTCGTCGACATAACCAATGGAGGCCCCGACGATATTGCTGGGGGTGAGCATCAGACGCAGTTCATGGTGATTGCGAAACAGcga                                        >  2:327769/1‑149 (MQ=255)
                                                                                                           aG‑GCCGCATCACGCTTGGCATCTTCGCCATAATGCTGGAACGGGCGCGAGGCCGTACCGACACGTTCGTCGACATAACCAATGGAGGCCCCGACGATATTGCTGGGGGTGAGCATCAGACGCAGTTCATGGTGATTGCGAAACAGcgac                                       <  3:337798/149‑1 (MQ=255)
                                                                                                           aG‑GCAGCATCACGCTTGGCATCTTCGCCATAATGCTTGAACGGGCGCTAGGCCGTACCGACACGTTAGTCGACATAACCAATGGAGGCCCCGACGATATTGCTGGGGGTGAGCATCAGACGCAGTTCATGGTGATTGCGAAACAGcgac                                       <  8:305389/149‑1 (MQ=255)
                                                                                                                      aCGCTTGGCATCTTCGCCATAATGCTGGAACGGGCGCGAGGCCGTACCGACACGTTCGTCGACATAACCAATGGAGGCCCCGACGATATTGCTGGGGGTGAGCATCAGACGCAGTTCATGGTGATTGCGAAACAGCGACGACAGCCCt                              <  3:292256/148‑1 (MQ=255)
                                                                                                                        gCTTGGCATCTTCGCCATAATGCTGGTACGGGCGCGAGGCCGTACCGACACGTTCGTCGACATAACCAATGGAGGCCCCGACGATATTGCTGGGGGTGAGCATCAGACGCAGTTCATGGTGATTGCGAAACAGCGACGACAGCCCTTGa                           <  3:111688/149‑1 (MQ=255)
                                                                                                                            ggCATCTTCGCCATAATGCTGGAACGGGCGCGAGGCCGTACCGACACGTTCGTCGACATAACCAATGGAGGCCCCGACGATATTGCTGGGGGTGAGCATCAGACGCAGTTCATGGTGATTGCGAAACAGCGACGACAGCCCTTGATaa                        <  1:327769/148‑1 (MQ=255)
                                                                                                                              cATCTTCGCCATAATGCTGGAACGGGCGCGAGGCCGTACCGACACGTTCGTCGACATAACCAATGGAGGCCCCGACGATATTGCTGGGGGTGAGCATCAGACGCAGTTCATGGTGATTGCGAAACAGCGACGACAGCCCTTGATAATTg                     <  3:153449/149‑1 (MQ=255)
                                                                                                                               aTCTTCGCCATAATGCTGGAACGGGCGCGAGGCCGTACCGACACGTTCGTCGACATAACCAATGGAGGCCCCGACGATATTGCTGGGGGTGAGCATCAGACGCAGTTCATGGTGATTGCGAAACAGCGACGACAGCCCTTGATAATTGa                    >  8:249182/1‑149 (MQ=255)
                                                                                                                                tCTTCGCCATAATGCTGGAACGGGCGCGAGGCCGTACCGACACGTTCGTCGACATAACCAATGGAGGCCCCGACGATATTGCTGGGGGTGAGCATCAGACGCAGTTCATGGTGATTGCGAAACAGCGACGACAGCCCTTGATAATTGAc                   >  7:156793/1‑149 (MQ=255)
                                                                                                                                    cGCCATAATGCTGGAACGGGCGCGAGGCCGTACCGACACGTTCGTCGACATAACCAATGGAGGCCCCGACGATATTGCTGGGGGTGAGCATCAGACGCAGTTCATGGTGATTGCGAAACAGCGACGACAGCCCTTGATAATTGAc                   <  8:156793/145‑1 (MQ=255)
                                                                                                                                    cGCCATAATGCTGGAACGGGCGCGAGGCCGTACCGACACGTTCGTCGACATAACCAATGGAGGCCCCGACGATATTGCTGGGGGTGAGCATCAGACGCAGTTCATGGTGATTGCGAAACAGCGACGACAGCCCTTGATAATTGACCaa                <  3:284443/148‑1 (MQ=255)
                                                                                                                                    cGCCATAATGCTGGAACGGGCGCGAGGCCGTACCGACACGTTCGTCGACATAACCAATGGAGGCCCCGACGATATTGCTGGGGGTGAGCATCAGACGCAGTTCATGGTGATTGCGAAACAGCGACGACAGCCCTTGATAATTGACCAAc               >  7:96241/1‑149 (MQ=255)
                                                                                                                                     gCCATAATGCTGGAACGGGCGCGAGGCCGTACCGACACGTTCGTCGACATAACCAATGGAGGCCCCGACGATATTGCTGGGGGTGAGCATCAGACGCAGTTCATGGTGATTGCGAAACAGCGACGACAGCCCTTGATAATTGACCAAcg              >  5:173191/1‑149 (MQ=255)
                                                                                                                                         tAATGCTGGAACGGGCGCGAGGCCGTACCGACACGTTCGTCGACATAACCAATGGAGGCCCCGACGATATTGCTGGGGGTGAGCATCAGACGCAGTTCATGGTGATTGCGAAACAGCGACGACAGCCCTTGATAATTGACCAACGCGAc          <  2:214580/149‑1 (MQ=255)
                                                                                                                                              cTGGAACGGGCGCGAGGCCGTACCGACACGTTCGTCGACATAACCAATGGAGGCCCCGACGATATTGCTGGGGGTGAGCATCAGACGCAGTTCATGGTGATTGCGAAACAGCGACGACAGCCCTTGATAATTGACCAACGCGACCGAGc     <  4:142892/149‑1 (MQ=255)
                                                                                                                                                 gAACGGGCGCGAGGCCGTACCGACACGTTCGTCGACATAACCAATGGAGGCCCCGACGATATTGCTGGGGGTGAGCATCAGACGCAGTTCATGGTGATTGCGAAACAGCGACGACAGCCCTTGATAATTGACCAACGCGACCGAGCCGa  <  5:133801/149‑1 (MQ=255)
                                                                                                                                                 caaCTGGGGCGTGGCCGTACGGACACGTTCGTCGACATAACCAATGGAGGCGCCGACGATATTGCTGGGGGGGAGCATCTGAGGCAGTTCATGGTTATTGCGAAACAGCGACGACAGCCCTTGATAATTGACCAACGCGACCGAGCCGa  <  1:90925/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GAGATTGGGCGTGGTATCGCGGTCATAGCCCAACACGCCGAAATGGTCTGCCCGGGCACTTTCCCCGACCACCAGCACGGTCAGCGACTTGCGCTCATGTTTTAGCCAG‑GCCGCATCACGCTTGGCATCTTCGCCATAATGCTGGAACGGGCGCGAGGCCGTACCGACACGTTCGTCGACATAACCAATGGAGGCCCCGACGATATTGCTGGGGGTGAGCATCAGACGCAGTTCATGGTGATTGCGAAACAGCGACGACAGCCCTTGATAATTGACCAACGCGACCGAGCCGA  >  NC_002947/5211173‑5211465

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