Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,366,231 C→G 5.7% E157Q (GAA→CAA) ‡ PP_1191 ← S4 domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,366,2310CG 5.7% 219.6 / 4.1 70E157Q (GAA→CAA) ‡PP_1191S4 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (34/32);  new base G (2/2);  total (36/34)
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.00e-01

AGCTCCAGGCGCAGGTAGGCGCATAGCTCGGAAATCTGTCCTGGTTGTGGGTTTTTCAACACCATGCGCAGGTGTGTCTCGTTCTGCCAGCTGGCCTTGGCCTTGGGCAACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGGGTGGTTTCGCCGCGCACTTCAATGATGAACTCCTGCTCCAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCA  >  NC_002947/1366083‑1366378
                                                                                                                                                    |                                                                                                                                                    
aGCTCCAGGCGCAGGTAGGCGCATAGCTCGGAAATCTGTCCTGGTTGTGGGTTTTTCAACACCATGCGCAGGTGTGTCTCGTTCTGCCAGCTGGCCTTGGCCTTGGGCAACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTc                                                                                                                                                      <  3:210418/149‑1 (MQ=255)
      aGGCGCAGGTAGGCGCATAGCTCGGAAATCTGTCCTGGTTGTGGGTTTTTCAACACCATGCGCAGGTGTGTCTCGTTCTGCCAGCTGGCCTTGGCCTTGGGCAACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAgcgc                                                                                                                                                >  4:56469/1‑149 (MQ=255)
         cgcAGGTAGGCGCATAGCTCGGAAATCTGTCCTGGTTGTGGGTTTTTCAACACCATGCGCAGGTGTGTCTCGTTCTGCCAGCTGGCCTTGGCCTTGGGCAACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTg                                                                                                                                             <  2:264172/149‑1 (MQ=255)
           cAGGTAGGCGCATAGCTCGGAAATCTGTCCTGGTTGTGGGTTTTTCAACACCATGCGCAGGTGTGTCTCGTTCTGCCAGCTGGCCTTGGCCTTGGGCAACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTggg                                                                                                                                           <  3:56469/149‑1 (MQ=255)
           cAGGTAGGCGCATAGCTCGGAAATCTGTCCTGGTTGTGGGTTTTTCAACACCATGCGCAGGTGTGTCTCGTTCTGCCAGCTGGCCTTGGCCTTGGGCAACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGCTGGAACGCCTggc                                                                                                                                           <  3:5242/149‑2 (MQ=255)
            aGGTAGGCGCATAGCTCGGAAATCTGTCCTGGTTGTGGGTTTTTCAACACCATGCGCAGGTGTGTCTCGTTCTGCCAGCTGGCCTTGGCCTTGGGCAACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTgggg                                                                                                                                          >  1:263106/1‑149 (MQ=255)
            aGGTAGGCGCATAGCTCGGAAATCTGTCCTGGTTGTGGGTTTTTCAACACCATGCGCAGGTGTGTCTCGTTCTGCCAGCTGGCCTTGGCCTTGGGCAACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTgggg                                                                                                                                          <  6:325325/149‑1 (MQ=255)
              gTAGGCGCATAGCTCGGAAATCTGTCCTGGTTGTGGGTTTTTCAACACCATGCGCAGGTGTGTCTCGTTCTGCCAGCTGGCCTTGGCCTTGGGCAACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTgggggg                                                                                                                                        <  4:78657/149‑1 (MQ=255)
                aGGCGCATAGCTCGGAAATCTGTCCTGGTTGTGGGTTTTTCAACACCATGCGCAGGTGTGTCTCGTTCTGCCAGCTGGCCTTGGCCTTGGGCAACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGggtg                                                                                                                                     >  2:213136/1‑149 (MQ=255)
                aGGCGCATAGCTCGGAAATCTGTCCTGGTTGTGGGTTTTTCAACACCATGCGCAGGTGTGTCTCGTTCTGCCAGCTGGCCTTGGCCTTGGGCAACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGggtg                                                                                                                                     >  7:330131/1‑149 (MQ=255)
                 ggCGCATAGCTCGGAAATCTGTCCTGGTTGTGGGTTTTTCAACACCATGCGCAGGTGTGTCTCGTTCTGCCAGCTGGCCTTGGCCTTGGGCAACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGggtg                                                                                                                                     <  7:94528/148‑1 (MQ=255)
                 ggCGCATAGCTCGGAAATCTGTCCTGGTTGTGGGTTTTTCAACACCATGCGCAGGTGTGTCTCGTTCTGCCAGCTGGCCTTGGCCTTGGGCAACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGggtg                                                                                                                                     <  2:263106/148‑1 (MQ=255)
                    gcATAGCTCGGAAATCTGTCCTGGTTGTGGGTTTTTCAACACCATGCGCAGGTGTGTCTCGTTCTGCCAGCTGGCCTTGGCCTTGGGCAACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGggtggt                                                                                                                                   >  4:260306/1‑147 (MQ=255)
                    gcATAGCTCGGAAATCTGTCCTGGTTGTGGGTTTTTCAACACCATGCGCAGGTGTGTCTCGTTCTGCCAGCTGGCCTTGGCCTTGGGCAACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGggtggt                                                                                                                                   <  3:260306/147‑1 (MQ=255)
                       tAGCTCGGAAATCTGTCCTGGTTGTGGGTTTTTCAACACCATGCGCAGGTGTGTCTCGTTCTGCCAGCTGGCCTTGGCCTTGGGCAACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGGGTGGTTTcgc                                                                                                                              <  2:343824/149‑1 (MQ=255)
                         gCTCGGAAATCTGTCCTGGTTGTGGGTTTTTCAACACCATGCGCAGGTGTGTCTCGTTCTGCCAGCTGGCCTTGGCCTTGGGCAACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGGGTGGTTTcgcc                                                                                                                             <  1:213136/148‑1 (MQ=255)
                             ggAAATCTGTCCTGGTTGTGGGTTTTTCAACACCATGCGCAGGTGTGTCTCGTTCTGCCAGCTGGCCTTGGCCTTGGGCAACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGGGTGGTTTCGCCgcgc                                                                                                                         <  8:330131/148‑1 (MQ=255)
                              gAAATCTGTCCTGGTTGTGGGTTTTTCAACACCATGCGCAGGTGTGTCTCGTTCTGCCAGCTGGCCTTGGCCTTGGGCAACTCGCGGTCATTGCGCATGGCCCCACGCCCCAGGCGCTGGAACGCCTggc                                                                                                                                           >  6:28597/1‑129 (MQ=255)
                              gAAATCTGTCCTGGTTGTGGGTTTTTCAACACCATGCGCAGGTGTGTCTCGTTCTGCCAGCTGGCCTTGGCCTTGGGCAACTCGCGGTCATTGCGCATGGCCCCACGCCCCAGGCGCTGGAACGCCTggc                                                                                                                                           <  5:28597/130‑2 (MQ=255)
                                    tGTCCTGGTTGTGGGTTTTTCAACACCATGCGCAGGTGTGTCTCGTTCTGCCAGCTGGCCTTGGCCTTGGGCAACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGGGTGGTTTCGCCGCGCACTTCAat                                                                                                                 >  6:163933/1‑149 (MQ=255)
                                      tCCTGGTTGTGGGTTTTTCAACACCATGCGCAGGTGTGTCTCGTTCTGCCAGCTGGCCTTGGCCTTGGGCAACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGGGTGGTTTCGCCGCGCACTTCaa                                                                                                                  <  6:309321/146‑1 (MQ=255)
                                      tCCTGGTTGTGGGTTTTTCAACACCATGCGCAGGTGTGTCTCGTTCTGCCAGCTGGCCTTGGCCTTGGGCAACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGGGTGGTTTCGCCGCGCACTTCAatga                                                                                                               >  8:42557/1‑149 (MQ=255)
                                        cTGGTTGTGGGTTTTTCAACACCATGCGCAGGTGTGTCTCGTTCTGCCAGCTGGCCTTGGCCTTGGGCAACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGGGTGGTTTCGCCGCGCACTTCAatgatg                                                                                                             <  1:303005/149‑1 (MQ=255)
                                             tgtgGGTTTTTCAACACCATGCGCAGGTGTGTCTCGTTCTGCCAGCTGGCCTTGGCCTTGGGCAACTCGCGGTCATTGCGCATGTCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGGGTGGTTTCGCCGCGCACTTCAatg                                                                                                                >  4:87617/1‑141 (MQ=255)
                                             tgtgGGTTTTTCAACACCATGCGCAGGTGTGTCTCGTTCTGCCAGCTGGCCTTGGCCTTGGGCAACTCGCGGTCATTGCGCATGTCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGGGTGGTTTCGCCGCGCACTTCAatg                                                                                                                <  3:87617/141‑1 (MQ=255)
                                               tgGGTTTTTCAACACCATGCGCAGGTGTGTCTCGTTCTGCCAGCTGGCCTTGGCCTTGGGCAACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGGGTGGTTTCGCCGCGCACTTCAATGATGAACTCCt                                                                                                      >  5:208765/1‑149 (MQ=255)
                                                 ggTTTTTCAACACCATGCGCAGGTGTGTCTCGTTCTGCCAGCTGGCCTTGGCCTTGGGCAACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGGGTGGTTTCGCCGCGCACTTCAATGATGAACTCCTGc                                                                                                    <  5:163933/149‑1 (MQ=255)
                                                  gTTTTTCAACACCATGCGCAGGTGTGTCTCGTTCTGCCAGCTGGCCTTGGCCTTGGGCAACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTcc                                                                                                                                                     >  5:96027/1‑100 (MQ=255)
                                                  gTTTTTCAACACCATGCGCAGGTGTGTCTCGTTCTGCCAGCTGGCCTTGGCCTTGGGCAACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTcc                                                                                                                                                     <  6:96027/100‑1 (MQ=255)
                                                    ttttCAACACCATGCGCAGGTGTGTCTCGTTCTGCCAGCTGGCCTTGGCCTTGGGCAACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGGGTGGTTTCGCCGCGCACTTCAATGATGAACTCCTGCTc                                                                                                  >  8:140015/1‑148 (MQ=255)
                                                         aaCACCATGCGCAGGTGTGTCTCGTTCTGCCAGCTGGCCTTGGCCTTGGGCGACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGGGTGGTTTCGCCGCGCACTTCAATGATGAACTCCTGCTCCAGGc                                                                                             >  6:87375/1‑148 (MQ=255)
                                                         aaCACCATGCGCAGGTGTGTCTCGTTCTGCCAGCTGGCCTTGGCCTTGGGCAACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGGGTGGTTTCGCCGCGCACTTCAATGATGAAc                                                                                                          <  6:301179/135‑1 (MQ=255)
                                                         aaCACCATGCGCAGGTGTGTCTCGTTCTGCCAGCTGGCCTTGGCCTTGGGCAACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGGGTGGTTTCGCCGCGCACTTCAATGATGAAc                                                                                                          >  5:301179/1‑135 (MQ=255)
                                                          acacCATGCGCAGGTGTGTCTCGTTCTGCCAGCTGGCCTTGGCCTTGGGCAACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGGGTGGTTTCGCCGCGCACTTCAATGATGAACTCCTGCTCCAGGCGa                                                                                           >  4:330311/1‑149 (MQ=255)
                                                            acCATGCGCAGGTGTGTCTCGTTCTGCCAGCTGGCCTTGGCCTTGGGCAACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCt                                                                                                                                              <  4:337415/97‑1 (MQ=255)
                                                            acCATGCGCAGGTGTGTCTCGTTCTGCCAGCTGGCCTTGGCCTTGGGCAACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCt                                                                                                                                              >  3:337415/1‑97 (MQ=255)
                                                             ccATGCGCAGGTGTGTCTCGTTCTGCCAGCTGGCCTTGGCCTTGGGCAACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGGGTGGTTTCGCCGCGCACTTCAat                                                                                                                 <  8:98057/124‑1 (MQ=255)
                                                             ccATGCGCAGGTGTGTCTCGTTCTGCCAGCTGGCCTTGGCCTTGGGCAACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGGGTGGTTTCGCCGCGCACTTCAat                                                                                                                 >  7:98057/1‑124 (MQ=255)
                                                              cATGCGCAGGTGTGTCTCGTTCTGCCAGCTGGCCTTGGCCTTGGGCAACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTggggg                                                                                                                                         >  2:183873/1‑100 (MQ=255)
                                                                  cgcAGGTGTGTCTCGTTCTGCCAGCTGGCCTTGGCCTTGGGCAACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGGGTGGTTTCGCCGCGCACTTCAATGATGAACTCCTGCTCCAGGCGACGCAGGt                                                                                    <  3:330311/148‑1 (MQ=255)
                                                                   gcAGGTGTGTCTCGTTCTGCCAGCTGGCCTTGGCCTTGGGCAACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGGGTGGTTTCGCCGCGCACTTCAATGATGAACTCCTGCTCCAGGCGACGCAGGTCg                                                                                  >  4:331094/1‑149 (MQ=255)
                                                                      ggTGTGTCTCGTTCTGCCAGCTGGCCTTGGCCTTGGGCAACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGGGTGGTTTCGCCGCGCACTTCAATGATGAACTCCTGCTCCAGGCGACGCAGGTCGGcc                                                                               <  7:140015/149‑1 (MQ=255)
                                                                                gTTCTGCCAGCTGGCCTTGGCCTTGGGCAACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGGGTGGTTTCGCCGCGCACTTCAATGATGAACTCCTGCTCCAGGCGACGCAGGTCGGCCTCGATCTTGc                                                                     <  3:331094/149‑1 (MQ=255)
                                                                                    tGCCAGCTGGCCTTGGCCTTGGGCAACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGGGTGGTTTCGCCGCGCACTTCAATGATGAACTCCTGCTCCAGGCGACGCAGGTCGGCCTCGATCTTGCGGGa                                                                 >  7:174200/1‑149 (MQ=255)
                                                                                     gCCAGCTGGCCTTGGCCTTGGGCAACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGGGTGGTTTCGCCGCGCACTTCAATGATGAACTCCTGCTCCAGGCGACGCAGGTCGGCCTCGATCTTGCGGGa                                                                 >  5:239319/1‑148 (MQ=255)
                                                                                           tggccttggcctTGGGCAACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGGGTGGTTTCGCCGCGCACTTCAATGATGAACTCCTGCTCCAGGCGa                                                                                           <  8:288970/116‑1 (MQ=255)
                                                                                           tggccttggcctTGGGCAACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGGGTGGTTTCGCCGCGCACTTCAATGATGAACTCCTGCTCCAGGCGa                                                                                           >  7:288970/1‑116 (MQ=255)
                                                                                           tggccttggcctTGGGCAACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGGGTGGTTTCGCCGCGCACTTCAATGATGAACTCCTGCTCCAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGc                                                          <  8:174200/149‑1 (MQ=255)
                                                                                              ccttggccttggGCAACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGGGTGGTTTCGCCGCGCACTTCAATGATGAACTCCTGCTCCAGGc                                                                                             >  1:288763/1‑111 (MQ=255)
                                                                                              ccttggccttggGCAACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGGGTGGTTTCGCCGCGCACTTCAATGATGAACTCCTGCTCCAGGc                                                                                             <  2:288763/111‑1 (MQ=255)
                                                                                                    ccgtggGCAACTCGCGGGCAGTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGGGTGGTTTCGCCGCGCACTTCAATGATGAACTCCTGCTCCAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTgg                                                 <  6:302222/146‑1 (MQ=255)
                                                                                                     cttggGCAACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGGGTGGTTTCGCCGCGCACTTCAATGATGAACTCCTGCTCCAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTggg                                                >  6:290546/1‑149 (MQ=255)
                                                                                                      ttggGCAACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGGGTGGTTTCGCCGCGCACTTCAATGATGAACTCCTGCTCCAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTggg                                                >  3:94495/1‑148 (MQ=255)
                                                                                                      ttggGCAACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGGGTGGTTTCGCCGCGCACTTCAATGATGAACTCCTGCTCCAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTgg                                                 <  5:244702/147‑1 (MQ=255)
                                                                                                      ttggGCAACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGGGTGGTTTCGCCGCGCACTTCAATGATGAACTCCTGCTCCAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTgg                                                 >  6:244702/1‑147 (MQ=255)
                                                                                                            aaCTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGGGTGGTTTCGCCGCGCACTTCAATGATGAACTCCTGCTCCAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTaaa                                            >  8:304339/1‑146 (MQ=255)
                                                                                                            aaCTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGGGTGGTTTCGCCGCGCACTTCAATGATGAACTCCTGCTCCAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTaaa                                            >  8:299555/1‑146 (MQ=255)
                                                                                                            aaCTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGGGTGGTTTCGCCGCGCACTTCAATGATGAACTCCTGCTCCAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTaaa                                            <  7:299555/146‑1 (MQ=255)
                                                                                                            aaCTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGGGTGGTTTCGCCGCGCACTTCAATGATGAACTCCTGCTCCAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTaaa                                            <  7:304339/146‑1 (MQ=255)
                                                                                                               tCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGGGTGGTTTCGCCGCGCACTTCAATGATGAACTCCTGCTCCAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGc                                                          <  4:231038/129‑1 (MQ=255)
                                                                                                               tCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGGGTGGTTTCGCCGCGCACTTCAATGATGAACTCCTGCTCCAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGc                                                          >  3:231038/1‑129 (MQ=255)
                                                                                                               tCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGGGTGGTTTCGCCGCGCACTTCAATGATGAACTCCTGCTCCAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTg                                      >  1:130851/1‑149 (MQ=255)
                                                                                                                 gcgGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGGGTGGTTTCGCCGCGCACTTCAATGATGAACTCCTGCTCCAGGCGACGCAGGTCGGCCTCGATCTTGc                                                                     <  7:232763/116‑1 (MQ=255)
                                                                                                                 gcgGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGGGTGGTTTCGCCGCGCACTTCAATGATGAACTCCTGCTCCAGGCGACGCAGGTCGGCCTCGATCTTGc                                                                     >  8:232763/1‑116 (MQ=255)
                                                                                                                 gcgGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGGGTGGTTTCGCCGCGCACTTCAATGATGAACTCCTGCTCCAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTg                                      >  6:81095/1‑147 (MQ=255)
                                                                                                                 gcgGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGGGTGGTTTCGCCGCGCACTTCAATGATGAACTCCTGCTCCAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGc                                     <  6:239319/148‑1 (MQ=255)
                                                                                                                            gcATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGGGTGGTTTCGCCGCGCACTTCAATGATGAACTCCTGCTCCAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGgcg                         >  7:137115/1‑149 (MQ=255)
                                                                                                                              aTGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGGGATGGTTTCGCCGCGCACTTCAATGATGAACTCCTGCTCCAGGCGACGCAGGTCGGCCCCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGgcgc                        >  1:64001/1‑149 (MQ=255)
                                                                                                                              aTGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGGGTGGTTTCGCCGCGCACTTCAATGATGAACTCCTGCTCCAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGgcgc                        >  7:20815/1‑148 (MQ=255)
                                                                                                                                    ccACGCGCCAGGCGTTCCAGCGCCTGGGGGGTGGTTTCGCCGCGCACTTCAATGATGAACTCCTGCTCCAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTc                 <  2:130851/149‑1 (MQ=255)
                                                                                                                                         cgcCAGGCGTTCCAGCGCCTGGGGGGTGGTTTCGCCGCGCACTTCAATGATGAACTCCTGCTCCAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCg            <  8:137115/149‑1 (MQ=255)
                                                                                                                                         ccccAGGCGCTGGAACGCCTGGGGGGTGGTTTCGCCGCGCACTTCAATGATGAACTCCTGCTCCAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCg            >  4:282820/3‑149 (MQ=255)
                                                                                                                                          gcCAGGCGTTCCAGCGCCTGGGGGGTGGTTTCGCCGCGCACTTCAATGATGAACTCCTGCTCCAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCgg           >  6:266641/1‑149 (MQ=255)
                                                                                                                                                    ccAGCGCCTGGGGGGTGGTTTCGCCGCGCACTTCAATGATGAACTCCTGCTCCAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCa  <  5:87375/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGCTCCAGGCGCAGGTAGGCGCATAGCTCGGAAATCTGTCCTGGTTGTGGGTTTTTCAACACCATGCGCAGGTGTGTCTCGTTCTGCCAGCTGGCCTTGGCCTTGGGCAACTCGCGGTCATTGCGCATGGCCCCACGCGCCAGGCGTTCCAGCGCCTGGGGGGTGGTTTCGCCGCGCACTTCAATGATGAACTCCTGCTCCAGGCGACGCAGGTCGGCCTCGATCTTGCGGGATACCCGCCAGTCCTGGGTAAACACCTGCAGGCCGCTGGCGCCACGCTCCAGCGGTGCCACGCA  >  NC_002947/1366083‑1366378

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: