Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,450,335 T→A 6.3% C137S (TGC→AGC)  nicA → nicotinate dehydrogenase subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,450,3350TA 6.3% 236.7 / 4.2 79C137S (TGC→AGC) nicAnicotinate dehydrogenase subunit A
Reads supporting (aligned to +/- strand):  ref base T (40/34);  new base A (2/3);  total (42/37)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.13e-01

CGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGATCGCATTCCCGAGCCCAGCGACGAACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCGGAGAAAGCCATGAACCACTCACAACAGGTGCCCAGCCGCGACCAACTGCTGGCCAAGACCGGCGTGTTGCTGATCGTCGACCAGATCACCCCGC  >  NC_002947/4450187‑4450481
                                                                                                                                                    |                                                                                                                                                  
cGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGATCGCATTCCCGAGCCCAGCGACGAACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCt                                                                                                                                                    >  3:251940/1‑149 (MQ=255)
 gTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGATCGCATTCCCGAGCCCAGCGACGAACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTg                                                                                                                                                   >  6:284197/1‑149 (MQ=255)
 gTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGATCGCATTCCCGAGCCCAGCGACGAACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTg                                                                                                                                                   >  5:182490/1‑149 (MQ=255)
  tGCAACAGGCGTTCATCGACTAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGATCGCATTCCCGAGCCCAGCGACGAACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGc                                                                                                                                                  <  5:393085/149‑1 (MQ=255)
    cAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGATCGCATTCCCGAGCCCAGCGACGAACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCg                                                                                                                                                 <  1:224604/148‑1 (MQ=255)
        aGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGATCGCATTCCCGAGCCCAGCGACGAACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCAGCGGCAca                                                                                                                                            <  2:300497/149‑2 (MQ=255)
         ggCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGATCGCATTCCCGAGCCCAGCGACGAACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCAGCGGCAcagg                                                                                                                                          <  8:368419/149‑4 (MQ=255)
            gTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGATCGCATTCCCGAGCCCAGCGACGAACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCAGCGGCAcagg                                                                                                                                          <  4:128700/147‑4 (MQ=255)
            gTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGATCGCATTCCCGAGCCCAGCGACGAACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCAGCGGCAcagg                                                                                                                                          >  3:128700/1‑144 (MQ=255)
                 tCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGATCGCATTCCCGAGCCCAGCGACGAACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGa                                                                                                                                    >  8:208158/1‑148 (MQ=255)
                   gacgaGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGATCGCATTCCCGAGCCCAGCGACGAACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGaaa                                                                                                                                  <  5:355521/148‑1 (MQ=255)
                    acgaGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGATCTCATTCCCGAGCCCAGCGACGAACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATc                                                                                                                                <  6:182490/149‑1 (MQ=255)
                    acgaGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGATCGCATTCCCGAGCCCAGCGACGAACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATc                                                                                                                                >  4:161249/1‑149 (MQ=255)
                    acgaGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGATCGCATTCCCGAGCCCAGCGACGAACAGATTCGCAACGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATc                                                                                                                                <  1:80912/149‑1 (MQ=255)
                     cgaGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGATCGCATTCCCGAGCCCAGCGACGAACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATcc                                                                                                                               >  8:44485/1‑149 (MQ=255)
                        gCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGATCGCATTCCCGAGCCCAGCGACGAACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTgc                                                                                                                            <  2:30157/149‑1 (MQ=255)
                           ggcggcCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGATCGCATTCCCGAGCCCAGCGACGAACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTgcgcg                                                                                                                         <  3:302829/149‑1 (MQ=255)
                                           tACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGATCGCATTCCCGAGCCCAGCGACGAACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGc                                                                                                         <  7:208158/149‑1 (MQ=255)
                                            aCTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGATCGCATTCCCGAGCCCAGCGACGAACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGc                                                                                                         <  6:330761/148‑1 (MQ=255)
                                             cTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGATCGCATTCCCGAGCCCAGCGACGAACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCg                                                                                                       <  6:259204/149‑1 (MQ=255)
                                              tGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGATCGCTTTCCCGAGCCCAGCGACGAACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGAGCTGcc                                                                                                        >  1:395697/1‑147 (MQ=255)
                                               gCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGATCGCATTCCCGAGCCCAGCGACGAACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGcc                                                                                                        >  3:44183/1‑146 (MQ=255)
                                                cATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGATCGCATTCCCGAGCCCAGCGACGATCAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTCCGCCGCGCTGCCgaga                                                                                                    >  7:114012/1‑149 (MQ=255)
                                                 gtGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGATCGCATTCCCGAGCCCAGCGACGAACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGAc                                                                                                   >  6:58685/2‑149 (MQ=255)
                                                 aTGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGATCGCATTCCCGAGCCCAGCGACGAACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGAc                                                                                                   >  4:148886/1‑149 (MQ=255)
                                                    aaCGGCATGATCATGACCGCCAAGGCCCTGCTCGATCGCATTCCCGAGCCCAGCGACGAACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGAccc                                                                                                 <  5:284197/148‑1 (MQ=255)
                                                     aCGGCATGATCATGACCGCCAAGGCCCTGCTCGATCGCATTCCCGAGCCCAGCGACGAACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCg                                                                                                >  2:193228/1‑148 (MQ=255)
                                                      cGGCATGATCATGACCGCCAAGGCCCTGCTCGATCGCATTCCCGAGCCCAGCGACGAACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCGga                                                                                              <  6:104606/149‑1 (MQ=255)
                                                      cGGCATGATCATGACCGCCAAGGCCCTGCTCGATCGCATTCCCGAGCCCAGCGACGAACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCGga                                                                                              <  5:58685/149‑1 (MQ=255)
                                                       ggCATGATCATGACCGCCAAGGCCCTGCTCGATCGCATTCCCGAGCCCAGCGACGAACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCg                                                                                                >  4:342431/1‑146 (MQ=255)
                                                          aTGATCATGACCGCCAAGGCCCTGCTCGATCGCATTCCCGAGCCCAGCGACGAACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCGgaga                                                                                            >  2:357493/1‑147 (MQ=255)
                                                          aTGATCATGACCGCCAAGGCCCTGCTCGATCGCATTCCCGAGCCCAGCGACGAACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCGGAGaaa                                                                                          <  7:44485/149‑1 (MQ=255)
                                                              tCATGACCGCCAAGGCCCTGCTCGATCGCATTCCCGAGCCCAGCGACGAACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCGGAGAAAg                                                                                         >  1:338290/1‑146 (MQ=255)
                                                              tCATGACCGCCAAGGCCCTGCTCGATCGCATTCCCGAGCCCAGCGACGAACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCGGAGAAAGCCa                                                                                      <  3:161249/149‑1 (MQ=255)
                                                                 tGACCGCCAAGGCCCTGCTCGATCGCATTCCCGAGCCCAGCGACGAACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCGGAGAAAGCCATg                                                                                    >  7:218955/1‑148 (MQ=255)
                                                                   aCCGCCAAGGCCCTGCTCGATCGCATTCCCGAGCCCAGCGACGAACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCAGCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCGGAGAAAGCCATGAAc                                                                                 <  8:101941/149‑1 (MQ=255)
                                                                   aCCGCCAAGGCCCTGCTCGATCGCATTCCCGAGCCCAGCGACGAACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCGGAGAAAGCCATGAAc                                                                                 <  8:218955/149‑1 (MQ=255)
                                                                     cGCCAAGGCCCTGCTCGATCGCATTCCCGAGCCCAGCGACGAACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCGGAGAAAGCCATGAACCa                                                                               >  7:219046/1‑149 (MQ=255)
                                                                      gCCAAGGCCCTGCTCGATCGCATTCCCGAGCCCAGCGACGAACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCGGAGAAAGCCATGAACCAc                                                                              <  3:182076/149‑1 (MQ=255)
                                                                       ccAAGGCCCTGCTCGATCGCATTCCCGAGCCCAGCGACGAACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCGGAGAAAGCCATGAACCACt                                                                             >  3:365790/1‑149 (MQ=255)
                                                                        cAAGGCCCTGCTCGATCGCATTCCCGAGCCCAGCGACGAACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCGGAGAAAGCCATGAACCACt                                                                             >  3:117205/1‑148 (MQ=255)
                                                                            gCCCTGCTCGATCGCATTCCCGAGCCCAGCGACGAACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCGGAGAAAGCCATGAACCACTCacaa                                                                        <  3:342431/149‑1 (MQ=255)
                                                                              ccTGCTCGATCGCATTCCCGAGCCCAGCGACGAACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCGGAGAAAGCCATGAACCACTCacaaca                                                                      <  4:251940/149‑1 (MQ=255)
                                                                                     gATCGCATTCCCGAGCCCAGCGACGAACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGATGCGGCACCCACGTCGAAATCCTGCGCGCCGAGCGCCGCGCTGCCGAGACCCGGAGAAAGCCATGAACCACTCACAACAGGTGcc                                                                >  7:75803/1‑148 (MQ=255)
                                                                                            ttCCCGAGCCCAGCGACGAACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCt                                                                                                                              <  2:346261/79‑1 (MQ=255)
                                                                                            ttCCCGAGCCCAGCGACGAACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCt                                                                                                                              >  1:346261/1‑79 (MQ=255)
                                                                                            ttCCCGAGCCCAGCGACGAACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCGGAGAAAGCCATGAACCACTCACAACAGGTGCCCAGCcgc                                                         >  4:240816/1‑148 (MQ=255)
                                                                                                cGAGCCCAGCGACGAACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCGGAGAAAGCCATGAACCACTCACAACAGGTGCCCAGCCGCGAcc                                                     >  2:90201/1‑148 (MQ=255)
                                                                                                    cccAGCGACGAACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCGGAGAAAGCCATGAACCACTCACAACAGGTGCCCAGCCGCGACCAActg                                                >  6:60589/1‑149 (MQ=255)
                                                                                                       aGCGACGAACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCGGAGAAAGCCATGAACCACTCACAACAGGTGCCCAGCCGCGACCAActgctg                                             <  1:357493/149‑1 (MQ=255)
                                                                                                        gCGACGAACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCGGAGAAAGCCATGAACCACTCACAACAGGTGCCCAGCCGCGACCAACTGCTgg                                            >  5:360712/1‑149 (MQ=255)
                                                                                                         cgacgaACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCGGAGAAAGCCATGAACCACTCACAACAGGTGCCCAGCCGCGACCAACTGCTGGc                                           <  8:219046/149‑1 (MQ=255)
                                                                                                           acgaACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCGga                                                                                              >  5:102607/1‑96 (MQ=255)
                                                                                                           acgaACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCGga                                                                                              <  6:102607/96‑1 (MQ=255)
                                                                                                            cgaACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCGGAGAAAGCCATGAACCACTCACAACAGGTGCCCAGCCGCGACCAACTGCTGGCCaa                                        <  3:148886/149‑1 (MQ=255)
                                                                                                             gaACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCGGAGAAAGCCATGAACCACTCACAACAGGTGCCCAGCCGCGACCAACTGCTGGCCAAg                                       >  4:276147/1‑149 (MQ=255)
                                                                                                             gaACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCGGAGAAAGCCATGAACCACTCACAACAGGTGCCCAGCCGCGACCAACTGCTGGCCAAg                                       >  1:151569/1‑149 (MQ=255)
                                                                                                              aaCAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCGGAGAAAGCCATGAACCACTCACAACAGGTGCCCAGCCGCGACCAACTGCTGGCCAAGa                                      >  4:320707/1‑149 (MQ=255)
                                                                                                               aCAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCGGAGAAAGCCATGAACCACTCACAACAGGTGCCCAGCCGCGACCAACTGCTGGCCAAGa                                      <  2:395697/148‑1 (MQ=255)
                                                                                                               aCAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCGGAGAAAGCCATGAACCACTCACAACAGGTGCCCAGCCGCGACCAACTGCTGGCCAAGAc                                     >  7:330806/1‑149 (MQ=255)
                                                                                                                cAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCGGAGAAAGCCATGAACCACTCACAACAGGTGCCCAGCCGCGACCAACTGCTGGCCAAGAcc                                    <  2:338290/149‑1 (MQ=255)
                                                                                                                    ttCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCGGAGAAAGCCATGAACCACTCACAACAGGTGCCCAGCCGCGACCAACTGCTGGCCAAGACCGGc                                 <  5:60589/148‑1 (MQ=255)
                                                                                                                       gCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCGGAGAAAGCCATGAACCACTCACAACAGGTGCCCAGCCGCGACCAACTGCTGGCCAAGACCGGCgtgt                             >  4:359971/1‑149 (MQ=255)
                                                                                                                       gCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCGGAGAAAGCCATGAACCACTCACAACAGGTGCCCAGCCGCGACCAACTGCTGGCCAAGACCGGCgtgt                             >  5:137901/1‑149 (MQ=255)
                                                                                                                        cAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCGGAGAAAGCCATGAACCACTCACAACAGGTGCCCAGCCGCGACCAACTGCTGGCCAAGACCGGCGTGtt                            >  2:206428/1‑149 (MQ=255)
                                                                                                                        cAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCGGAGAAAGCCATGAACCACTCACAACAGGTGCCCAGCCGCGACCAACTGCTGGCCAAGACCGGCGTGtt                            <  1:193228/149‑1 (MQ=255)
                                                                                                                             aGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCGGAGAAAGCCATGAACCACTCACAACAGGTGCCCAGCCGCGACCAACTGCTGGCCAAGACCGGCGTGTTGCTGa                       <  7:7612/149‑1 (MQ=255)
                                                                                                                              gCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCGGAGAAAGCCATGAACCACTCACAACAGGTGCCCAGCCGCGACCAACTGCTGGCCAAGACCGGCGTGTTGCTGAt                      >  3:125825/1‑149 (MQ=255)
                                                                                                                                  tCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCGGAGAAAGCCATGAACCACTCACAACAGGTGCCCAGCCGCGACCACCTGCTGGCCAAGACCGGCGTGTTGCTGAtcgtc                  >  2:190801/1‑149 (MQ=255)
                                                                                                                                  tCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCGGAGAAAGCCATGAACCACTCACAACAGGTGCCCAGCCGCGACCAACTGCTGGCCAAGACCGGCGTGTTGCTGAtcgtc                  >  5:325393/1‑149 (MQ=255)
                                                                                                                                  tCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCGGAGAAAGCCATGAACCACTCACAACAGGTGCCCAGCCGCGACCAACTGCTGGCCAAGACCGGCGTGTTGCTGAtcgtc                  <  3:359971/149‑1 (MQ=255)
                                                                                                                                       acgtgggtgCCGCAGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCGGAGAAAGCCATGAACCACTCACAACAGGTGCACAGCCGCGACCAACTGCTGGCCAAGACCGGCGTGTTGCTGATCGTCGAc               >  1:14208/7‑147 (MQ=255)
                                                                                                                                         aCCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCGGAGAAAGCCATGAACCACTCACAACAGGTGCCCAGCCGCGACCAACTGCTGGCCAAGACCGGCGTGTTGCTGATCGTCGACCAg            >  3:352919/1‑148 (MQ=255)
                                                                                                                                          ccTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCGGAGAAAGCCATGAACCACTCACAACAGGTGCCCAGCCGCGACCAACTGCTGGCCAAGACCGGCGTGTTGCTGATCGTCGAc               >  1:64557/1‑144 (MQ=255)
                                                                                                                                                ccGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCGGAGAAAGCCATGAACCACTCACAACAGGTGCCCAGCCGCGACCAACTGCTGGCCAAGACCGGCGTGTTGCTGATCGTCGACCAGATCAcccc    >  4:95666/1‑149 (MQ=255)
                                                                                                                                                ccGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCGGAGAAAGCCATGAACCACTCACAACAGGTGCCCAGCCGCGACCAACTGCTGGCCAAGACCGGCGTGTTGCTGATCGTCGACCAGATCAcccc    <  3:240816/149‑1 (MQ=255)
                                                                                                                                                ccGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCGGAGAAAGCCATGAACCACTCACAACAGGTGCCCAGCCGCGACCAACTGCTGGCCAAGACCGGCGTGTTGCTGATCGTCGACCAGATCAcccc    >  2:226214/1‑149 (MQ=255)
                                                                                                                                                 cGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCGGAGAAAGCCATGAACCACTCACAACAGGTGCCCAGCCGCGACCAACTGCTGGCCAAGACCGGCGTGTTGCTGATCGTCGACCAGATCACCCCg   <  2:64557/149‑1 (MQ=255)
                                                                                                                                                  gCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCGGAGAAAGCCATGAACCACTCACAACAGGTGCCCAGCCGCGACCAACTGCTGGCCAAGACCGGCGTGTTGCTGATCGTCGACCAGATCACCCCGc  <  3:95666/149‑1 (MQ=255)
                                                                                                                                                  gCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCCGAGAAAGCCATGAACCACTCACAACAGGTGCCCAGCCGCGACCAACTGCTGGCCAAGACCGGCGTGTTGCTGATCGTCGACCAGATCACCCCGc  >  1:300638/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGATCGCATTCCCGAGCCCAGCGACGAACAGATTCGCAATGAGCTGTCCGCCAACCTGTGCCGCTGCGGCACCCACGTCGAAATCCTGCGCGCCGTGCGCCGCGCTGCCGAGACCCGGAGAAAGCCATGAACCACTCACAACAGGTGCCCAGCCGCGACCAACTGCTGGCCAAGACCGGCGTGTTGCTGATCGTCGACCAGATCACCCCGC  >  NC_002947/4450187‑4450481

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: