Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,978,708 C→G 12.7% 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,7080CG12.7% 173.4 / 16.6 63F45L (TTC→TTG) cmkcytidylate kinase
Reads supporting (aligned to +/- strand):  ref base C (29/26);  new base G (4/4);  total (33/30)
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.84e-01

AGTCGTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTCAGCAACGTCATCCG  >  NC_002947/1978570‑1978854
                                                                                                                                          |                                                                                                                                                  
aGTCGTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGTCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTGAATGCCAGCa                                                                                                                                          <  7:278961/149‑1 (MQ=255)
aGTCGTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCa                                                                                                                                          <  3:355482/149‑1 (MQ=255)
 ggcGTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCGCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTGAATGCCAGCaa                                                                                                                                         <  3:236983/147‑1 (MQ=255)
 gTCGTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTGAATGCCAGCaa                                                                                                                                         <  8:265726/149‑1 (MQ=255)
     tGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCAc                                                                                                                                     >  3:14005/1‑149 (MQ=255)
       aaTTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACgg                                                                                                                                   >  8:186309/1‑149 (MQ=255)
         ttCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTCTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCg                                                                                                                                 >  7:213817/1‑149 (MQ=255)
            aaTGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTc                                                                                                                               >  6:379121/1‑148 (MQ=255)
              tGGCCCCGTGTATCACCATCGACGGACCAAGCGGCGCGGGCAAGGTCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGGGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGAc                                                                                                                            <  8:191143/149‑1 (MQ=255)
              tGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGAc                                                                                                                            <  2:10450/149‑1 (MQ=255)
                 cACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCt                                                                                                                          >  3:85399/1‑148 (MQ=255)
                  aCCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGa                                                                                                                        <  5:379121/149‑1 (MQ=255)
                      gTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTGAATGCCAGCaagc                                                                                                                                       >  6:78210/1‑128 (MQ=255)
                      gTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTGAATGCCAGCaagc                                                                                                                                       <  5:78210/130‑3 (MQ=255)
                         aTCACCATCGCCGGACCATGCGGGTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGGTGGTTTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGGAACCACGGCGTCTACCTGACCAACg                                                                                                                  <  6:27879/148‑1 (MQ=255)
                         aTCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACg                                                                                                                  >  8:300462/1‑148 (MQ=255)
                          tCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTTCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGaa                                                                                                                <  2:3803/149‑1 (MQ=255)
                          tCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGaa                                                                                                                <  8:13730/149‑1 (MQ=255)
                              cATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAgc                                                                                                            >  8:2854/1‑149 (MQ=255)
                               aTCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAgct                                                                                                           >  6:338717/1‑149 (MQ=255)
                               aTCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGACGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAgct                                                                                                           >  1:71062/1‑149 (MQ=255)
                                    cGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGa                                                                                                      <  7:2854/149‑1 (MQ=255)
                                     ggACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGAc                                                                                                     >  5:312439/1‑149 (MQ=255)
                                            gCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGaa                                                                                                                <  5:180387/131‑1 (MQ=255)
                                            gCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGaa                                                                                                                >  6:180387/1‑131 (MQ=255)
                                            gCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCtt                                                                                              <  5:338717/149‑1 (MQ=255)
                                                cTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTgccg                                                                                          <  6:312439/149‑1 (MQ=255)
                                                 tCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCTGGAGCTGTGCTGGAATCTGCTGGACTCGGGCGCGCTGTTTCTCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCTAGGTAGAGGTGCTGACCTAGCTTgccga                                                                                         <  6:272763/149‑2 (MQ=255)
                                                 tCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCATGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTgccgc                                                                                         >  4:99834/1‑149 (MQ=255)
                                                 tCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTgccgc                                                                                         >  5:122251/1‑149 (MQ=255)
                                                   gggCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTgccgcc                                                                                        >  7:71185/1‑148 (MQ=255)
                                                    ggCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTGAATGCCAGCaag                                                                                                                                        >  8:119130/1‑98 (MQ=255)
                                                       aaGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTCTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCt                                                                                   <  8:213817/149‑1 (MQ=255)
                                                         gggCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTgg                                                                                 <  8:71185/149‑1 (MQ=255)
                                                         gggCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTg                                                                                  >  1:70985/1‑148 (MQ=255)
                                                          ggCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGa                                                                                >  7:25043/1‑149 (MQ=255)
                                                           gCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGAt                                                                               <  7:186309/149‑1 (MQ=255)
                                                            cACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGAtg                                                                              <  8:25043/149‑1 (MQ=255)
                                                               ggTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGAt                                                                               <  1:214973/145‑1 (MQ=255)
                                                               ggTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGAt                                                                               >  2:214973/1‑145 (MQ=255)
                                                                    ccGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGtt                                                                       >  2:75234/1‑148 (MQ=255)
                                                                      gggCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCAt                                                                    <  7:300462/149‑1 (MQ=255)
                                                                                cccGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAg                                                          <  4:85399/149‑1 (MQ=255)
                                                                                        ctgggctggAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGt                                                    >  2:84704/1‑147 (MQ=255)
                                                                                             ctggAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTa                                              >  7:366341/1‑148 (MQ=255)
                                                                                              tggAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTac                                             >  7:207423/1‑148 (MQ=255)
                                                                                                gAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGa                                                           <  1:296183/132‑1 (MQ=255)
                                                                                                gAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGa                                                           >  2:296183/1‑132 (MQ=255)
                                                                                                   gctgctGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCATGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGa                                       <  3:99834/149‑1 (MQ=255)
                                                                                                          gACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCTACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCt                                >  3:391154/1‑149 (MQ=255)
                                                                                                           aCTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCtt                               <  2:71062/149‑1 (MQ=255)
                                                                                                               cGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAgg                           >  3:323937/1‑149 (MQ=255)
                                                                                                                           tCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTCa               <  1:84704/149‑1 (MQ=255)
                                                                                                                           tCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTCa               <  4:14005/149‑1 (MQ=255)
                                                                                                                            cGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGa                                       <  6:212852/124‑1 (MQ=255)
                                                                                                                            cGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGa                                       >  5:212852/1‑124 (MQ=255)
                                                                                                                               ttGCTGGCATTGAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTCAGCaa           >  1:190593/1‑149 (MQ=255)
                                                                                                                               ttGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTCAGCaa           >  4:36649/1‑149 (MQ=255)
                                                                                                                                   tGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGAtgt                                                                             <  3:340814/79‑1 (MQ=255)
                                                                                                                                   tGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGAtgt                                                                             >  4:340814/1‑79 (MQ=255)
                                                                                                                                     gCATTGAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTCAGCAACGTCAt     >  6:177982/1‑149 (MQ=255)
                                                                                                                                     gCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTCAGCAACGTCAt     >  1:377330/1‑149 (MQ=255)
                                                                                                                                        ttCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTCAGCAACGTCATCCg  >  2:133910/1‑149 (MQ=255)
                                                                                                                                          |                                                                                                                                                  
AGTCGTGAATTCAATGGCACCGGTTATCACCATCGACGGACCAAGCGGCTCGGGCAAGGGCACGGTCGCCGGGCTGATTGCCCGCGAGCTGGGCTGGAAGCTGCTGGACTCCGGCGCGCTGTATCGCTTGCTGGCATTCAATGCCAGCAACCACGGCGTCGACCTGACCAACGAAGAGCTGCTGACCAAGCTTGCCGCCCATCTGGATGTGCAGTTCATCGCTGCCGAGCCCGGTAAGCTACAACAGATCATCCTTGAGGGTGAGGACGTCAGCAACGTCATCCG  >  NC_002947/1978570‑1978854

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