Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,962,391 A→G 5.5% L13P (CTG→CCG)  recC ← exonuclease V (RecBCD complex), gamma chain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,962,3910AG 5.5% 229.1 / 3.0 73L13P (CTG→CCG) recCexonuclease V (RecBCD complex), gamma chain
Reads supporting (aligned to +/- strand):  ref base A (29/40);  new base G (2/2);  total (31/42)
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.25e-01

AAATCAATGTTTGCCGCAATACCAAACTTTTGCGACAGGGTCATTTGCAGCCACTGTGCCATACCGGTACTTTGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGA  >  NC_000913/2962244‑2962536
                                                                                                                                                   |                                                                                                                                                 
aaaTCAATGTTTGCCGCAATACCAAACTTTTGCGACAGGGTCATTTGCAGCCACTGTGCCATACCGGTACTTTGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAg                                                                                                                                                  <  4:81853/149‑1 (MQ=255)
aaaTCAATGTTTGCCGCAATACCAAACTTTTGCGACAGGGTCATTTGCAGCCACTGTGCCATACCGGTACTTTGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAg                                                                                                                                                  <  6:302813/149‑1 (MQ=255)
aaaTCAATGTTTGCCGCAATACCAAACTTTTGCGACAGGGTCATTTGCAGCCACTGTGCCATACCGGTACTTTGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAg                                                                                                                                                  <  5:202780/149‑1 (MQ=255)
 aaTCAATGTTTTCCGCAATACCACACTTTTGCGACAGGGTCATTTGCAGCCAGTGTGCCATACCGGTACTTTGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGc                                                                                                                                                 <  1:98042/149‑1 (MQ=255)
 aaTCAATGTTTGCCGGAATTCCAAACTTTTGCGTCAGGGTCATTTGCAGCGACTGTGCCATACCGGTAATTTGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGc                                                                                                                                                 <  8:169057/149‑1 (MQ=255)
     aaTGTTTGCCGCAATACCAAACTTTTGCGACAGGGTCATTTGCAGCCACTGTGCCATACCGGTACTTTGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACg                                                                                                                                              <  1:162049/148‑1 (MQ=255)
         tttGCCGCAATACCAAACTTTTGCGACAGGGTCATTTGCAGCCACTGTGCCATACCGGTACTTTGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTc                                                                                                                                            >  5:318312/1‑146 (MQ=255)
             ccGCTATAGCATTCTTTTGCGTCAGGGTCTTTTGCAGACACTGTGCCATACCGGTACTTTGCACCAGAATCAGCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCTCGGCCAGACGa                                                                                                                                     <  6:172690/149‑1 (MQ=255)
             ccGCAATACCAAACTTTTGCGACAGGGTCATTTGCAGCCACTGTGCCATACCGGTACTTTGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGa                                                                                                                                     >  8:298845/1‑149 (MQ=255)
                cAATACCAAACTTTTGCGACAGGGTCATTTGCAGCCACTGTGCCATACCGGTACTTTGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTg                                                                                                                                  <  7:142715/149‑1 (MQ=255)
                cAATACCAAACTTTTGCGACAGGGTCATTTGCAGCCACTGTGCCATACCGGTACTTTGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTg                                                                                                                                  >  8:339137/1‑149 (MQ=255)
                   tACCAAACTTTTGCGACAGGGTCATTTGCAGCCACTGTGCCATACCGGTACTTTGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGaa                                                                                                                               <  2:229348/149‑1 (MQ=255)
                          cTTTTGCGACAGGGTCATTTGCAGCCACTGTGCCATACCGGTACTTTGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGa                                                                                                                        <  4:10786/149‑1 (MQ=255)
                           ttttGCGACAGGGTCATTTGCAGCCACTGTGCCATACCGGTACTTTGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGAc                                                                                                                       >  3:39212/1‑149 (MQ=255)
                             ttGCGACAGGGTCATTTGCAGCCACTGTGCCATACCGGTACTTTGCACCAGAATCAGCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATACACGCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGcccc                                                                                                                     >  7:358713/1‑149 (MQ=255)
                              tGCGACAGGGTCATTTGCAGCCACTGTGCCATACCGGTACTTTGCACCAGAATCATCTCTGGTTCGAAAGGAGCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGAccc                                                                                                                     <  4:145761/148‑1 (MQ=255)
                                      ggTCATTTGCAGCCACTGTGCCATACCGGTACTTTGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGCCAATACACTCCAGCAACGCTTCGAGCACGTCCAGCCGATTGGAATGGTAGACCCTTAACATAg                                                                                                            >  6:169371/1‑149 (MQ=255)
                                       gTCATTTGCAGCCACTGTGCCATACCGGTACTTTGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGc                                                                                                           <  6:318312/149‑1 (MQ=255)
                                           tttGCAGCCACTGTGCCATACCGGTACTTTGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGAGTGGAATGGTAGACCCTTAACATAGCGGCt                                                                                                       >  2:30551/1‑149 (MQ=255)
                                           tttGCAGCCACTGTGCCATACCGGTACTATGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGAAGACCCTAAAATAGCGGcac                                                                                                      >  6:8160/1‑147 (MQ=255)
                                                 gCCACTGTGCCATACCGGTACTTTGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGa                                                                                                  >  4:240705/1‑148 (MQ=255)
                                                  ccACTGTGCCATACCGGTACTTTGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGac                                                                                                 >  1:51557/1‑148 (MQ=255)
                                                    aCTGTGCCATACCCGTACTTTGCACAAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGactac                                                                                              <  6:35004/149‑1 (MQ=255)
                                                        tgCCATACCGGTACTTTGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGAc                                                                                          >  3:116130/1‑149 (MQ=255)
                                                        tgCCATACCGGTACTTTGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGAc                                                                                          <  7:48301/149‑1 (MQ=255)
                                                             tACCGGTACTTTGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCGGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGa                                                                                           >  1:110509/1‑143 (MQ=255)
                                                             tACCGGTACTTTGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCGGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGa                                                                                           <  2:110509/143‑1 (MQ=255)
                                                               ccGGTACTTTGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGAtt                                                                                                                                   >  6:284131/1‑101 (MQ=255)
                                                               ccGGTACTTTGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGAtt                                                                                                                                   <  5:284131/101‑1 (MQ=255)
                                                               ccGGTACTTTGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGAc                                                                                                                       >  1:43573/1‑113 (MQ=255)
                                                               ccGGTACTTTGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGAc                                                                                                                       <  2:43573/113‑1 (MQ=255)
                                                                cGGTACTTTGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCggg                                                                                  <  4:116130/149‑1 (MQ=255)
                                                                cGGTACTTTGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCggg                                                                                  <  5:254483/149‑1 (MQ=255)
                                                                   tACTTTGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCaa                                                                               <  3:201315/149‑1 (MQ=255)
                                                                     cTTTGCACCAGAATCATCTCTGGTTCGAACGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATg                                                                             <  3:65939/149‑1 (MQ=255)
                                                                            ccAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGc                                                                        >  1:51383/1‑147 (MQ=255)
                                                                            ccAGAATCATCTCTGGGTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCgg                                                                      <  6:288204/149‑1 (MQ=255)
                                                                             cAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCggg                                                                     >  1:297501/1‑149 (MQ=255)
                                                                                  tcatcTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCAtc                                                                >  8:125441/1‑149 (MQ=255)
                                                                                    atcTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCAtctc                                                              <  2:161152/149‑1 (MQ=255)
                                                                                     tctcTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGc                                                                                 >  4:320664/1‑129 (MQ=255)
                                                                                     tctcTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGc                                                                                 <  3:320664/129‑1 (MQ=255)
                                                                                     tctcTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCg                                                                       <  7:168545/139‑1 (MQ=255)
                                                                                     tctcTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCg                                                                       >  8:168545/1‑139 (MQ=255)
                                                                                       tctGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCgg                                                                      <  5:301842/138‑1 (MQ=255)
                                                                                       tctGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCgg                                                                      >  6:301842/1‑138 (MQ=255)
                                                                                              cGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTctct                                                     <  7:118526/148‑1 (MQ=255)
                                                                                                aaaGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGcc                                                  <  1:262564/149‑1 (MQ=255)
                                                                                                  aGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGAACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTATCTTcccc                                                >  8:11518/1‑149 (MQ=255)
                                                                                                    gATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCGGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCg                                              >  7:308965/1‑149 (MQ=255)
                                                                                                     aTCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCgg                                             <  7:125441/149‑1 (MQ=255)
                                                                                                     aTCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCgg                                             <  7:298845/149‑1 (MQ=255)
                                                                                                        gTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGt                                           <  2:297501/148‑1 (MQ=255)
                                                                                                           cAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCGGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGaaa                                       <  8:308965/149‑1 (MQ=255)
                                                                                                           cAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGaaa                                       <  8:36297/149‑1 (MQ=255)
                                                                                                           cAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGaaa                                       <  2:51557/149‑1 (MQ=255)
                                                                                                                gttcgcgttcgACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCgg                                             <  1:283572/138‑1 (MQ=255)
                                                                                                                gttcgcgttcgACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCgg                                             >  2:283572/1‑138 (MQ=255)
                                                                                                                gttcgcgtgcgACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGc                                  <  5:270844/149‑1 (MQ=255)
                                                                                                                gagcgcgttcgACAATAAACTCCCTCAACGCTTCCACCACGTCCCGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGCCTACTGACGATTCGGGCACTGCAGGCGGGTCAGCGCGCCTTCTCTGCCCCCGGGTGAAACTAGc                                  <  5:196983/146‑1 (MQ=255)
                                                                                                                    gcgttcgACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTaa                              <  7:339137/149‑1 (MQ=255)
                                                                                                                     cgttcgACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAAc                             >  4:94253/1‑149 (MQ=255)
                                                                                                                      gttctACAATAATCTCCATCAACGCTTCCAGCACGTCCAGACGTTTGGATTGGTAGACCCTAACATATCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGc                           <  5:8160/149‑1 (MQ=255)
                                                                                                                      gttcgACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACg                            >  4:119921/1‑149 (MQ=255)
                                                                                                                         cgACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCt                          <  2:131326/148‑1 (MQ=255)
                                                                                                                          gACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTg                         <  3:240705/148‑1 (MQ=255)
                                                                                                                             aaTAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGc                     >  6:204155/1‑149 (MQ=255)
                                                                                                                               tAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTACTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCt                    >  5:148311/1‑148 (MQ=255)
                                                                                                                                    tCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGCcaca              >  5:252169/1‑149 (MQ=255)
                                                                                                                                          aaCGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGc        >  7:100307/1‑149 (MQ=255)
                                                                                                                                               ttCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGc    >  1:48961/1‑148 (MQ=255)
                                                                                                                                               ttCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGc    >  1:43521/1‑148 (MQ=255)
                                                                                                                                                 ccAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGa  <  4:39212/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
AAATCAATGTTTGCCGCAATACCAAACTTTTGCGACAGGGTCATTTGCAGCCACTGTGCCATACCGGTACTTTGCACCAGAATCATCTCTGGTTCGAAAGGATCGTCCAGCCGTTCGCGTTCGACAATAAACTCCATCAACGCTTCCAGCACGTCCAGACGATTGGAATGGTAGACCCTTAACATAGCGGCTCCTGACTACTGACGATTCGGGCAATGCAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGA  >  NC_000913/2962244‑2962536

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: