Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,978,708 C→G 9.0% 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 9.0% 230.7 / 11.6 78F45L (TTC→TTG) cmkcytidylate kinase
Reads supporting (aligned to +/- strand):  ref base C (35/36);  new base G (3/4);  total (38/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.07e-01

AAGAGGGCAAGTCGTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTCAGCAACGT  >  NC_002947/1978561‑1978848
                                                                                                                                                    |                                                                                                                                            
aagaGGGCAAGTCGTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTc                                                                                                                                              >  3:171949/1‑148 (MQ=255)
  gagGGCAAGTCGTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGGCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAAt                                                                                                                                           <  5:128124/149‑1 (MQ=255)
   agGGCAAGTCGTGAATTGATTGGCACCGGTTATCACGATCGAGGGACCAAGCGGCTGGGGCAAGGGCTCGGTCGCGGGGCTGATTGCCCGCGAGCTGGGCTGGTAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATg                                                                                                                                          <  2:123969/149‑1 (MQ=255)
     ggCAAGTCGTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGc                                                                                                                                         >  5:67007/1‑148 (MQ=255)
      gCAAGTCGTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCa                                                                                                                                       <  8:318564/149‑1 (MQ=255)
      gCAAGTCGGGAATTCAATGGCAGCGGTTATCACCATAGACGGACCAAGCTGCTCGGGCAAGGGCTCGGTCGCCGTGCTGATTGCCCGCTAGCTGTGCTGGAAGCTGCT‑GGACTCCTGCGCGCTGTATCGCTTGCTGGCATTCAATGCCa                                                                                                                                       <  7:349742/149‑1 (MQ=255)
        aaGTCGTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTGAATGCCAGc                                                                                                                                     <  1:26867/149‑1 (MQ=255)
          gTCGTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCT‑GGACTCCAGCGCGCTGTATCGCTTGCTGGCATTGAATGCCAGCaa                                                                                                                                   <  3:97829/149‑1 (MQ=255)
             gTGAATTCAATGGCACCGGTTGTCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCa                                                                                                                                >  2:114057/1‑149 (MQ=255)
             gTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTGAATGCCAGCaagc                                                                                                                                 <  8:151902/148‑3 (MQ=255)
             gTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCa                                                                                                                                >  8:268581/1‑149 (MQ=255)
             gTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCa                                                                                                                                >  2:14210/1‑149 (MQ=255)
             gTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCa                                                                                                                                >  3:70867/1‑149 (MQ=255)
              tGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCAc                                                                                                                               >  5:133781/1‑149 (MQ=255)
              tGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCAc                                                                                                                               <  4:171949/149‑1 (MQ=255)
                  ttCAATGGCACCGGTTATCACCATCGACGTACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGCCTGGAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATCCCAGCAACCACGGCg                                                                                                                           >  5:81999/1‑149 (MQ=255)
                          cACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCt                                                                                                                    >  8:79692/1‑148 (MQ=255)
                             cGGTTTTCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCTAGCTGGGCTGGAAGCTGCTCTGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGAc                                                                                                                 <  1:114057/149‑1 (MQ=255)
                             cGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGAc                                                                                                                 >  6:238542/1‑148 (MQ=255)
                                   tCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGaa                                                                                                          <  6:129355/149‑1 (MQ=255)
                                   tCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGaa                                                                                                          <  1:205850/149‑1 (MQ=255)
                                     ctcATCGACGGACCAAGGGGCGCGGGCACTGGCACGGTCGCCGGGCTGATGGCCCGCGAGCTGGGCGGGAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAaga                                                                                                        <  8:333124/147‑1 (MQ=255)
                                     aCCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAaga                                                                                                        <  6:133781/149‑1 (MQ=255)
                                      ccATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCa                                                                                                                                       >  6:167212/1‑117 (MQ=255)
                                      ccATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCa                                                                                                                                       <  5:167212/117‑1 (MQ=255)
                                        aTCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGCCCTGACCAACGAaga                                                                                                        >  7:333124/1‑146 (MQ=255)
                                         tCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAgctg                                                                                                    >  1:371205/1‑149 (MQ=255)
                                            aCGGACCAAGCTGCTCGGGCAAGGGCAGGGTCGCCGGGATGATTGCCCGCGAGCTGGGCTGGAAGCTGCT‑GGTCTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGctgctg                                                                                                 <  6:89498/149‑1 (MQ=255)
                                              ggACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGAc                                                                                               <  4:115817/149‑1 (MQ=255)
                                                   aaGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGCCCAAGc                                                                                          >  7:205702/1‑149 (MQ=255)
                                                    aGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCt                                                                                         >  7:96581/1‑149 (MQ=255)
                                                    aGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCt                                                                                         >  7:70889/1‑149 (MQ=255)
                                                    aGCGGCTCGGGCAAGGGAACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCt                                                                                         >  8:26630/1‑149 (MQ=255)
                                                       ggCGCGGGCAAGGGCACGGTCGCCGGGCGTATTGCCCGCGAGCTGGGCTGGAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTgc                                                                                      <  6:201008/149‑1 (MQ=255)
                                                        gCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTgcc                                                                                     >  2:70538/1‑149 (MQ=255)
                                                          tCGTGCAAGGTCACGGTCTCCGGGCTTATTGCCCGCTAGCTGGGATGGAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGAcc                                                                                              <  7:256533/138‑1 (MQ=255)
                                                          tCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTGAATGCCAGCaagcga                                                                                                                               >  7:226604/1‑101 (MQ=255)
                                                          tCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTGAATGCCAGCaagcga                                                                                                                               <  8:226604/105‑5 (MQ=255)
                                                          tCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGAATTCAATGCCAGAAACCACGGCGTCGACCTGACAAACGAAGAGCTGATGACAAAgatcgtaa                                                                                    >  8:256533/1‑141 (MQ=255)
                                                            gggCAAGGGCAGGGTCGGCGGGCTGTTTGCCGGCGTGCTGGGGTGGAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGccc                                                                                 <  2:371205/149‑1 (MQ=255)
                                                            gggCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTgccgc                                                                                   <  6:67007/147‑1 (MQ=255)
                                                             ggCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCa                                                                                <  7:79692/149‑1 (MQ=255)
                                                              gCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCAt                                                                               >  2:209824/1‑149 (MQ=255)
                                                               cAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATc                                                                              <  1:14210/149‑1 (MQ=255)
                                                                       cGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGAtgtg                                                                      <  5:332999/149‑1 (MQ=255)
                                                                        ggTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGAtgtg                                                                      >  3:208819/1‑148 (MQ=255)
                                                                        ggTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGc                                                                     <  8:96581/149‑1 (MQ=255)
                                                                                ggCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATc                                                             <  5:166345/149‑1 (MQ=255)
                                                                                   tGATTGCCCGCGAGCTGGGCTGGAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTg                                                                            <  8:212768/131‑1 (MQ=255)
                                                                                   tGATTGCCCGCGAGCTGGGCTGGAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTg                                                                            >  7:212768/1‑131 (MQ=255)
                                                                                   tGATTGCCCGCGAGCTGGGCTGGAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCt                                                          <  8:70889/149‑1 (MQ=255)
                                                                                     aTTGCCCGCGAGCTGGGCTGGAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGc                                                        <  4:70867/149‑1 (MQ=255)
                                                                                         cccGCGAGCTGGGCTGGAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAg                                                    <  1:209824/149‑1 (MQ=255)
                                                                                           cgcgAGCTGGGCTGGAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGcc                                                  <  7:268581/149‑1 (MQ=255)
                                                                                                 ctgggctggAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTaa                                            <  3:230974/149‑1 (MQ=255)
                                                                                                  tgggctggAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGAt                                                                         <  7:180977/119‑1 (MQ=255)
                                                                                                  tgggctggAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGAt                                                                         >  8:180977/1‑119 (MQ=255)
                                                                                                    ggctggAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGc                                          <  4:208819/148‑1 (MQ=255)
                                                                                                          aaGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTacaa                                     >  4:337104/1‑147 (MQ=255)
                                                                                                                    gACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCt                          >  2:138508/1‑149 (MQ=255)
                                                                                                                      cTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCTATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCtt                         >  6:230671/1‑148 (MQ=255)
                                                                                                                         cGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGa                       <  5:238542/147‑1 (MQ=255)
                                                                                                                         cGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCACCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAgg                     >  6:89997/1‑149 (MQ=255)
                                                                                                                         cGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTACCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAgg                     <  7:101341/149‑1 (MQ=255)
                                                                                                                                 tGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACGACAGATCATCCTTGAGGGTGAGGAc             >  3:29799/1‑149 (MQ=255)
                                                                                                                                 tGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGa                 >  7:62955/1‑145 (MQ=255)
                                                                                                                                 tGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGa                 <  8:62955/145‑1 (MQ=255)
                                                                                                                                 tGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTGGAGGGTGAGGAc             >  7:179345/1‑149 (MQ=255)
                                                                                                                                    aTCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCtt                         >  5:39060/1‑134 (MQ=255)
                                                                                                                                    aTCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCtt                         <  6:39060/134‑1 (MQ=255)
                                                                                                                                     tCGCTTGCTGGCTTTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCTCTGCCGAGCCCGGTAAGCTACAACTGATCATCCTTGAGGGTGAGGACGTc          <  7:26630/148‑1 (MQ=255)
                                                                                                                                     tCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGat                                >  5:179956/1‑126 (MQ=255)
                                                                                                                                     tCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGat                                <  6:179956/126‑1 (MQ=255)
                                                                                                                                        gttGCTGGCATTGAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTCAGc       >  1:366237/2‑148 (MQ=255)
                                                                                                                                          tGCTGGCATTGAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTCGGCAAc    >  7:174608/1‑149 (MQ=255)
                                                                                                                                           gCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTCAGCAACg   >  1:313360/1‑149 (MQ=255)
                                                                                                                                            cTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTCAGCAACGt  >  3:15350/1‑149 (MQ=255)
                                                                                                                                             tGGCATTCAATGCCAGCAACCCCGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTCAGCAACGt  >  5:82444/1‑148 (MQ=255)
                                                                                                                                             tGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGAAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTCAGCACCGt  >  2:342649/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
AAGAGGGCAAGTCGTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCT‑GGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTCAGCAACGT  >  NC_002947/1978561‑1978848

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