Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,125,857 G→T 5.1% G218W (GGG→TGG)  PP_5374 → choline/carnitine/betaine transporter family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,125,8570GT 5.1% 253.0 / 3.2 79G218W (GGG→TGG) PP_5374choline/carnitine/betaine transporter family protein
Reads supporting (aligned to +/- strand):  ref base G (45/30);  new base T (2/2);  total (47/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.42e-01

CCGCTGAGCATGCGCTCGATCCTTTACCCGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCGGAATGGAATATGCTGTTGTCGGT  >  NC_002947/6125710‑6126005
                                                                                                                                                   |                                                                                                                                                    
ccGCTGAGCATGCGCTCGATCCTTTACCCGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCg                                                                                                                                                      <  5:78614/148‑1 (MQ=255)
  gCTGAGCATGCGCTCGATCCTTTACCCGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGc                                                                                                                                                   >  5:410638/1‑149 (MQ=255)
      aGCATGCGCTCGATCCTTTACCCGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGa                                                                                                                                                <  6:410638/148‑1 (MQ=255)
      aGCATGCGCTCGATCCTTAACCCGCTGATTGGGGAGCGCATCTATGGGCCCATCCGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGGCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGa                                                                                                                                                <  2:127119/148‑1 (MQ=255)
          tGCGCTCGATCCTTTACCCGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAg                                                                                                                                           <  1:404337/149‑1 (MQ=255)
             gcTCGATCCTTTACCCGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCTCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGgtg                                                                                                                                        <  4:166772/149‑1 (MQ=255)
             gcTCGATCCTTTACCCGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGgtg                                                                                                                                        >  2:316242/1‑149 (MQ=255)
                 gATCCTTTACCCGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCTGGCTCAGCCCGGTGtt                                                                                                                                      <  6:92056/147‑1 (MQ=255)
                 gATCCTTTACCCGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGttt                                                                                                                                     <  6:378862/148‑1 (MQ=255)
                     cTTTACCCGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCTGGCTCAGCCCGGTGttcatttg                                                                                                                                <  6:165186/149‑7 (MQ=255)
                        tACCCGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCGGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGc                                                                                                                              <  5:304633/148‑1 (MQ=255)
                          cccGCTGATTGGTGAGCACATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGccc                                                                                                                            >  5:104606/1‑148 (MQ=255)
                           ccGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCAt                                                                                                                          >  5:381908/1‑149 (MQ=255)
                             gCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTCATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCAt                                                                                                                          >  2:231855/1‑147 (MQ=255)
                               tGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGc                                                                                                                      >  1:44713/1‑149 (MQ=255)
                                 aTTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGc                                                                                                                      <  2:93280/147‑1 (MQ=255)
                                 aTTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCt                                                                                                                    <  8:143504/149‑1 (MQ=255)
                                    ggTGAGCGCATCTATGGGCCCATCGGACATGTGGTTGACATGCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTgg                                                                                                                  <  4:149779/148‑1 (MQ=255)
                                     gTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCTGGCTCAGCCCGGTGttcat                                                                                                                                   >  1:172820/1‑127 (MQ=255)
                                     gTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTgggg                                                                                                                <  2:328878/149‑1 (MQ=255)
                                     gTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTgggg                                                                                                                >  5:232023/1‑149 (MQ=255)
                                     gTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGAGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTgggg                                                                                                                <  5:61016/149‑1 (MQ=255)
                                     gTGAGCGCATCTATGGGCCCATCGGACATCAGGTGGACAACCATGCCGTGGTGATTACCGCCTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATCAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCAGCAGCCTggtg                                                                                                                >  7:228562/1‑147 (MQ=255)
                                          cgcATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCa                                                                                                           <  2:112694/149‑1 (MQ=255)
                                           gcATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAg                                                                                                          <  6:232023/149‑1 (MQ=255)
                                            cATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTg                                                                                                             >  1:272177/1‑145 (MQ=255)
                                             aTCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCGCATCAGCCTGGGGGTg                                                                                                             >  4:222283/1‑144 (MQ=255)
                                                  tGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGTCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCAcc                                                                                                    <  7:258969/148‑1 (MQ=255)
                                                     gCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGa                                                                                                <  1:231855/149‑1 (MQ=255)
                                                     gCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGa                                                                                                >  2:205069/1‑149 (MQ=255)
                                                     gCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGa                                                                                                >  5:174216/1‑149 (MQ=255)
                                                       ccATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATc                                                                                              <  1:316242/149‑1 (MQ=255)
                                                        cATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCg                                                                                             <  1:229631/149‑1 (MQ=255)
                                                         aTCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGt                                                                                            >  7:231574/1‑149 (MQ=255)
                                                          tCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTg                                                                                           <  5:131177/149‑1 (MQ=255)
                                                              aCATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCa                                                                                       >  4:72853/1‑149 (MQ=255)
                                                                aTGTGGTGGACATCCTTGCCGTTGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCa                                                                                                                                            <  5:350366/94‑1 (MQ=255)
                                                                aTGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCa                                                                                                                                            >  6:350366/1‑94 (MQ=255)
                                                                  gtggtgGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTc                                                                                        >  7:114627/1‑144 (MQ=255)
                                                                   tggtggACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTAc                                                                                   >  2:410535/1‑148 (MQ=255)
                                                                   tggtggACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATAAGCCTGGGGGTGCAGATCACCTTGATCGTGCTAATTAcc                                                                                  >  5:94589/1‑149 (MQ=255)
                                                                     gtggACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCt                                                                                >  2:47202/1‑149 (MQ=255)
                                                                        gACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTgggg                                                                                                                >  2:94224/1‑114 (MQ=255)
                                                                              cTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTtacta                                                                        <  2:44713/148‑1 (MQ=255)
                                                                                     tggtgATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCg                                                                <  1:47202/149‑1 (MQ=255)
                                                                                     tggtgATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCg                                                                >  8:43334/1‑149 (MQ=255)
                                                                                      ggtgATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCgg                                                               >  1:302514/1‑149 (MQ=255)
                                                                                        tgATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCAt                                                                                      <  5:232569/124‑1 (MQ=255)
                                                                                        tgATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCAt                                                                                      >  6:232569/1‑124 (MQ=255)
                                                                                         gATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGa                                                            >  7:362196/1‑149 (MQ=255)
                                                                                          aTTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGa                                                            <  7:43334/148‑1 (MQ=255)
                                                                                          aTTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGAt                                                           >  1:287931/1‑149 (MQ=255)
                                                                                               cgcgTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTAcc                                                                      >  1:67686/1‑133 (MQ=255)
                                                                                               cgcgTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTAcc                                                                      <  2:67686/133‑1 (MQ=255)
                                                                                               cgcgTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACAGTATCGGTGATGGCCg                                                      >  8:295630/1‑149 (MQ=255)
                                                                                                     cGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTg                                                             >  5:11273/1‑136 (MQ=255)
                                                                                                     cGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTg                                                             <  6:11273/136‑1 (MQ=255)
                                                                                                         gtATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTcgcg                                            >  8:173074/1‑149 (MQ=255)
                                                                                                              ccAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCa                                       <  1:410535/149‑1 (MQ=255)
                                                                                                              ccAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCa                                       >  4:146828/1‑149 (MQ=255)
                                                                                                               cAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTc                                               >  5:238793/1‑140 (MQ=255)
                                                                                                               cAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTc                                               <  6:238793/140‑1 (MQ=255)
                                                                                                               cAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCAt                                      >  7:45714/1‑149 (MQ=255)
                                                                                                                       gggcctgggcGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGcc                              >  8:343334/1‑149 (MQ=255)
                                                                                                                          cctgggcGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCttt                           >  5:50790/1‑149 (MQ=255)
                                                                                                                          cctgggcGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCttt                           >  8:108494/1‑149 (MQ=255)
                                                                                                                           ctgggcGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTc                          <  6:104606/149‑1 (MQ=255)
                                                                                                                              ggcGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCgg                        >  4:87661/1‑148 (MQ=255)
                                                                                                                               gcGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCggaa                      >  2:364318/1‑149 (MQ=255)
                                                                                                                                     tgCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGa                                    >  1:198188/1‑129 (MQ=255)
                                                                                                                                     tgCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGa                                    <  2:198188/129‑1 (MQ=255)
                                                                                                                                      gCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCGGAATGGAata               >  6:7065/1‑149 (MQ=255)
                                                                                                                                        aaaTGAACACCGGGCTGAGCCAGGAGTTAGACCTCACCATCAAACTTGGGGTCCAGATCACCTTGATCGAGCTCATTACCATCGTTAATAACGTATCGGTGATTGCAGCCGTGTCCAGTCCCATGAAGCCCCTATCGGAATTCAATATg             >  8:171612/1‑149 (MQ=255)
                                                                                                                                         aaTGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCATATCACCTTGATCGAGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCGGAATGGAAGATGc            >  3:256276/1‑149 (MQ=255)
                                                                                                                                         aaTGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGACTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCAGTCGTGTCGCGTGGCATGAAGCGCCTTTCGGAATGTAATATGc            >  4:75763/1‑149 (MQ=255)
                                                                                                                                          aTGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCGGAATGGAATATGCt           <  3:72853/149‑1 (MQ=255)
                                                                                                                                          aTGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCGGAATGGAATATGCt           >  2:164740/1‑149 (MQ=255)
                                                                                                                                             aaCACCTGGCTCAGCCCGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCGGAATGGAATATGCtgtt        >  5:255546/1‑149 (MQ=255)
                                                                                                                                                  tgggCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCGGAATGGAATATGCTTTTGTCgg   >  2:427265/2‑149 (MQ=255)
                                                                                                                                                  cGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCGGAATGGAATATGCTGTTGTCgg   >  4:356087/1‑149 (MQ=255)
                                                                                                                                                  cGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCGGAATGGAATATGCTGTTGTCgg   >  8:53839/1‑149 (MQ=255)
                                                                                                                                                   gggCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCGGAATGGAATATGCTGTTGTCggt  <  8:37841/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCGCTGAGCATGCGCTCGATCCTTTACCCGCTGATTGGTGAGCGCATCTATGGGCCCATCGGACATGTGGTGGACATCCTTGCCGTGGTGATTACCGCGTTCGGTGTATCCCAGTCGCTGGGCCTGGGCGTGGTGCAAATGAACACCGGGCTGAGCCAGGTGTTTGACCTGCCCATCAGCCTGGGGGTGCAGATCACCTTGATCGTGCTCATTACCCTCGTTACTACCGTATCGGTGATGGCCGGCGTGTCGCGTGGCATGAAGCGCCTTTCGGAATGGAATATGCTGTTGTCGGT  >  NC_002947/6125710‑6126005

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: