Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,440,363 T→C 8.2% T155A (ACC→GCC)  amn ← AMP nucleosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,440,3630TC 8.2% 168.9 / 5.2 61T155A (ACC→GCC) amnAMP nucleosidase
Reads supporting (aligned to +/- strand):  ref base T (33/23);  new base C (2/3);  total (35/26)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.42e-01

GGATCCACGGCTGTACATGGCGCCAGTCGCTGCCGGTGTAGTGCACCAGGCGGCGTAGCGAGAAGTCGACGCGCGCGGCATCGAACAACGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGTAAGGGATATTCTGCTCGCTGACGCCGACTTCTACCGTAACGGTGTAGT  >  NC_002947/5440216‑5440509
                                                                                                                                                   |                                                                                                                                                  
ggATCCACGGCTGTACATGGCGCCAGTCGCTGCCGGTGTAGTGCACCAGGCGGCGTAGCGAGAAGTCGACGCGCGCGGCATCGAACAACGCCAGTGGCATGATTTCAGCACTGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTg                                                                                                                                                   <  4:323931/149‑1 (MQ=255)
 gATCCACGGCTGTACATGGCGCCAGTCGCTGCCGGTGTAGTGCACCAGGCGGCGTAGCGAGAAGTCGACGCGCGCGGCATCGAACAACGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTgg                                                                                                                                                  >  6:256289/1‑149 (MQ=255)
      aCGGCTGTACATGGCGCCAGTCGCTGCCGGTGTAGTGCACCAGGCGGCGTAGCGAGAAGTCGACGCGCGCGGCATCGAACAACGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGgccgcc                                                                                                                                             >  2:193970/1‑149 (MQ=255)
             tACATGGCGCCAGTCGCTGCCGGTGTAGTGCACCAGGCGGCGTAGCGAGAAGTCGACGCGCGCGGCATCGAACAACGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGCGGCCACCGACaaca                                                                                                                                      <  1:404292/149‑4 (MQ=255)
               cATGGCGCCAGTCGCTGCCGGTGTAGTGCACCAGGCGGCGTAGCGAGAAGTCGACGCGCGCGGCATCGAACAACGCCAGGGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGCGGCCACCGACaacatc                                                                                                                                    <  4:16158/149‑6 (MQ=255)
                     gcCAGTCGCTGCCGGTGTAGTGCACCAGGCGGCGTAGCGAGAAGTCGACGCGCGCGGCATCGAACAACGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTg                                                                                                                              >  4:128536/1‑149 (MQ=255)
                      ccAGTCGCTGCCGGTGTAGTGCACCAGGCGGCGTAGCGAGAAGTCGACGCGCGCGGCATCGAACAACGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTgg                                                                                                                             <  3:235001/149‑1 (MQ=255)
                       cAGTCGCTGCCGGTGTAGTGCACCAGGCGGCGTAGCGAGAAGTCGACGCGCGCGGCATCGAACAACGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTggg                                                                                                                            <  5:272096/149‑1 (MQ=255)
                                ccGGTGTAGTGCACCAGGCGGCGTAGCGAGAAGTCGACGCGCGCGGCATCGAACAACGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGAGCTGGGGAACACCcg                                                                                                                   >  3:301220/1‑149 (MQ=255)
                                        gTGCACCAGGCGGCGTAGCGAGAAGTCGACGCGCGCGGCATCGAACAACGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCt                                                                                                           >  8:307839/1‑149 (MQ=255)
                                        gTGCACCAGGCGGCGTAGCGAGAAGTCGACGCGCGCGGCATCGAACAACGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCt                                                                                                           >  3:369893/1‑149 (MQ=255)
                                         tGCACCAGGCGGCGTAGCGAGAAGTCGACGCGCGCGGCATCGAACAACGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTg                                                                                                          >  2:248383/1‑149 (MQ=255)
                                             ccAGGCGGCGTAGCGAGAAGTCGACGCGCGCGGCATCGAACAACGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTgcgc                                                                                                       >  1:64754/1‑148 (MQ=255)
                                             ccAGGCGGCGTAGCGAGAAGTCGACGCGCGCGGCATCGAACAACGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTgcgc                                                                                                       <  2:64754/148‑1 (MQ=255)
                                                 gcggcgTAGCGAGAAGTCGACGCGCGCGGCATCGAACAACGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCcgc                                                                                                                  <  5:138927/133‑1 (MQ=255)
                                                 gcggcgTAGCGAGAAGTCGACGCGCGCGGCATCGAACAACGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCcgc                                                                                                                  >  6:138927/1‑133 (MQ=255)
                                                       tAGCGAGAAGTCGACGCGCGCGGCATCGAACAACGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGc                                                                                            >  3:104224/1‑149 (MQ=255)
                                                           gagaAGTCGACGCGCGCGGCATCGAACAACGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACAcc                                                                                                                     >  7:46219/1‑120 (MQ=255)
                                                           gagaAGTCGACGCGCGCGGCATCGAACAACGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACAcc                                                                                                                     <  8:46219/120‑1 (MQ=255)
                                                               aGTCGACGCGCGCGGCATCGAACAACGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCg                                                                                    <  2:347197/149‑1 (MQ=255)
                                                                        gcgcgGCATCGAACAACGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCtcg                                                                           <  7:307839/149‑1 (MQ=255)
                                                                        gcgcgGCATCGAACAACGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCtcg                                                                           >  4:143051/1‑149 (MQ=255)
                                                                                   aacaacGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTc                                                                >  7:36797/1‑149 (MQ=255)
                                                                                        cGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGGCca                                                           <  3:143051/149‑1 (MQ=255)
                                                                                          ccAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGAcc                                                            >  2:52700/1‑146 (MQ=255)
                                                                                          ccAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGAcc                                                            <  1:52700/146‑1 (MQ=255)
                                                                                            aGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGACCAGGGCGAcgagctg                                                       <  1:193970/149‑7 (MQ=255)
                                                                                                 cATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAgggg                                                   >  6:386015/1‑148 (MQ=255)
                                                                                                  aTGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGt                                                  >  4:326477/1‑148 (MQ=255)
                                                                                                   tGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGt                                                  >  6:284490/1‑147 (MQ=255)
                                                                                                    gATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGACGCCCTGGTCGACCACATAGACGTAAg                                               >  8:278089/1‑149 (MQ=255)
                                                                                                           gCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGTAAGGGATAtt                                        <  3:166074/149‑1 (MQ=255)
                                                                                                            cACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGGCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGTAAGGGATATTc                                       >  1:285698/1‑149 (MQ=255)
                                                                                                             aCGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGTAAGGGATATTCt                                      <  1:248383/149‑1 (MQ=255)
                                                                                                                 tCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGTAAGGGATATTCTGccc                                  >  5:165357/1‑147 (MQ=255)
                                                                                                                 tCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGTAAGGGATATTCTGCTc                                  >  7:158177/1‑149 (MQ=255)
                                                                                                                  cccAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGTAAGGGATATTCTGCTCg                                 >  5:418236/1‑149 (MQ=255)
                                                                                                                   ccAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACca                                                           >  6:263931/1‑122 (MQ=255)
                                                                                                                   ccAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACca                                                           <  5:263931/122‑1 (MQ=255)
                                                                                                                   ccAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGTAAGGGATATTCTGCTCGc                                >  4:212538/1‑149 (MQ=255)
                                                                                                                          tACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGTAAGGGATATTCTGCTCGCTGACGcc                         >  1:270999/1‑149 (MQ=255)
                                                                                                                            cAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGTAAGGGATATTCTGCTCGCTGACGCCGa                       <  3:99887/149‑1 (MQ=255)
                                                                                                                               gCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGTAAGGGATATTCTGCTCGCTGACGCCGACtt                    >  5:205154/1‑149 (MQ=255)
                                                                                                                                 cATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGTAAGGGATATTCTGCTCGCTGACGCCGACTTCt                  <  8:158177/149‑1 (MQ=255)
                                                                                                                                  aTCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGTAAGGGATATTCTGCTc                                  <  5:227829/132‑1 (MQ=255)
                                                                                                                                  aTCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGTAAGGGATATTCTGCTc                                  >  6:227829/1‑132 (MQ=255)
                                                                                                                                   tCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTCGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCAACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGTAAGGGATATTCTGCTCGCTGACGCCGACTTCTAc                >  8:127693/1‑149 (MQ=255)
                                                                                                                                     aGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGTAAGGGATATTCTGCTCGCTGACGCCGACTTCTAcc               <  3:326477/148‑1 (MQ=255)
                                                                                                                                     aGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGTAAGGGATATTCTGCTCGCTGACGCCGACTTCTAcc               <  3:212538/148‑1 (MQ=255)
                                                                                                                                        gATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGTAAGGGATATTCTGCTCGCTGACGCCGACTTCTACCGTaa           <  3:128536/149‑1 (MQ=255)
                                                                                                                                        cgacctgtCGGCGGCCACCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGTAAGGGATATTCTGCTCGCTGACGCCGACTTCTAc                <  6:358013/139‑1 (MQ=255)
                                                                                                                                        cgacctgtCGGCGGCCACCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGTAAGGGATATTCTGCTCGCTGACGCCGACTTCTAc                >  5:358013/6‑144 (MQ=255)
                                                                                                                                          tgttgtCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGTAAGGGATATTCTGCTCGCTGACGCCGACTTCtaccgtaccg         >  1:385983/1‑149 (MQ=255)
                                                                                                                                          tgttgtCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGTAAGGGATATTCTGCTCGCTGACGCCGACTTCTACCGTAACg         <  2:270999/149‑1 (MQ=255)
                                                                                                                                          acctgtCGGCGGCCACCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGTAAGGGATATTCTGCTCGCTGACGCCGACTTCTACCGTAACg         >  2:3518/4‑149 (MQ=255)
                                                                                                                                               tCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAggg                                                    >  1:254665/1‑101 (MQ=255)
                                                                                                                                               tCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAggg                                                    <  2:254665/101‑1 (MQ=255)
                                                                                                                                               tCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGTAAGGGATATTCTGCTCGCTGACGCCGACTTCTACCGTAACGgtg      >  7:297319/1‑147 (MQ=255)
                                                                                                                                                  gTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTACTGCCGCTACCGGCCGGCTCGTCGCCGTGGTCGACCAAATAGGGGTAAGGGATATTCTGCTCGCTGACGCCGACTTCTTCCGTAACGGTGTAGt  >  5:4033/1‑148 (MQ=255)
                                                                                                                                                  gTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGTAAGGGATATTCTGCTCGCTGACGc                          >  3:253193/1‑124 (MQ=255)
                                                                                                                                                  gTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGTAAGGGATATTCTGCTCGCTGACGc                          <  4:253193/124‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GGATCCACGGCTGTACATGGCGCCAGTCGCTGCCGGTGTAGTGCACCAGGCGGCGTAGCGAGAAGTCGACGCGCGCGGCATCGAACAACGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGTAAGGGATATTCTGCTCGCTGACGCCGACTTCTACCGTAACGGTGTAGT  >  NC_002947/5440216‑5440509

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