Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,978,708 C→G 7.1% F45L (TTC→TTG)  cmk → cytidylate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,978,7080CG 7.1% 216.1 / 3.4 67F45L (TTC→TTG) cmkcytidylate kinase
Reads supporting (aligned to +/- strand):  ref base C (29/33);  new base G (3/2);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAGAGGGCAAGTCGTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAA‑GCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTCAGCAACGTCATCC  >  NC_002947/1978560‑1978853
                                                                                                                                                     |                                                                                                                                                 
gaagaGGGCAAGTCGTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAA‑GCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTc                                                                                                                                                   <  4:115169/149‑1 (MQ=255)
      ggCAAGTCGTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAA‑GCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGcc                                                                                                                                             <  4:346163/149‑1 (MQ=255)
       gCAAGTCGTGAATTCAATGGCACCGGTTATCACCATCGACGGACCACGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAA‑GCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCa                                                                                                                                            <  8:268991/149‑1 (MQ=255)
       gCAAGTCGTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAA‑GCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGc                                                                                                                                              >  8:119126/1‑147 (MQ=255)
                 aaTTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAA‑GCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTGAATGCCAGCaag                                                                                                                                       >  5:341085/1‑143 (MQ=255)
                 aaTTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAA‑GCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTGAATGCCAGCaag                                                                                                                                       <  6:341085/144‑2 (MQ=255)
                   ttCAATGGCACCGGTTATCACCATCGACGGACCATGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAACGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGc                                                                                                                                 <  1:106862/149‑1 (MQ=255)
                      aaTGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAA‑GCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCg                                                                                                                             >  3:291901/1‑149 (MQ=255)
                      aaTGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAA‑GCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCg                                                                                                                             >  4:41521/1‑149 (MQ=255)
                           cACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAA‑GCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCt                                                                                                                         >  6:45887/1‑148 (MQ=255)
                            aCCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAA‑GCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGa                                                                                                                       >  4:92997/1‑149 (MQ=255)
                                     cACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAA‑GCTGCTGGACTCCGGCGCGCTCTATCGCTTGCTGGAATTCAATGCCAGCAACCACGGCGTCGACATGACCAACg                                                                                                                 >  7:124594/1‑146 (MQ=255)
                                        cATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAA‑GCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAgc                                                                                                           >  3:293997/1‑149 (MQ=255)
                                         aTCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAA‑GCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAgct                                                                                                          >  6:74360/1‑149 (MQ=255)
                                             aCGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAA‑GCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAag                                                                                                              >  6:104735/1‑141 (MQ=255)
                                             aCGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAA‑GCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAag                                                                                                              <  5:104735/141‑1 (MQ=255)
                                             aCGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAA‑GCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAgctgct                                                                                                       >  3:313515/1‑148 (MQ=255)
                                              cGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAA‑GCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGa                                                                                                     <  4:313515/149‑1 (MQ=255)
                                              cGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAA‑GCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGa                                                                                                     <  7:112903/149‑1 (MQ=255)
                                                gACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAA‑GCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGAcc                                                                                                   <  7:28619/149‑1 (MQ=255)
                                                 aCCAAGCGGCACCGGCAAGGGCACGGTCGCAGGGCTGATTGCCCGCGAGATGGGCTGGAA‑GCTGCTGGACTCAGGCGCGCTGTATAGATTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCAGACAAACGAAGAGCTGCTGAcc                                                                                                   >  6:241120/1‑148 (MQ=255)
                                                    aaGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAA‑GCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGc                                                                                               <  5:75820/149‑1 (MQ=255)
                                                        ggCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAA‑GCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTgc                                                                                           <  3:41521/149‑1 (MQ=255)
                                                        ggCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAA‑GCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTgc                                                                                           <  2:7138/149‑1 (MQ=255)
                                                         gCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAA‑GCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTgcc                                                                                          <  3:6442/149‑1 (MQ=255)
                                                         gCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAA‑GCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTgcc                                                                                          >  5:71250/1‑149 (MQ=255)
                                                                 aaGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAA‑GCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCt                                                                                  <  6:291349/149‑1 (MQ=255)
                                                                   gggCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAA‑GCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTg                                                                                 >  7:302441/1‑148 (MQ=255)
                                                                      cACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAA‑GCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGAtg                                                                             >  6:384186/1‑149 (MQ=255)
                                                                         ggTCGCGGGGCTGATTGCCCGCGAGCTTGGCTGGAA‑GCTGCTGGACTCCGGCGGGCTGTATCGGTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGAtgtg                                                                           <  4:348673/148‑1 (MQ=255)
                                                                               cGGGCTGATTGCCCGCGAGCTGGGCTGGAA‑GCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCa                                                                    <  4:291901/149‑1 (MQ=255)
                                                                               cGGGCTGATTGCCCGCGAGCTGGGCTGGAA‑GCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCa                                                                    <  8:302441/149‑1 (MQ=255)
                                                                                gggCTGATTGCCCGCGAGCTGGGCTGGAA‑GCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGCGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCAt                                                                   >  7:317714/1‑149 (MQ=255)
                                                                                 ggCTGATTGCCCGCGAGCTGGGCTGGAA‑GCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATc                                                                  <  8:45169/149‑1 (MQ=255)
                                                                                    tgtttGCCCGCGAGCTGGTCTGGAA‑GCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTGAATGCCAGCAACCACGGCGTCGACGTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCt                                                               <  3:114227/146‑1 (MQ=255)
                                                                                      aTTGCCCGCGAGCTGGGCTGTAA‑GCGGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGACAGCAACCACGGCGTCGACCTGACCAACGAAGATCTGCTGACCAAGCTTGCCGCCCATCTGTATGTGCAGTTCATCGCTGc                                                             <  8:124594/149‑1 (MQ=255)
                                                                                       ttGCCCGCGAGCTGGGCTGGAA‑GCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCAGCCCATCTGGATGTGCAGTTCATCGCTGcc                                                            >  7:247988/1‑149 (MQ=255)
                                                                                        tGCCCGCGAGCTGGGCTGGAA‑GCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGcc                                                            <  8:167702/148‑1 (MQ=255)
                                                                                          cccGCGAGCTGGGCTGGAA‑GCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGa                                                          <  4:293997/148‑1 (MQ=255)
                                                                                          cccGCGAGCTGGGCTGGAA‑GCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAg                                                         >  7:101046/1‑149 (MQ=255)
                                                                                             gcgAGCTGGGCTGGAA‑GCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCTGCAACCACGGCGTCGACCTGACCa                                                                                                                    <  7:233390/87‑1 (MQ=255)
                                                                                             gcgAGCTGGGCTGGAA‑GCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCa                                                                                                                    >  8:233390/1‑87 (MQ=255)
                                                                                             gcgAGCTGGGCTGGAA‑GCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGccc                                                      <  7:119126/149‑1 (MQ=255)
                                                                                              cgAGCTGGGCTGGAA‑GCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCg                                                     <  5:74360/149‑1 (MQ=255)
                                                                                                  ctgggctggAA‑GCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGt                                                   >  3:207371/1‑147 (MQ=255)
                                                                                                    gggctggAA‑GCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGc                                                        <  4:207541/140‑1 (MQ=255)
                                                                                                    gggctggAA‑GCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGc                                                        >  3:207541/1‑140 (MQ=255)
                                                                                                      gctggAA‑GCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTa                                             <  6:71250/149‑1 (MQ=255)
                                                                                                       ctggAA‑GCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTa                                             >  7:289027/1‑148 (MQ=255)
                                                                                                       ctggAA‑GCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTa                                             >  6:122953/1‑148 (MQ=255)
                                                                                                               ctgctgGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTacaa                                          >  5:252041/1‑144 (MQ=255)
                                                                                                                      aCTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCtt                              <  5:384186/149‑1 (MQ=255)
                                                                                                                        tCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGcc                                                       >  1:286527/1‑122 (MQ=255)
                                                                                                                        tCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGcc                                                       <  2:286527/122‑1 (MQ=255)
                                                                                                                           ggCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAggg                         >  6:257006/1‑149 (MQ=255)
                                                                                                                             cgcgcTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTg                       <  3:92997/149‑1 (MQ=255)
                                                                                                                                    tATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGt                <  5:45887/149‑1 (MQ=255)
                                                                                                                                     aTCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGt                >  2:236738/1‑148 (MQ=255)
                                                                                                                                      tCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAgg                          >  6:184262/1‑137 (MQ=255)
                                                                                                                                      tCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTCa              <  8:289027/149‑1 (MQ=255)
                                                                                                                                      tCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCCCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGATGGTGAGTACGTc               >  8:302331/1‑148 (MQ=255)
                                                                                                                                       cggttGGTGGGATTCTATGCCAGCAACCTCGGAGTGGACCTGTCCAACTAAGAGCTGCTGACCAAGATTGCCGCCCTTCTGGATGTGCAGTTCATCGCTGGCGAGCCCGGTAAGCTACAAAAGATCATCGTTGAGGGTGAGGACGTCa              <  5:241120/145‑1 (MQ=255)
                                                                                                                                        gCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTctg             <  5:122953/148‑3 (MQ=255)
                                                                                                                                           tGCTGGCATTGAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTCAGCAAc         >  8:47215/1‑149 (MQ=255)
                                                                                                                                             cacgCACTCAAGGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTCAGCAACGt       <  8:247988/146‑1 (MQ=255)
                                                                                                                                                  aTTGAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGCCCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTCAGCAACGTCATcc  >  2:282258/1‑149 (MQ=255)
                                                                                                                                                    tCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAg                     <  6:146035/128‑1 (MQ=255)
                                                                                                                                                    tCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAg                     >  5:146035/1‑128 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GAAGAGGGCAAGTCGTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAA‑GCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTCAGCAACGTCATCC  >  NC_002947/1978560‑1978853

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