Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,631,143 T→A 12.9% L203Q (CTG→CAG)  lon‑II → DNA‑binding ATP‑dependent protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,631,1430TA12.9% 152.5 / 14.2 62L203Q (CTG→CAG) lon‑IIDNA‑binding ATP‑dependent protease
Reads supporting (aligned to +/- strand):  ref base T (23/31);  new base A (4/4);  total (27/35)
Fisher's exact test for biased strand distribution p-value = 7.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.20e-01

ACGAGCCAGGGCGCCTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATGAAGTCGAAGAGCTGAAAAA  >  NC_002947/2631002‑2631290
                                                                                                                                             |                                                                                                                                                   
acgaGCCAGGGCGCCTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTc                                                                                                                                              <  1:160909/149‑1 (MQ=255)
 cgaGCCAGGGCGCCTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTc                                                                                                                                              >  4:101990/1‑148 (MQ=255)
  gaGCCAGGGCGCCTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCg                                                                                                                                             >  7:213032/1‑148 (MQ=255)
   aGCCAGGGCGCCTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCTGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGAGCTCCTGCAGGTCGa                                                                                                                                            >  3:29560/1‑148 (MQ=255)
   aGCCAGGGCGCCTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGaa                                                                                                                                           <  5:167399/149‑1 (MQ=255)
     ccAGGGCGCCTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGc                                                                                                                                                   >  3:109474/1‑139 (MQ=255)
     ccAGGGCGCCTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGc                                                                                                                                                   <  4:109474/139‑1 (MQ=255)
       aGGGCGCCTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTGCAGGAGGTCGatttcg                                                                                                                                       <  3:45563/149‑6 (MQ=255)
            gcCTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCAtt                                                                                                                                  >  8:13158/1‑149 (MQ=255)
                     acacCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCgg                                                                                                                         <  2:69025/149‑1 (MQ=255)
                     acacCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCgg                                                                                                                         <  2:73537/149‑1 (MQ=255)
                       acCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGt                                                                                                                       <  7:7258/149‑1 (MQ=255)
                        ccATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGt                                                                                                                       >  7:292950/1‑148 (MQ=255)
                        ccATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATACCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTc                                                                                                                      <  8:222545/149‑1 (MQ=255)
                          aTGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTc                                                                                                                      >  7:293591/1‑147 (MQ=255)
                             gCCGCGCACATTGCCCTGAAGATCGAGCAGAAGCATGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATTCCGAAATCTACCTCCTGCAGGTCTAAAAGCGCATTCGCGGCCGGGTCaaga                                                                                                                  <  2:160875/148‑1 (MQ=255)
                             gCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGTGTCaaga                                                                                                                  <  4:213796/148‑1 (MQ=255)
                                      aTGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGa                                                                                                        <  7:82412/149‑1 (MQ=255)
                                              gAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATg                                                                                                          >  6:326715/1‑139 (MQ=255)
                                              gAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATg                                                                                                          <  5:326715/139‑1 (MQ=255)
                                                     gAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGa                                                                                         <  3:101990/149‑1 (MQ=255)
                                                     gAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGa                                                                                         <  5:128738/149‑1 (MQ=255)
                                                      aGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGc                                                                                                               <  6:209852/126‑1 (MQ=255)
                                                      aGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGc                                                                                                               >  5:209852/1‑126 (MQ=255)
                                                          gAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGa                                                                                                        <  4:157280/129‑1 (MQ=255)
                                                          gAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGa                                                                                                        >  3:157280/1‑129 (MQ=255)
                                                            aGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGtactac                                                                                  <  8:93576/149‑1 (MQ=255)
                                                                 gcaaTTCTTGAATTCCTCGAGTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGATGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCCCATTCGCGGCCGGGTCAAGAACCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGaa                                                                             <  4:29560/147‑1 (MQ=255)
                                                                 gAAATTCTTGAAATTGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGaa                                                                             >  4:143194/1‑149 (MQ=255)
                                                                       cTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGc                                                                        <  7:286026/148‑1 (MQ=255)
                                                                        ttGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAg                                                                         >  4:255343/1‑146 (MQ=255)
                                                                                 tcgACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGcc                                                             <  8:293591/149‑1 (MQ=255)
                                                                                    aCTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCaagaa                                                                                                                 <  3:188047/94‑1 (MQ=255)
                                                                                    aCTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCaagaa                                                                                                                 >  4:188047/1‑94 (MQ=255)
                                                                                        gACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGc                                                                        >  6:133345/1‑131 (MQ=255)
                                                                                        gACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGc                                                                        <  5:133345/131‑1 (MQ=255)
                                                                                        gACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCa                                                        <  6:54789/147‑1 (MQ=255)
                                                                                         aCTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCaga                                                      >  1:280908/1‑148 (MQ=255)
                                                                                              ccGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATTGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATg                                                                                                          >  7:130984/1‑91 (MQ=255)
                                                                                              ccGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATTGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATg                                                                                                          <  8:130984/91‑1 (MQ=255)
                                                                                              ccGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAAc                                                >  4:86736/1‑149 (MQ=255)
                                                                                                 tggcGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTg                                             <  3:169514/146‑1 (MQ=255)
                                                                                                     gAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGa                                         >  1:264732/1‑149 (MQ=255)
                                                                                                       gCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATg                                       <  2:264732/149‑1 (MQ=255)
                                                                                                           gTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGa                                   <  8:213032/149‑1 (MQ=255)
                                                                                                                 gcgcTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCagaaag                                                   >  4:289121/1‑127 (MQ=255)
                                                                                                                 gcgcTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCagaaag                                                   <  3:289121/127‑1 (MQ=255)
                                                                                                                   cctgcagGAGGTCGATTTCGGCATCCAGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGc                            >  5:21953/2‑148 (MQ=255)
                                                                                                                       ctgGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCcaca                        <  2:280908/148‑1 (MQ=255)
                                                                                                                       ctgGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGACCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACaa                       >  4:140231/1‑149 (MQ=255)
                                                                                                                        tgGATGCCGAAATCGTCCTCCTTCAGGTCGAAAAGCGCATTCGCGGCCGTGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTTAATGAGCAGATGAAGGCCATTCAGAAAGAACtt                                              <  8:295474/125‑1 (MQ=255)
                                                                                                                        tgGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACtt                                              >  7:295474/1‑125 (MQ=255)
                                                                                                                          gATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGAATTCGCGGCCGGGTCATGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACATCAATGAGCAGATGAAGGCAATTCAGATAGAACTTGGCGATGGTGACGAAGGCCACAATg                     >  6:319969/1‑148 (MQ=255)
                                                                                                                               gcttttCGACCTGCAGGAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATGAAGTc                >  2:89926/6‑148 (MQ=255)
                                                                                                                               gcttttCGACCTGCAGGAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATGAAGTc                <  1:89926/143‑1 (MQ=255)
                                                                                                                               cGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAAt                      <  6:313875/142‑1 (MQ=255)
                                                                                                                               cGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAAt                      >  5:313875/1‑142 (MQ=255)
                                                                                                                                 ttttCGACCTGCAGGAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATGAAGTCg               >  8:152044/4‑147 (MQ=255)
                                                                                                                                 ttttCGACCTGCAGGAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATGAAGTCg               >  8:156550/4‑147 (MQ=255)
                                                                                                                                  aaTCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATGAAGTCGAag            <  8:6439/149‑1 (MQ=255)
                                                                                                                                   ttCGACCTGCAGGAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATGAAGTCGAag            <  7:156550/147‑1 (MQ=255)
                                                                                                                                   ttCGACCTGCAGGAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATGAAGTCGAag            <  7:152044/147‑1 (MQ=255)
                                                                                                                                            ctGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATGAAGTCGAAGAGCTGaaaaa  <  3:143194/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
ACGAGCCAGGGCGCCTGGTCGACACCATGGCCGCGCACATGGCCCTGAAGATCGAGCAGAAGCAGGAAATTCTTGAAATCGTCGACTTGACTACCCGTGTCGAGCACGTGCTGGCGCTGCTGGATGCCGAAATCGACCTCCTGCAGGTCGAAAAGCGCATTCGCGGCCGGGTCAAGAAGCAGATGGAGCGCAGCCAGCGCGAGTACTACCTGAATGAGCAGATGAAGGCCATTCAGAAAGAACTTGGCGATGGTGACGAAGGCCACAATGAAGTCGAAGAGCTGAAAAA  >  NC_002947/2631002‑2631290

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: