Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,440,368 G→A 6.8% A153V (GCG→GTG)  amn ← AMP nucleosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,440,3680GA 6.8% 173.6 / 2.0 59A153V (GCG→GTG) amnAMP nucleosidase
Reads supporting (aligned to +/- strand):  ref base G (28/27);  new base A (2/2);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.14e-01

CTGTACATGGCGCCAGTCGCTGCCGGTGTAGTGCACCAGGCGGCGTAGCGAGAAGTCGACGCGCGCGGCATC‑GAACAACGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGTAAGGGATATTCTGCTCGCTGACGCCGACTTCTACCGTAACGGTGTAGTCATG  >  NC_002947/5440226‑5440513
                                                                                                                                               |                                                                                                                                                 
cTGTACATGGCGCCAGTCGCTGCCGGTGTAGTGCACCAGGCGGCGTAGCGAGAAGTCGACGCGCGCGGCATC‑GAACAACGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACa                                                                                                                                             <  6:322133/149‑1 (MQ=255)
     cATGGCGCCAGTCGCTGCCGGTGTAGTGCACCAGGCGGCGTAGCGAGAAGTCGACGCGCGCGGCATC‑GAACAACGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCg                                                                                                                                        >  8:103584/1‑149 (MQ=255)
        ggCGCCAGTCGCTGCCGGTGTAGTGCACCAGGCGGCGTAGCGAGAAGTCGACGCGCGCGGCATC‑GAACAACGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTg                                                                                                                                     >  3:36552/1‑149 (MQ=255)
         gcgcCAGTCGCTGCCGGTGTAGTGCACCAGGCGGCGTAGCGAGAAGTCGACGCGCGCGGCATC‑GAACAACGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTg                                                                                                                                     <  4:36552/148‑1 (MQ=255)
          cgcCAGTCGCTGCCGGTGTAGTGCACCAGGCGGCGTAGCGAGAAGTCGACGCGCGCGGCATC‑GAACAACGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGCGGCCACCGACaacat                                                                                                                                         >  1:170734/1‑139 (MQ=255)
          cgcCAGTCGCTGCCGGTGTAGTGCACCAGGCGGCGTAGCGAGAAGTCGACGCGCGCGGCATC‑GAACAACGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGCGGCCACCGACaacat                                                                                                                                         <  2:170734/143‑5 (MQ=255)
             cAGTCGCTGCCGGTGTAGTGCACCAGGCGGCGTAGCGAGAAGTCGACGCGCGCGGCATC‑GAACAACGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTggg                                                                                                                                <  5:207806/149‑1 (MQ=255)
             cAGTCGCTGCAGGTGTAGTGCACCAGGCGTCGTAGCGAGAAGTCGAGGCGCGCGGCATC‑GAACAACGCGAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTggg                                                                                                                                <  6:376772/149‑1 (MQ=255)
                tCGCTGCCGGTGTAGTGCACCAGGCGGCGTAGCGAGAAGTCGACGCGCGCGGCATC‑GAACAACGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGaa                                                                                                                             >  8:143234/1‑149 (MQ=255)
                           gtAGTGCACCAGGCGGCGTAGCGAGAAGTCGACGCGCGCGGCATC‑GAACAACGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCgcgc                                                                                                                    >  7:142079/1‑147 (MQ=255)
                              gTGCACCAGGCGGCGTAGCGAGAAGTCGACGCGCGCGGCATC‑GAACAACGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCt                                                                                                               >  6:280148/1‑149 (MQ=255)
                              gTGCACCAGGCGGCGTAGCGAGAAGTCGACGCGCGCGGCATC‑GAACAACGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCt                                                                                                               >  2:202651/1‑149 (MQ=255)
                                 cACCAGGCGGCGTAGCGAGAAGTCGACGCGCGCGGCATC‑GAACAACGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTgc                                                                                                             <  8:142079/148‑1 (MQ=255)
                                     aGGCGGCGTAGCGAGAAGTCGACGCGCGCGGCATCTGAACAACGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTgg                                                                                                                                 >  1:56941/1‑125 (MQ=255)
                                     aGGCGGCGTAGCGAGAAGTCGACGCGCGCGGCATC‑GAACAACGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTgg                                                                                                                                 <  2:56941/124‑1 (MQ=255)
                                       gcggcgTAGCGAGAAGTCGACGCGCGCGGCATC‑GAACAACGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCgg                                                                                                      <  7:103584/149‑1 (MQ=255)
                                                 gagaAGTCGACGCGCGCGGCATC‑GAACAACGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGc                                                                                             >  5:359931/1‑148 (MQ=255)
                                                   gaAGTCGACGCGCGCGGCATC‑GAACAACGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTAc                                                                                          >  3:240280/1‑149 (MQ=255)
                                                         gACGCGCGCGGCATC‑GAACAACGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGc                                                                                      >  6:28258/1‑147 (MQ=255)
                                                         gACGCGCGCGGCATC‑GAACAACGCCAGTGGCATGATTTCAGCACGGTCCCAGGCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGc                                                                                      >  8:3001/1‑147 (MQ=255)
                                                              gcgcgGCATC‑GAACAACGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCtcg                                                                               >  7:97774/1‑149 (MQ=255)
                                                                gcgGCATC‑GAACAACGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCtcgtc                                                                             <  5:28258/149‑1 (MQ=255)
                                                                gcgGCATC‑GAACAACGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCtcgtc                                                                             <  6:359931/149‑1 (MQ=255)
                                                                     aTC‑GAACAACGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCt                                                                        >  4:332853/1‑149 (MQ=255)
                                                                      tC‑GAACAACGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGAGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTg                                                                       >  1:265923/1‑149 (MQ=255)
                                                                         gAACAACGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTgg                                                                      <  5:232149/148‑1 (MQ=255)
                                                                         gAACAACGCCAGTGGCATGATTTCAGCACGGTCCCAGGCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTgg                                                                      <  7:3001/148‑1 (MQ=255)
                                                                            caacGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGa                                                                  <  3:61458/149‑1 (MQ=255)
                                                                             aacGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGAc                                                                 <  1:32646/149‑1 (MQ=255)
                                                                               cGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACca                                                               >  6:123026/1‑149 (MQ=255)
                                                                                    gTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCCGCTCGTCGCCCTGGTCGACCACATAg                                                          <  1:202651/149‑1 (MQ=255)
                                                                                    gTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAg                                                          <  1:162317/149‑1 (MQ=255)
                                                                                     tGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAgg                                                         <  4:192468/149‑1 (MQ=255)
                                                                                     tGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAgg                                                         <  5:280148/149‑1 (MQ=255)
                                                                                       gCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAggg                                                        <  3:332853/148‑1 (MQ=255)
                                                                                        cATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAgggg                                                       >  1:47680/1‑148 (MQ=255)
                                                                                            aTTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGTAAgg                                                  <  8:97774/149‑1 (MQ=255)
                                                                                            aTTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGTAAg                                                   <  2:47680/148‑1 (MQ=255)
                                                                                            aTTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGAGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGTAAg                                                   <  2:265923/148‑1 (MQ=255)
                                                                                             tttCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGTAAgg                                                  >  1:370310/1‑148 (MQ=255)
                                                                                                 aGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCt                                                                        >  5:27637/1‑122 (MQ=255)
                                                                                                 aGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCt                                                                        <  6:27637/122‑1 (MQ=255)
                                                                                                        tCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGTAAGGGATATTCTGCTc                                      <  2:370310/149‑1 (MQ=255)
                                                                                                         cccAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGTAAGGGATATTCTGCTCg                                     >  5:347155/1‑149 (MQ=255)
                                                                                                         cccAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGTAAGGGATATTCTGCTCg                                     >  3:72044/1‑149 (MQ=255)
                                                                                                              tCGTACAGGCCATCAGTGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAggg                                                        <  5:231757/125‑1 (MQ=255)
                                                                                                              tCGTACAGGCCATCAGTGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAggg                                                        >  6:231757/1‑125 (MQ=255)
                                                                                                                 tACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGTAAGGGATATTCTGCTCGCTGACGcc                             >  2:161482/1‑149 (MQ=255)
                                                                                                                 tACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGT‑ATGCCGCGACCGGCCAGCAACGCCCCCGGGGACCCCCCATAAGGGGAAGGGGTTTTTAGCTAGCTGACGCCg                            >  4:113750/1‑149 (MQ=255)
                                                                                                                    aGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGTAAGGGATATTCTGCTCGCTGACGCCGAc                          <  7:143234/149‑1 (MQ=255)
                                                                                                                            aGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGTAAGGGATATTCTGCTCGCTGACGCCGACTTCTACCg                  >  3:254492/1‑149 (MQ=255)
                                                                                                                             gCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGTAAGGGATATTCTGCTCGCTGACGCCGACTTCTACCg                  >  4:85291/1‑148 (MQ=255)
                                                                                                                               gATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGTAAGGGATATTCTGCTCGCTGACGCCGACTTCTACCGTaa               >  8:230760/1‑149 (MQ=255)
                                                                                                                                 tgttgtCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGTAAGGGATATTCTGCTCGCTGACGCCGACTTCTACCGTAACg             <  4:72044/149‑1 (MQ=255)
                                                                                                                                   ctgtCGGCGGCCACCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGTAAGGGATATTCTGCTCGCTGACGCCGAc                          >  6:289433/2‑134 (MQ=255)
                                                                                                                                   ctgtCGGCGGCCACCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGTAAGGGATATTCTGCTCGCTGACGCCGAc                          <  5:289433/133‑1 (MQ=255)
                                                                                                                                        ggTGGACGCCGACAGGTAGGTGATGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGTAAGGGATATTCTGCTCGCTGACGCCGACTTCTACCGTAACGGTGTAGt      >  2:301768/1‑149 (MQ=255)
                                                                                                                                         gTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGTAAGGGATATTCTGCTCGCTGACGCCGACTTCTACCGTAACGGTGTAGt      >  3:332210/1‑148 (MQ=255)
                                                                                                                                         gTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGTAAGGGATATTCTGCTCGCTGACGCCGACTTCTACCGTAACGGTGTAGTc     <  4:254492/149‑1 (MQ=255)
                                                                                                                                             ccgccgACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGTAAGGGATATTCTGCTCGCTGACGCCGACTTCTACCGTAACGGTGTAGTCATg  >  4:277226/1‑148 (MQ=255)
                                                                                                                                             ccgccgACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGTAAGGGATATTCTGCTCGCTGACGCCGACTTCTACCGTAACGGTGTAGTCATg  >  5:213801/1‑148 (MQ=255)
                                                                                                                                             ccgccgACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAAAGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGTAAGGGATATTCTGCTCGCTGACGCCGACTTCTACCGTAACGGTGTAGTCATg  >  2:367967/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
CTGTACATGGCGCCAGTCGCTGCCGGTGTAGTGCACCAGGCGGCGTAGCGAGAAGTCGACGCGCGCGGCATC‑GAACAACGCCAGTGGCATGATTTCAGCACGGTCCCAGTCGTACAGGCCATCAGCGATGTTGTCGGTGGCCGCCGACAGGTCGGTGCTGGGGAACACCCGCGCCAGCTGCGCAGCGGTAATGCCGCTACCGGCCAGCTCGTCGCCCTGGTCGACCACATAGGGGTAAGGGATATTCTGCTCGCTGACGCCGACTTCTACCGTAACGGTGTAGTCATG  >  NC_002947/5440226‑5440513

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