Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,286,952 C→A 7.1% L118L (CTG→CTT)  ggt ← gamma‑glutamyltranspeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,286,9520CA 7.1% 211.2 / 4.9 70L118L (CTG→CTT) ggtgamma‑glutamyltranspeptidase
Reads supporting (aligned to +/- strand):  ref base C (29/36);  new base A (2/3);  total (31/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01

GTACTGCTTTTCCGCCACCTTGAAGCCGTTTTTCGCATAGCCGATGGCCGGGGTCAGCAACTCGCTCCAGGGCAGCTTGCCGAACTTCTGGTGCGCCTCCCACAGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGTTCTCGATGACCTCGCCCTTGTCGTCCAGGTACATGTTGCGCGTGGCGGCCTTGGGCGCGACTTCGCGGTAGTCGAGGAAGTAAGGCTTGCCATCGACGAACAGGGTCATGAAGCCGCCACCGCCGATATTACCGGCCT  >  NC_002947/5286805‑5287100
                                                                                                                                                   |                                                                                                                                                    
gtactgCTTTTCCGCCACCTTGAAGCCGTTTTTCGCATAGCCGATGGCCGGGGTCAGCAACTCGATCCAGGGCAGCTTGCCGAACTTCTGGTGCGCCTCCCACAGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCa                                                                                                                                                     >  7:192848/1‑149 (MQ=255)
  actgCTTTTCCGCCACCTTGAAGCCGTTTTTCGCATAGCCGATGGCCGGGGTCAGCAACTCGCTCCAGGGCAGCTTGCCGAACTTCTGGTGCGCCTCCCACAGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCAgg                                                                                                                                                   >  5:177133/1‑149 (MQ=255)
      cTTTTCCGCCACCTTGAAGCCGTTTTTCGCATAGCCGATGGCCGGGGTCAGCAACTCGCTCCAGGGCAGCTTGCCGAACTTCTGGTGCGCCTCCCACAGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCAGGCt                                                                                                                                                 >  7:83673/1‑147 (MQ=255)
      cTTTTCCGCCACCTTGAAGCCGTTTTTCGCATAGCCGATGGCCGGGGTCAGCAACTCGCTCCAGGGCAGCTTGCCGAACTTCTGGTGCGCCTCCCACAGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTg                                                                                                                                                <  3:12483/148‑1 (MQ=255)
       ttttCCGCCACCTTGAAGCCGTTTTTCGCATAGCCGATGGCCGGGGTCAGCAACTCGCTCCAGGGCAGCTTGCCGAACTTCTGGTGCGCCTCCCACAGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACAAgggcgag                                                                                                                                              >  4:252168/1‑144 (MQ=255)
         ttCCGCCACCTTGAAGCCGTTTTTCGCATAGCCGATGGCCGGGGTCAGCAACTCGCTCCAGGGCAGCTTGCCGAACTTCTGGTGCGCCTCCCACAGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAg                                                                                                                                              <  3:246946/147‑1 (MQ=255)
              ccACCTTGAAGCCGTTTTTCGCATAGCCGATGGCCGGGGTCAGCAACTCGCTCCAGGGCAGCTTGCCGAACTTCTGGTGCGCCTCCCACAGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGTtctc                                                                                                                                        >  8:285279/1‑148 (MQ=255)
              ccACCTTGAAGCCGTTTTTCGCATAGCCGATGGCCGGGGTCAGCAACTCGCTCCAGGGCAGCTTGCCGAACTTCTGGTGCGCCTCCCACAGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACAAGGGCGAGGtcatc                                                                                                                                        >  5:119740/1‑144 (MQ=255)
                  cTTGAAGCCGTTTTTCGCATAGCCGATGGCCGGGGTCAGCAACTCGCTCCAGGGCAGCTTGCCGAACTTCTGGTGCGCCTCCCACAGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGTTCTCGATGa                                                                                                                                   >  5:42287/1‑149 (MQ=255)
                   ttGAAGCCGTTTTTCGCATAGCCGATGGCCGGGGTCAGCAACTCGCTCCAGGGCAGCTTGCCGAACTTCTGGTGCGCCTCACACAGGCCCATTACCGTCCCGGGTACGACTGCGGCGCGTGCCCCGACCACGCTGAGGTTCTCGATGAc                                                                                                                                  >  6:118042/1‑149 (MQ=255)
                       aGCCGTTTTTCGCATAGCCGATGGCCGGGGTCAGCAACTCGCTCCAGGGCAGCTTGCCGAACTTCTGGTGCGCCTCCCACAGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGTTCTCGATGACCTCg                                                                                                                              >  4:62965/1‑149 (MQ=255)
                        gCCGTTTTTCGCATAGCCGATGGCCGGGGTCAGCAACTCGCTCCAGGGCAGCTTGCCGAACTTCTGGTGCGCCTCCCACAGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGTTCTCGATGACCTCGc                                                                                                                             >  6:224087/1‑149 (MQ=255)
                         ccGTTTTTCGCATAGCCGATGGCCGGGGTCAGCAACTCGCTCCAGGGCAGCTTGCCGAACTTCTGGTGCGCCTCCCACAGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGTTCTCGATGACCTCGcc                                                                                                                            >  8:23238/1‑149 (MQ=255)
                          cGTTTTTCGCATAGCCGATGGCCGGGGTCAGCAACTCGCTCCAGGGCAGCTTGCCGAACTTCTGGTGCGCCTCCCACAGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGTTCTCGATGACCTCGccc                                                                                                                           >  1:292194/1‑149 (MQ=255)
                              tttCGCAAAGCCGATGGCCGGGGTCAGCAACTCGCTCCAGGGCAGCTTGCCGAACTTCTGGTGCGCCTCCCACAGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGTTCTCGATGACCTCGCCCTTg                                                                                                                        <  7:131894/148‑1 (MQ=255)
                               cgagcttAGCCGATGGCCGGGGTCAGCAACTCGCTCCAGGGCAGCGTGCCGAACTTCTGGTGCGCCTCCCACAGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGTTCTCGATGACCTCGCCCTTgtc                                                                                                                      <  3:78138/143‑1 (MQ=255)
                                  gCATAGCCGATGGCCGGGGTCAGCAACTCGCTCCAGGGCAGCTTGCCGAACTTCTGGTGCGCCTCCCACAGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGTTCTCGATGACCTCGCCCTTgtcgtc                                                                                                                   <  5:313758/149‑1 (MQ=255)
                                  gCATAGCCGATGGCCGGGGTCAGCAACTCGCTCCAGGGCAGCTTGCCGAACTTCTGGTGCGCCTCCCACAGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACAAGGGCGAGGTCATCGAGAACCTCAGCCTGgtcggg                                                                                                                   <  5:315674/149‑3 (MQ=37)
                                  gCATAGCCGATGGCCGGGGTCAGCAACTCGCTCCAGGGCAGCTTGCCGAACTTCTGGTGCGCCTCCCACAGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACAAGGGCGAGGTCATCGAGAACCTCAGCCTGgtcggg                                                                                                                   <  3:252168/149‑3 (MQ=37)
                                  gCATAGCCGATGGCCGGGGTCAGCAACTCGCTCCAGGGCAGCTTGCCGAACTTCTGGTGCGCCTCCCACAGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACAAGGGCGAGGTCATCGAGAACCTCAGCCTGgtcgg                                                                                                                    <  6:119740/148‑2 (MQ=37)
                                  gCATAGCCGATGGCCGGGGTCAGCAACTCGCTCCAGGGCAGCTTGCCGAACTTCTGGTGCGCCTCCCACAGGCCCAATACCGTACCCCGTACCCCTGCGGCGCGTGCCCCGACCAGGCTGAGGATCTCGATGACCTCCCCCTTgtcgtc                                                                                                                   >  1:197607/1‑149 (MQ=255)
                                   cATAGCCGATGGCCGGGGTCAGCAACTCGCTCCAGGGCAGCTTGCCGAACTTCTGGTGCGCCTCCCACAGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGTTCTCGATGACCTCGCCCTTgtcgtc                                                                                                                   >  3:181946/1‑148 (MQ=255)
                                     tAGCCGATGGCCGGGGTCAGCAACTCGCTCCAGGGCAGCTTGCCGAACTTCTGGTGCGCCTCCCACAGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGTTCTCGATGACCTCGCCCTTGTCGTCCAg                                                                                                                <  8:83673/149‑1 (MQ=255)
                                      aGCCGATGGCCGGGGTCAGCAACTCGCTCCAGGGCAGCTTGCCGAACTTCTGGTGCGCCTCCCACAGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCAGGCGGAGGGTCTCGATGACCTCGCCCCTgtc                                                                                                                      >  4:78138/1‑142 (MQ=255)
                                          gATGGCCGGGGTCAGCAACTCGCTCCAGGGCAGCTTGCCGAACTTCTGGTGCGCCTCCCACAGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGTTCTCGATGACCTCGCCCTTg                                                                                                                        <  2:46830/136‑1 (MQ=255)
                                            tGGCCGGGGTCAGCAACTCGCTCCAGGGCAGCTTGCCGAACTTCTGGTGCGCCTCCCACAGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGTTCTCGATGACCTCGCCCTTGTCGTCCAGGTACATg                                                                                                         <  5:394106/149‑1 (MQ=255)
                                               ccGGGGTCAGCAACTCGCTCCAGGGCAGCTTGCCGAACTTCTGGTGCGCCTCCCACAGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGTTCTCGATGACCTCGCCCTTGTCGTCCAGGTACATGTTg                                                                                                      >  8:31369/1‑149 (MQ=255)
                                               ccGGGGTCAGCAACTCGCTCCAGGGCAGCTTGCCGAACTTCTGGTGCGCCTCCCACAGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGTTCTCGATGACCTCGCCCTTGTCGTCCAGGTACATGTTg                                                                                                      <  7:285279/149‑1 (MQ=255)
                                                cgtgttCAGCAACTCGCTCCAGGGCAGCTTGCCGAACTTCTGGTGCGCCTCCCACAGGCACATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGTTCTCGATGACCTCGCCCTTGTCGTCCAGGTACATGTTgc                                                                                                     <  8:192848/144‑1 (MQ=255)
                                                       aGCAACTCGCTCCAGGGCAGCTTGCCGAACTTCTGGTGCGCCTCCCACAGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGTTCTCGATGACCTCGCCCTTGTCGTCCAGGTACATGTTGCGCGTggc                                                                                              >  1:134007/1‑149 (MQ=255)
                                                             tCGCTCCAGGGCAGCTTGCCGAACTTCTGGTGCGCCTCCCACAGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGTTCTCGATGACCTCGCCCTTGTCGTCCAGGTACATGTTGCGCGTGGCGGCCtt                                                                                        >  1:287540/1‑149 (MQ=255)
                                                               gCTCCAGGGCAGCTTGCCGAACTTCTGGTGCGCCTCCCACAGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGTTCTCGATGACCTCGCCCTTGTCGTCCAGGTACATGTTGCGCGTGGCGGCCtt                                                                                        <  3:200588/147‑1 (MQ=255)
                                                               gCTCCAGGGCAGCTTGCCGAACTTCTGGTGCGCCTCCCACAGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGTTCTCGATGACCTCGCCCTTGTCGTCCAGGTACATGTTGCGCGTGGCGGCCtt                                                                                        <  3:263097/147‑1 (MQ=255)
                                                               gCTCCAGGGCAGCTTGCCGAACTTCTGGTGCGCCTCCCACAGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGTTCTCGATGACCTCGCCCTTGTCGTCCAGGTACATGTTGCGCGTGGCGGCCTTgg                                                                                      <  7:23238/149‑1 (MQ=255)
                                                               gCTCCAGGGCAGCTTGCCGAACTTCTGGTGCGCCTCCCACAGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGTTCTCGATGACCTCGCCCTTGTCGTCCAGGTACATGTTGCGCGTGGCGGCCTTgg                                                                                      <  2:283944/149‑1 (MQ=255)
                                                                   cAGGGCAGCTTGCCGAACTTCTGGTGCGCCTCCCACAGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGTTCTCGATGACCTCGCCCTTGTCGTCCAGGTACATGTTGCGCGTGGCGGCCTTGGgc                                                                                    <  7:31369/147‑1 (MQ=255)
                                                                    aGGGCAGCTTGCCGAACTTCTGGTGCGCCTCCCACAGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGTTCTCGATGACCTCGCCCTTGTCGTCCAGGTACATGTTGCGCGTGGCGGCCTTGGgcgc                                                                                  >  6:87754/1‑148 (MQ=255)
                                                                                cGAACTTCTGGTGCGCCTCCCACAGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGTTCTCGATGACCTCGCCCTTGTCGTCCAGGTACATGTTGCGCGTGGCGGCCTTGGGCGCGACTTCGCGGTAg                                                                     >  5:165055/1‑149 (MQ=255)
                                                                                  aaCTTCTGGTGCGCCTCCCACAGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGTTCTCg                                                                                                                                       >  8:83116/1‑81 (MQ=255)
                                                                                  aaCTTCTGGTGCGCCTCCCACAGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGTTCTCg                                                                                                                                       <  7:83116/81‑1 (MQ=255)
                                                                                  aaCTTCTGGTGCGCCTCCCACAGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGTTCTCGATGACCTCGCCCTTGTCGTCCAGGTACATGTTGCGCGTGGCGGCCTTGGGCGCGACTTCGCGGTAGTc                                                                   >  3:377674/1‑149 (MQ=255)
                                                                                       cTGGTGCGCCTCCCACAGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGTTCTCGATGACCTCGCCCTTGTCGTCCAGGTACATGTTGCGCGTGGCGGCCTTGGGCGCGACTTCGCGGTAGTCGAgg                                                               <  6:165055/148‑1 (MQ=255)
                                                                                          gTGCGCCTCCCACAGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGTTCTCGATGACCTCGCCCTTGTCGTCCAGGTACATGTTGCGCGTGGCGGCCTTGGGCGCGACTTCGCGGTAGTCGAGGaa                                                             <  3:62965/147‑1 (MQ=255)
                                                                                            gcgcCTCCCCCAGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGGTCTCGATGACCTCGCCCTGGTCGTCCAGGTACCTGTTGAGCGTGGCGGCCTTGGGCGCGACATCGCTGTAGTCGAGGAAGt                                                           >  2:316527/1‑147 (MQ=255)
                                                                                            gcgcCTCCCACAGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGTTCTCGATGACCTCGCCCTTGTCGTCCAGGTACATGTTGCGCGTGGCGGCCTTGGGCGCGACTTCGCGGTAGTCGAGGAAGTaa                                                         <  3:298325/149‑1 (MQ=255)
                                                                                              gcCTCCCACAGGCCCATTACCGTACCGTGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGTTCTCGATGACCTCGCCCTTGTCGTCCAGGTACATGTTGCGCGTGGCGGCCTTGGGCGCGACTTCGCGGTAGTCGAGGAAGTAAgg                                                       <  5:180430/149‑1 (MQ=255)
                                                                                              gcCTCCCACAGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGTTCTCGATGACCTCGCCCTTGTCGTCCAGGTACATGTTGCGCGTGGCGGCCTTGGGCGCGACTTCGCGGTAGTCGAGGAAGTAAgg                                                       <  5:224087/149‑1 (MQ=255)
                                                                                               ccTCCCACAGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGTTCTCGATGACCTCGCCCTTGTCGTCCAGGTACATGTTGCGCGTGGCGGCCTTGGGCGCGACTTCGCGGTAGTCGAGGAAGTAAGGc                                                      <  2:30918/149‑1 (MQ=255)
                                                                                                  cccACAGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGTTCTCGATGACCTCGCCCTTGTCGTCCAGGTACATGTTGCGCGTGGCGGCCTTGGGCGCGACTTCGCGGTAGTCGAGGAAGTAAGGCTTg                                                   <  2:292194/149‑1 (MQ=255)
                                                                                                   ccAGAGGCACATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGGTCTCGATGACCTCGCCCTTGTCGTCCAGGTACATGTTGCGCGTGGCGGCGTTGGGCGCGACTTCGCGGTAGTCGAGGAAGTAAGGCTTg                                                   <  5:118042/148‑1 (MQ=255)
                                                                                                   ccACAGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGTTCTCGATGACCTCGCCCTTGTCGTCCAGGTACATGTTGCGCGTGGCGGCCTTGGGCGCGACTTCGCGGTAGTCGAGGAAGTAAGGCTTg                                                   <  4:201589/148‑1 (MQ=255)
                                                                                                   ccACAGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGTTCTCGATGACCTCGCCCTTGTCGTCCAGGTACATGTTGCGCGTGGCGGCCTTGGGCGCGACTTCGCGGTAGTCGAGGAAGTAAGGCTTGc                                                  >  3:65269/1‑149 (MQ=255)
                                                                                                    cacaGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGTTCTCGATGACCTCGCCCTTGTCGTCCAGGTACATGTTGCGCGTGTCGGCCTTGGGCGCGACTTCGCGGTAGTCGAGGAAGTAAGGCTTGcc                                                 <  2:197607/149‑1 (MQ=255)
                                                                                                     acaGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGTTCTCGATGACCTCGCCCTTGTCGTCCAGGTACATGTTGCGAGTGGCGGCCTTGGGCGCGACTTCGCGGTAGTCGAGGAAGTAAGGCTTGc                                                  >  8:254538/1‑147 (MQ=255)
                                                                                                      caGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGTTCTCGATGACCTCGCCCTTGTCGTCCAGGTACATGTTGCGCGTGGCGGCCTTGGGCGCGACTTCGCGGTAGTCGAGGAAGTAAGGCTTGCCAt                                               <  6:42287/149‑1 (MQ=255)
                                                                                                         gaccATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGTTCTCGATGACCTCGCCATTGTCGTCCAGGTACATGTTGCGCGTGGCTGCCTTGGGCGCGACTTCGCGTTAGTCGAGGAAGTAAGGCTTGCCATcga                                            <  1:391416/147‑1 (MQ=255)
                                                                                                           ccATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGTTCTCGATGACCTCGCCCTTGTCGTCCAGGTACATGTTGCGCGTGGCGGCCTTGGGCGCGACTTCGCGGTAGTCGAGGAAGTAAGGCTTGCCATcgac                                           >  3:162010/1‑148 (MQ=255)
                                                                                                              ttACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGTTCTCGATGACCTCGCCCTTGTCGTCCAGGTACATGTTGCGCGTGGCGGCCTTGGGCGCGACTTCGCGGTAGTCGAGGAAGTAAGGCTTGCCATCGACGAAc                                       <  5:87754/149‑1 (MQ=255)
                                                                                                                 ccgtaccgGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGTTCTCGATGACCTCGCCCTTGTCGTCCAGGTACATGTTGCGCGTGGCGGCCTTGGGCGCGACTTCGCGGTAGTCGAGGAAGTAAGGCTTGCCATCGACGAACAgg                                    <  6:177133/149‑1 (MQ=255)
                                                                                                                 ccgtaccgGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGTTCTCGATGACCTCGCCCTTGTCGTCCAGGTACATGTTGCGCGTGGCGGCCTTGGGCGCGACTTCGCGGTAGTCGAGGAAGTAAGGCTTGCCATCGACGAACAgg                                    <  4:181946/149‑1 (MQ=255)
                                                                                                                    taccgGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGTTCTCGAGGACCTCGCCGTTGTCGTCCAGGTCCATGTTGCGCGTGGCGGCGTGGGGCGCGACTTCGCGGTAGTCGAGGAGGTAAGGCTTGCCAGCGACGGACAGGGTc                                 >  4:211693/1‑149 (MQ=255)
                                                                                                                       cgGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGTTCTCGATGACCTCGCCCTTGTCGTCCAGGTACATGTTGCGCGTGGCGTCCTTGGGCGCGACTTCGCGGTAGTCGAGGATGTAAGGCTTGCCATCGACGAACAGGGTCATg                              <  1:349980/149‑1 (MQ=255)
                                                                                                                         ggTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGTTCTCGATGACCTCGCCCTTGTCGTCCAGGTACATGTTGCGCGTGGCGGCCTTGGGCGCGACTTCGCGGTAGTCGAGGAAGTAAGGCTTGCCATCGACGAACAGGGTCATGaa                            <  2:287540/149‑1 (MQ=255)
                                                                                                                                cTGCGGCGCGTGCCCCGACCAGGCTGAGGTTCTCGATGACCTCGCCCTTGTCGTCCAGGTACATGTTGCGCGTGGCGGCCTTGGGCGCGACTTCGCGGTAGTCGAGGAAGTAAGGCTTGCCATCGACGAACAGGGTCATGAAGCCGCCa                     >  4:354983/1‑149 (MQ=255)
                                                                                                                                       gcgTGCCCCGACCAGGCTGAGGTTCTCGATGACCTCGCCCTTGTCGTCCAGGTACATGTTGCGCGGGGCGGCCTTGGGCGCGACTTCGCGGTAGTCGAGGAAGTAAGGCTTGCCATCGACGAACAGGGTCATGAAGCCGCCACCGCCGa              <  3:354983/149‑1 (MQ=255)
                                                                                                                                          tGCCCCGACCAGGCTGAGGTTCTCGATGACCTCGCCCTTGTCGTCCAGGTACATGTTGCGCGTGGCGGCCTTGGGCGCGACTTCGCGGTAGTCGAGGAAGTAAGGCTTGCCATCGACGAACAGGGTCATGAAGCCGCCACCGCCGatat           >  7:220107/1‑149 (MQ=255)
                                                                                                                                            ccccGACCAGGCTGAGGTTCTCGATGACCTCGCCCTTGTCGTCCAGGTACATGTTGCGCGTGGCGGCCTTggg                                                                                     <  6:73586/73‑1 (MQ=255)
                                                                                                                                            ccccGACCAGGCTGAGGTTCTCGATGACCTCGCCCTTGTCGTCCAGGTACATGTTGCGCGTGGCGGCCTTggg                                                                                     >  5:73586/1‑73 (MQ=255)
                                                                                                                                            ccccGACCAGGCTGAGGTTCTCGATGACCTCGCCCTTGTCGTCCAGGTACATGTTGCGCGTGGCGGCCTTGGGCGCGACTTCGCGGTAGTCGAGGAAGTAAGGCTTGCCATCGACGAACAGGGTCATGAAGCCGCCACCGCCGATATTa         >  6:52293/1‑149 (MQ=255)
                                                                                                                                               cGACCAGGCTGAGGTTCTCGATGACCTCGCCCTTGTCGTCCAGGTACATGTTGCGCGTGGCGGCCTTGGGCGCGACTTCGCGGTAGTCGAGGAAGTAAGGCTTGCCATCGACGAACAGGGTCATGAAGCCGCCACCGCCGATATTACCg      <  4:162010/149‑1 (MQ=255)
                                                                                                                                                   cAGGCTGGGGTTCTCGATGACCTCGCCCTTGTCGTCCAGGTACATGTTGCGCGTGGCGGCCTTGGGCGCGACTTCGCGGTATTCGAGGAAGTAAGGCTTGCCATCGACGAACAGGGTCATGAAGCCGCCACCGCCGATATTACCGGCct  >  2:128355/1‑149 (MQ=255)
                                                                                                                                                   cAGGCTGAGGTTCTCGATGACCTCGCCCTTGTCGTCCAGGTACATGTTGCGCGTGGCGGCCTTGGGCGCGACTTCGCGGTAGTCGAGGAAGTAAGGCTTGCCATCGACGAACAGGGTCATGAAGCCGCCACCGCCGATATTACCGGCct  >  6:195429/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GTACTGCTTTTCCGCCACCTTGAAGCCGTTTTTCGCATAGCCGATGGCCGGGGTCAGCAACTCGCTCCAGGGCAGCTTGCCGAACTTCTGGTGCGCCTCCCACAGGCCCATTACCGTACCGGGTACGCCTGCGGCGCGTGCCCCGACCAGGCTGAGGTTCTCGATGACCTCGCCCTTGTCGTCCAGGTACATGTTGCGCGTGGCGGCCTTGGGCGCGACTTCGCGGTAGTCGAGGAAGTAAGGCTTGCCATCGACGAACAGGGTCATGAAGCCGCCACCGCCGATATTACCGGCCT  >  NC_002947/5286805‑5287100

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