Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 932,435 G→A 6.1% A4396V (GCC→GTC)  PP_0806 ← surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947932,4350GA 6.1% 251.3 / 3.8 82A4396V (GCC→GTC) PP_0806surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base G (42/35);  new base A (3/2);  total (45/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.02e-01

CAGGTTGTCCGGCGTGGCCGGTGCCAGCGCGCCCGGCGCCGTCACGTTTTGCGCAGGGGAGACATTCTGTGCGGCATCGGCCTGGACCACGCTCAGGCTTTCACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGC  >  NC_002947/932287‑932580
                                                                                                                                                    |                                                                                                                                                 
cAGGTTGTCCGGCGTGGCCGGTGCCAGCGCGCCCGGCGCCGTCACGTTTTGCGCAGGGGAGACATTCTGTGCGGCATCGGCCTGGACCACGCTCAGGCTTTCACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTcgg                                                                                                                                                   >  7:170783/1‑149 (MQ=255)
         cGGCGTGGCCGGTGCCAGCGCGCCCGGCGCCGTCACGTTTTGCGCAGGGGAGACATTCTGTGCGGCATCGGCCTGGACCACGCTCAGGCTTTCACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTg                                                                                                                                          >  8:95615/1‑149 (MQ=255)
           gCGTGGCCGGTGCCAGCGCGCCCGGCGCCGTCACGTTTTGCGCAGGGGAGACATTCTGTGCGGCATCGGCCTGGACCACGCTCAGGCTTTCACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTa                                                                                                                                        <  1:96067/149‑1 (MQ=255)
             gTGGCCGGTGCCAGCGCGCCCGGCGCCGGCACGTTTTGCGCAGGGGAGACATTCTGTGCGGCATCGGCCTGGACCACGCTCAGGCTTTCACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGc                                                                                                                                      >  4:299219/1‑149 (MQ=255)
                  cGGTGCCAGCGCGCCCGGCGCCGTCACGTTTTGCGCAGGGGAGACATTCTGTGCGGCATCGGCCTGGACCACGCTCAGGCTTTCACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGc                                                                                                                                 >  1:241140/1‑149 (MQ=255)
                   ggTGCCAGCGCGCCCGGCGCCGTCACGTTTTGCGCAGGGGAGACATTCTGTGCGGCATCGGCCTGGACCACGCTCAGGCTTTCACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGcc                                                                                                                                >  3:86389/1‑149 (MQ=255)
                        cAGCGCGCCCGGCGCCGTCACGTTTTGCGCAGGGGAGACATTCTGTGCGGCATCGGCCTGGACCACGCTCAGGCTTTCACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCa                                                                                                                           >  6:192493/1‑149 (MQ=255)
                         aGCGCGCCCGGCGCCGTCACGTTTTGCGCAGGGGAGACATTCTGTGCGGCATCGGCCTGGACCACGCTCAGGCTTTCACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCAc                                                                                                                          >  4:131576/1‑149 (MQ=255)
                         aGCGCGCCCGGCGCCGTCACGTTTTGCGCAGGGGAGACATTCTGTGCGGCATCGGCCTGGACCACGCTCAGGCTTTCACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCAc                                                                                                                          >  7:117539/1‑149 (MQ=255)
                          gcgcgcCCGGCGCCGTCACGTTTTGCGCAGGGGAGACATTCTGTGCGGCATCGGCCTGGACCACGCTCAGGCTTTCACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACg                                                                                                                         >  6:95884/1‑149 (MQ=255)
                             cgcCCGGCGCCGTCACGTTTTGCGCAGGGGAGACATTCTGAGCGGCATCGGCCTGGACCACGCTCAGGCTTTCACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTa                                                                                                                      <  4:47805/149‑1 (MQ=255)
                              gcCCGGCGCCGTCACGTTTTGCGCAGGGGAGACATTCTGTGCGGCATCGGCCTGGACCACGCTCAGGCTTTCACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTAc                                                                                                                     >  3:279821/1‑149 (MQ=255)
                               cccGGCGCCGTCACGTTTTGCGCAGGGGAGACATTCTGTGCGGCATCGGCCTGGACCACGCTCAGGCTTTCACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGcc                                                                                                                                >  3:55082/1‑137 (MQ=255)
                               cccGGCGCCGTCACGTTTTGCGCAGGGGAGACATTCTGTGCGGCATCGGCCTGGACCACGCTCAGGCTTTCACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGcc                                                                                                                                <  4:55082/137‑1 (MQ=255)
                                        gTCACGTTTTGCGCAGGGGAGACATTCTGTGCGGCATCGGCCTGGACCACGCTCAGGCTTTCACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCCGCACGGTACCGTGCGGTGc                                                                                                           >  8:129476/1‑149 (MQ=255)
                                           aCGTTTTGCGCAGGGGAGACATTCTGTGCGGCATCGGCCTGGACCACGCTCAGGCTTTCACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTc                                                                                                         <  3:112423/148‑1 (MQ=255)
                                                   cgcAGGGGAGACATTCTGTGCGGCATCGGCCTGGACCACGCTCAGGCTTTCACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTc                                                                                                 <  2:359649/148‑1 (MQ=255)
                                                      aGGGGAGACATTCTGTGCGGCATCGGCCTGGACCACGCTCAGGCTTTCACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACgg                                                                                             <  5:168614/149‑1 (MQ=255)
                                                          gagaCATTCTGTGCGGCATCGGCCTGGACCACGCTCAGGCTTTCACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTgc                                                                                         <  7:95615/149‑1 (MQ=255)
                                                           agaCATTCTGTGCGGCATCGGCCTGGACCACGCTCAGGCTTTCACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTgcg                                                                                        <  8:117539/149‑1 (MQ=255)
                                                            gaCATTCTGTGCGGCATCGGCCTGGACCACGCTCAGGCTTTCACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTgcgc                                                                                       <  5:95884/149‑1 (MQ=255)
                                                             aCATTCTGTGCGGCATCGGCCTGGACCACGCTCAGGCTTTCACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGcc                                                                                      >  7:94536/1‑149 (MQ=255)
                                                                     tgCGGCATCGGCCTGGACCACGCTCAGGCTTTCACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGc                                                                              >  3:143252/1‑149 (MQ=255)
                                                                          cATCGGCCTGGACCACGCTCAGGCTTTCACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCg                                                                         >  2:240399/1‑149 (MQ=255)
                                                                           aTCGGCCTGGACCACGCTCAGGCTTTCACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTc                                                                                >  7:23195/1‑141 (MQ=255)
                                                                           aTCGGCCTGGACCACGCTCAGGCTTTCACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTc                                                                                <  8:23195/141‑1 (MQ=255)
                                                                           aTCGGCCTGGACCACGCTCAGGCTTTCACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGc                                                                        >  8:22233/1‑149 (MQ=255)
                                                                               gCCTGGACCACGCTCAGGCTTTCACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTc                                                                                                                 >  3:237567/1‑104 (MQ=255)
                                                                               gCCTGGACCACGCTCAGGCTTTCACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTc                                                                                                                 <  4:237567/104‑1 (MQ=255)
                                                                                   ggACCACGCTCAGGCTTTCACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAgg                                                                >  5:215469/1‑149 (MQ=255)
                                                                                   ggACCACGCTCAGGCTTTCACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACAGTCCGGTGCTCGTACCGTAACGGTTGCGCCTGCTTCGCCATAGCCTGTCAgg                                                                >  3:144702/1‑149 (MQ=255)
                                                                                      ccACGCTCAGGCTTTCACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTg                                                             >  8:255248/1‑149 (MQ=255)
                                                                                            tCAGGCTTTCACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACg                                                                                              >  6:56748/1‑110 (MQ=255)
                                                                                            tCAGGCTTTCACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACg                                                                                              <  5:56748/110‑1 (MQ=255)
                                                                                            tCAGGCTTTCACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCa                                                       >  4:98803/1‑149 (MQ=255)
                                                                                            tCAGGCTTTCACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCa                                                       >  7:208737/1‑149 (MQ=255)
                                                                                              aGGCTTTCACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCAtt                                                     <  6:215469/149‑1 (MQ=255)
                                                                                              aGGCTTTCACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCAtt                                                     <  6:220486/149‑1 (MQ=255)
                                                                                              aGGCTTTCACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCAtt                                                     >  7:227421/1‑149 (MQ=255)
                                                                                              aGGCTTTCACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCAtt                                                     <  8:170783/149‑1 (MQ=255)
                                                                                              aGGCTTTCACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCAtt                                                     <  2:241140/149‑1 (MQ=255)
                                                                                                gCTTTCACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCGGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTgg                                                   >  4:201301/1‑149 (MQ=255)
                                                                                                gCTTTCACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTgcgc                                                                                       <  6:140416/113‑1 (MQ=255)
                                                                                                gCTTTCACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTgcgc                                                                                       >  5:140416/1‑113 (MQ=255)
                                                                                                gCTTTCACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTgg                                                   >  1:312053/1‑149 (MQ=255)
                                                                                                  tttCACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCg                                                 >  2:190933/1‑149 (MQ=255)
                                                                                                   ttCACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCg                                                 >  2:49459/1‑148 (MQ=255)
                                                                                                    tCACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCgg                                                <  4:86389/148‑1 (MQ=255)
                                                                                                     cACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGga                                               <  3:131576/148‑1 (MQ=255)
                                                                                                     cACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGga                                               <  7:255248/148‑1 (MQ=255)
                                                                                                     cACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGga                                               <  4:143252/148‑1 (MQ=255)
                                                                                                      aCCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGgaga                                             >  8:20856/1‑149 (MQ=255)
                                                                                                          gTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGgaga                                             <  8:139960/145‑1 (MQ=255)
                                                                                                          gTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGgaga                                             >  7:139960/1‑145 (MQ=255)
                                                                                                                gCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTg                                                    <  3:116538/132‑1 (MQ=255)
                                                                                                                gCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTg                                                    >  4:116538/1‑132 (MQ=255)
                                                                                                                 cTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTg                                  <  4:279821/149‑1 (MQ=255)
                                                                                                                    ccGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGc                                <  5:147081/148‑1 (MQ=255)
                                                                                                                    ccGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGc                                <  7:22233/148‑1 (MQ=255)
                                                                                                                     cGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTGCCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCa                              >  3:143232/1‑149 (MQ=255)
                                                                                                                       gTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTAGTACCGTCACCGTTGCGCCTGCCTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCAcg                            >  4:43007/1‑149 (MQ=255)
                                                                                                                        tGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGctcacg                           >  5:373207/1‑144 (MQ=255)
                                                                                                                         gACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCAcggc                          >  6:124019/1‑149 (MQ=255)
                                                                                                                           cAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCAcggcgg                        <  5:192493/149‑1 (MQ=255)
                                                                                                                              ggACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCAcggcggccg                     >  4:276689/1‑149 (MQ=255)
                                                                                                                              ggACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCAcggcggccg                     >  7:330695/1‑149 (MQ=255)
                                                                                                                                    gCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGgcggcg               >  3:19712/1‑149 (MQ=255)
                                                                                                                                     ctacaGTGGCCGCCGACGGTCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCgt              >  6:136341/5‑149 (MQ=255)
                                                                                                                                     cTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCgt              <  3:276689/149‑1 (MQ=255)
                                                                                                                                        caGTGGCCGCCGACGGTCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAg                                                                 <  4:81697/94‑1 (MQ=255)
                                                                                                                                        caGTGGCCGCCGACGGTCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAg                                                                 >  3:81697/2‑95 (MQ=255)
                                                                                                                                        caGTGGCCGCCGACGGTCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCgtggt           >  3:162353/2‑149 (MQ=255)
                                                                                                                                        aaGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGc                                         >  4:141792/1‑119 (MQ=255)
                                                                                                                                        aaGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGc                                         <  3:141792/119‑1 (MQ=255)
                                                                                                                                        aaGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCgtggt           >  7:275918/1‑149 (MQ=255)
                                                                                                                                        aaGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCgtg             <  5:124019/147‑1 (MQ=255)
                                                                                                                                        aaGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGCCGCCGGCCGCCgtgg            <  1:240399/148‑1 (MQ=255)
                                                                                                                                          gTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGgtgt         <  6:206900/149‑1 (MQ=255)
                                                                                                                                            gACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACa                                           <  4:53282/113‑1 (MQ=255)
                                                                                                                                            gACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACa                                           >  3:53282/1‑113 (MQ=255)
                                                                                                                                             gccgccgACGGTCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCg      <  4:162353/148‑1 (MQ=255)
                                                                                                                                                gTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTg   <  2:313684/149‑1 (MQ=255)
                                                                                                                                                 tCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGc  <  8:227421/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CAGGTTGTCCGGCGTGGCCGGTGCCAGCGCGCCCGGCGCCGTCACGTTTTGCGCAGGGGAGACATTCTGTGCGGCATCGGCCTGGACCACGCTCAGGCTTTCACCCGTGATCGCTGCCGGTGACAGGGACACGCTGAAGTGACCGTCGGCGGCCACTGTAGCGTTGCCCAGCACGGTACCGTCCGGTGCTCGTACCGTCACGGTTGCGCCTGCTTCGCCATCGCCTGTCAGGGTGGCGCCATTGGCGGAGACAGCCAGGTTGGCCACGGCGGCCGGCGGCGTGGTGTCCGCTGC  >  NC_002947/932287‑932580

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