Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,782,535 C→A 8.7% Q514H (CAG→CAT) ‡ glnD ← uridylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,782,5350CA 8.7% 276.3 / 13.6 92Q514H (CAG→CAT) ‡glnDuridylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (41/43);  new base A (3/5);  total (44/48)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.09e-01

GTCGGACAAGTCCTTGCGCTGAGCGGTGGTCGACATCACCAGGTGGTTCTGCACCAGCCAGACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGAATTTCTCGGACACCGGCGTGTATTGCAGCTTGCGCAGGTGCTTGATGAGGTTGAG  >  NC_002947/1782388‑1782681
                                                                                                                                                   |                                                                                                                                                     
gTCGGACAAGTCCTTGCGCTGAGCGGTGGTCGACATCACCAGGTGGTTCTGCACGAGCCAGACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCt                                                                                                                                                      <  3:102061/149‑1 (MQ=255)
   ggACAAGTCCTTGCGCTGAGCGGTGGTCGACATCACCAGGTGGTTCTGCACCAGCCAGACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTgg                                                                                                                                                    <  3:261265/148‑1 (MQ=255)
       aaGTCCTTGCGCTGAGCGGTGGTCGACATCACCAGGTGGTTCTGCACCAGCCAGACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGAcc                                                                                                                                               >  4:86448/1‑149 (MQ=255)
         gTCCTTGCGCTGAGCGGTGGTCGACCTCACCAGGTGGTTCTGCACCAGCCAGACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTACTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCt                                                                                                                                            >  8:134558/1‑149 (MQ=255)
         gTCCTTCCGCTCAGCGGTGGTCCACATCACCAGGACCTTCTGCACCAGCCAGACGATCAGGCGGCTGTCACAGGCCGGCAACTTCAGCCGTTAGCAGAACTACTGACCGTCCACTCAACCGATTACCGACTGGTCGCCCTGGACAccc                                                                                                                                              >  8:139325/1‑148 (MQ=255)
           ccTTGCGCTGAGCGGTGGTCGACATCACCAGGTGGTTCTGCACCAGCCAGACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCAGGGCGAccact                                                                                                                                            <  4:3334/148‑4 (MQ=255)
            cTTGCGCTGAGCGGTGGTCGACATCACCAGGTTGTTCTGCACCAGCCAGACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCAGGGCGAccactcg                                                                                                                                          <  6:195269/149‑6 (MQ=255)
            cTTGCGCTGAGCGGTGGTCGACATCACCAGGTGGTTCTGCACCAGCCAGACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGc                                                                                                                                         <  3:86448/149‑1 (MQ=255)
             ttGCGCTGAGCGGTGGTCGACATCACCAGGTGGTTCTGCACCAGCCAGACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGcc                                                                                                                                        >  1:29246/1‑149 (MQ=255)
             ttGCGCTGAGCGGTGGTCGACATCACCAGGTGGTTCTGCACCAGCCAGACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGc                                                                                                                                         >  7:215480/1‑148 (MQ=255)
             ttGCGCTGAGCGGTGGTCGACATCACCAGGTGGTTCTGCACCAGCAAGACGATCAGGCGGCTGTCACAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGc                                                                                                                                         >  5:131316/1‑148 (MQ=255)
               gcgcTGAGCGGTGGTCGACATCACCAGGTGGTTCTGCACCAGCCAGACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGc                                                                                                                                         <  8:215480/146‑1 (MQ=255)
               gcgcTGAGCGGTGGTCGACATCACCAGGTGGTTCTGCACCAGCCAGACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGc                                                                                                                                         <  6:131316/146‑1 (MQ=255)
               gcgcTGAGCGGTGGTCGACATCACCAGGTGGTTCTGCACCAGCCAGACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGa                                                                                                                                      >  2:14337/1‑149 (MQ=255)
                cgcTGAGCGGTGCTCGACATCACCAGGGGGTTCTGCACCAGCCAGACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGAt                                                                                                                                     <  6:329323/149‑1 (MQ=255)
                     aGCGGTGGTCGACATCACCAGGTGGTTCTGCACCAGCCAGACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGt                                                                                                                                <  5:118251/149‑1 (MQ=255)
                            gtCGACATCACCAGGTGGTTCTGCACCAGCCAGACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACa                                                                                                                          <  1:161643/148‑1 (MQ=255)
                            gtCGACATCACCAGGTGGTTCTGCACCAGCCAGACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAg                                                                                                                         >  4:99634/1‑149 (MQ=255)
                                aCATCACCAGGTGGTTCTGCACCAGCCAGACGATCAGGCGGCTGTGCCAGGCCGGCAACTGGTGCCGTTCGCAGCACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCg                                                                                                                     >  5:222285/1‑149 (MQ=255)
                                 cATCACCAGGTGGTTCTGCACCAGCCAGACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCgg                                                                                                                    >  5:336019/1‑149 (MQ=255)
                                       cAGGTGGTTCTGCACCAGCCAGACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGc                                                                                                                                                   >  6:77942/1‑113 (MQ=255)
                                            ggtTCTGCACCAGCCAGACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGAt                                                                                                          >  8:12771/1‑148 (MQ=255)
                                                 tGCACCAGCCAGACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAg                                                                                                                >  5:146959/1‑137 (MQ=255)
                                                 tGCACCAGCCAGACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAg                                                                                                                <  6:146959/137‑1 (MQ=255)
                                                   cACCAGCCAGACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCgg                                                                                                 >  6:222730/1‑149 (MQ=255)
                                                    aCCAGCCAGACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCAGGGCGAccactcg                                                                                                                                          >  8:95888/1‑104 (MQ=255)
                                                    aCCAGCCAGACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCAGGGCGAccactcg                                                                                                                                          <  7:95888/109‑6 (MQ=255)
                                                      cagccagACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGt                                                                                                              >  2:62769/1‑134 (MQ=255)
                                                       agccagACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTTCCg                                                                                                                                       <  2:55324/108‑1 (MQ=255)
                                                       agccagACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCg                                                                                                                                       >  1:55324/1‑108 (MQ=255)
                                                         ccagACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGGCGTTAGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTg                                                                                            <  1:82672/149‑1 (MQ=255)
                                                         ccagACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGc                                                                                                      >  1:346445/1‑139 (MQ=255)
                                                         ccagACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGc                                                                                                      <  2:346445/139‑1 (MQ=255)
                                                          cagACAATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTg                                                                                                                               >  7:272652/1‑113 (MQ=255)
                                                          cagACAATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTg                                                                                                                               <  8:272652/113‑1 (MQ=255)
                                                           agACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTggg                                                                                          >  5:252159/1‑148 (MQ=255)
                                                           agACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTggg                                                                                          >  8:8106/1‑148 (MQ=255)
                                                             aCGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAg                                                                                       <  6:336019/149‑1 (MQ=255)
                                                               gATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGgc                                                                                     >  8:51940/1‑149 (MQ=255)
                                                                 tCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGgcgc                                                                                   <  3:99634/149‑1 (MQ=255)
                                                                   aGGCGGCTGTCCCAGGCCGGCAACTGGTGCCTTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGcc                                                                                  <  6:252159/148‑1 (MQ=255)
                                                                        gcTGTTCCAGGCCGGCAACTGGTGCCGGTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGa                                                                           <  6:222285/149‑1 (MQ=255)
                                                                          tGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCTGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGc                                                                         >  2:83551/1‑149 (MQ=255)
                                                                              ccAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGc                                                                     >  3:310622/1‑149 (MQ=255)
                                                                               cAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGGACCGATTTCCGAGTGGTCGCCCTGTCGACCCTTGCCTATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCt                                                                    <  2:238112/149‑1 (MQ=255)
                                                                               cAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCt                                                                    <  7:12771/149‑1 (MQ=255)
                                                                                        aaCTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCgg                                                           >  6:220874/1‑149 (MQ=255)
                                                                                         aCTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCgg                                                           >  3:324341/1‑148 (MQ=255)
                                                                                            ggTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGt                                                                                                              >  2:248144/1‑96 (MQ=255)
                                                                                            ggTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGAAGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGt                                                                                                              <  1:248144/96‑1 (MQ=255)
                                                                                             gTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGccc                                                                                 >  5:78954/1‑123 (MQ=255)
                                                                                             gTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGccc                                                                                 <  6:78954/123‑1 (MQ=255)
                                                                                             gTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCACGAGCTTGCTGGCCAGCGGGAAt                                                       <  7:51940/148‑1 (MQ=255)
                                                                                              tGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGAAttt                                                     <  5:222730/149‑1 (MQ=255)
                                                                                                 cGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTGGCCCTGGGGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGAATTtctc                                                  <  7:8106/149‑1 (MQ=255)
                                                                                                  gTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGAATTTCTCg                                                 >  1:321872/1‑149 (MQ=255)
                                                                                                  gTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGAATTTCTCg                                                 >  2:266360/1‑149 (MQ=255)
                                                                                                   ttCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCTGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGAATTTCTCg                                                 <  1:83551/149‑1 (MQ=255)
                                                                                                     cGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCt                                                                    >  6:210506/1‑127 (MQ=255)
                                                                                                     cGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCt                                                                    <  5:210506/127‑1 (MQ=255)
                                                                                                      gCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTATATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTg                                                                   >  7:235138/1‑127 (MQ=255)
                                                                                                      gCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGAATTTCTCg                                                 <  5:132074/145‑1 (MQ=255)
                                                                                                      gCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGAATTTCTCg                                                 >  6:132074/1‑145 (MQ=255)
                                                                                                      gCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAAATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTg                                                                   <  8:235138/127‑1 (MQ=255)
                                                                                                        aGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGACGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGAATTTCTCGGACAcc                                           >  5:334522/1‑149 (MQ=255)
                                                                                                           aCTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGAATTTCTCGGACACCGGc                                        <  2:29246/149‑1 (MQ=255)
                                                                                                             ttCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGAATTTCTCGGACACCGGCgt                                      >  4:45218/1‑149 (MQ=255)
                                                                                                                  cgcgTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGAATTTCTCGGACACCGGCgtgt                                    <  7:111090/146‑1 (MQ=255)
                                                                                                                   gcgTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGAATTtctc                                                  >  6:91552/1‑131 (MQ=255)
                                                                                                                   gcgTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGAATTtctc                                                  <  5:91552/131‑1 (MQ=255)
                                                                                                                   gcgTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGAATTTCTCGGACACCGGCGTGTATTg                                <  2:321872/149‑1 (MQ=255)
                                                                                                                     gTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGAATTTCTCGGACACTGGCGTGTATTGCa                              <  5:220874/149‑1 (MQ=255)
                                                                                                                     gTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGAATTTCTCGGACACCGGCGTGTATTGCa                              <  4:310622/149‑1 (MQ=255)
                                                                                                                     gTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGAATTTCTCGGACACCGGCGTGTATTGCa                              <  4:324341/149‑1 (MQ=255)
                                                                                                                      tCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGAATTTCTCGGACACCGGCGTGTATTGCAg                             >  7:40368/1‑149 (MQ=255)
                                                                                                                         aCTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGAATTTCTCGGACACCGGCGTGTATTGCAGCtt                          >  8:150167/1‑149 (MQ=255)
                                                                                                                             cACCGCTTTCCGAGTGGGCGCCCTGGCGTCCCTTGCCGTTGTCGTGGTAAAGCCCGGCCATGTAGATCAGCTGGGGTTTGGGCAGGCGCCCCATCGAGCTTGCTGGCCAGCGGGAATTTCTCGGACACCGGCGTGTATTGCAGCTTgc                        <  1:266360/148‑1 (MQ=255)
                                                                                                                             cACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGAATTTCTCGGACACCGGCGTGTATTGCAGCTTgcg                       >  7:281392/1‑148 (MQ=255)
                                                                                                                             cACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGAATTTCTCGGACACCGGCGTGTATTGCAGCTTgcg                       >  8:106813/1‑148 (MQ=255)
                                                                                                                               gctaTTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGTGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGAATTTCTCGTACACCGGCGTGTATTGCAGCTTGCGCAg                    <  7:139325/146‑1 (MQ=255)
                                                                                                                               ccGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGAATTTCTCGGACACCGGCGTGTATTGCAGCTTGCGCAg                    <  7:134558/149‑1 (MQ=255)
                                                                                                                                cGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGAATTTCTCGGACACCGGCGTGTATTGCAGCTTGCGCAgg                   <  8:281392/149‑1 (MQ=255)
                                                                                                                                       gcaagggTCGCCAGGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGAATTTCTCGGACACCGGCGTGTATTGCAGCTTGCGCAGGTGCt               <  5:282914/141‑1 (MQ=255)
                                                                                                                                       gcaagggTCGCCAGGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGAATTTCTCGGACACCGGCGTGTATTGCAGCTTGCGCAGGTGCt               >  6:282914/6‑146 (MQ=255)
                                                                                                                                       gcaagggTCGCCAGGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGAATTTCTCGGACACCGGCGTGTATTGCAGCTTGCGCAGGTGCt               <  5:278006/141‑1 (MQ=255)
                                                                                                                                       gcaagggTCGCCAGGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGAATTTCTCGGACACCGGCGTGTATTGCAGCTTGCGCAGGTGCt               >  6:278006/6‑146 (MQ=255)
                                                                                                                                        gAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGAATTTCTCGGACACCGGCGTGTATTGCAGCTTGCGCAGGTGCTtgat           <  6:334522/149‑1 (MQ=255)
                                                                                                                                          gTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGAATTTCTCGGACACCGGCGTGTATTGCAGCTTGCGCAGGTGCTtgatga         <  1:14337/149‑1 (MQ=255)
                                                                                                                                           tGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGAATTTCTCGGACACCGGCGTGTATTGCAGCTTGCGCAGGTGCTTGATGAg        >  2:175953/1‑149 (MQ=255)
                                                                                                                                             gTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGAATTTCTCGGACACCGGCGTGTATTGCAGCTTGCGCAGGTGCTTGATGAGGt      <  2:170493/149‑1 (MQ=255)
                                                                                                                                                gCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGAATTTCTCGGACACCGGCGTGTATTGCAGCTTGCGCAGGTGCTtgat           >  2:320781/1‑141 (MQ=255)
                                                                                                                                                gCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGAATTTCTCGGACACCGGCGTGTATTGCAGCTTGCGCAGGTGCTtgat           <  1:320781/141‑1 (MQ=255)
                                                                                                                                                 cccTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGAATTTCTCGGACACCGGCGTGTATTGCAGCTTGCGCAGGTGCTTGATGAGGTTGAg  <  1:175953/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GTCGGACAAGTCCTTGCGCTGAGCGGTGGTCGACATCACCAGGTGGTTCTGCACCAGCCAGACGATCAGGCGGCTGTCCCAGGCCGGCAACTGGTGCCGTTCGCAGAACTTCTGCGCGTCCACTGCACCGATTTCCGAGTGGTCGCCCTGGCGACCCTTGCCGATGTCGTGGTACAGCCCGGCCAGGTAGATCAGCTCGGGTTTGGGCAGGCGCCCCATGAGCTTGCTGGCCAGCGGGAATTTCTCGGACACCGGCGTGTATTGCAGCTTGCGCAGGTGCTTGATGAGGTTGAG  >  NC_002947/1782388‑1782681

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: