Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,366,548 T→G 7.0% E51A (GAG→GCG)  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,5480TG 7.0% 222.2 / 7.1 72E51A (GAG→GCG) PP_1191S4 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (31/36);  new base G (3/2);  total (34/38)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.04e-01

AGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCG‑TTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCGG  >  NC_002947/1366404‑1366695
                                                                                                                                                |                                                                                                                                                    
aGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGAc                                                                                                                                                  <  8:47313/149‑1 (MQ=255)
 gggCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGAcg                                                                                                                                                 >  2:330491/1‑149 (MQ=255)
   gcgcgcACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGAcgcg                                                                                                                                               >  3:349728/1‑149 (MQ=255)
      cgcACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCGATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCgctg                                                                                                                                            <  6:380452/149‑1 (MQ=255)
        cACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGTGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCgctgct                                                                                                                                          <  6:377630/149‑1 (MQ=255)
           gCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCgctgctgct                                                                                                                                       <  4:349728/149‑1 (MQ=255)
           gCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCgctgctgct                                                                                                                                       <  4:16665/149‑1 (MQ=255)
             cTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTc                                                                                                                                      >  8:58492/1‑148 (MQ=255)
             cTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCGCGTCGAgctgctgccgg                                                                                                                                     <  7:36345/149‑4 (MQ=255)
                 gcgGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCt                                                                                                                                 <  6:193605/149‑1 (MQ=255)
                             gCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCgg                                                                                                                       <  7:58492/147‑1 (MQ=255)
                             gCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGgc                                                                                                                      <  4:186924/148‑1 (MQ=255)
                               tGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGgctgc                                                                                                                   >  3:163897/1‑149 (MQ=255)
                                 aGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCTCGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGgctgc                                                                                                                   <  7:74732/147‑1 (MQ=255)
                                   gCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAAcc                                                                                                                                  <  1:222608/130‑1 (MQ=255)
                                   gCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAAcc                                                                                                                                  >  2:222608/1‑130 (MQ=255)
                                    cTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGac                                                                                                              <  7:266753/149‑1 (MQ=255)
                                      gCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGacca                                                                                                            >  2:2808/1‑149 (MQ=255)
                                        ccATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGaccac                                                                                                           >  5:377127/1‑148 (MQ=255)
                                        ccATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGaccac                                                                                                           >  8:101612/1‑148 (MQ=255)
                                         cATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCCCGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCGTGAACTGCGGATGCTCGaccacaa                                                                                                         >  8:343111/1‑147 (MQ=255)
                                                    tGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGAcc                                                                                              >  1:74587/1‑149 (MQ=255)
                                                    tGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGAcc                                                                                              >  7:136345/1‑149 (MQ=255)
                                                      gcgcGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGAcc                                                                                              <  1:210732/147‑1 (MQ=255)
                                                      gcgcGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGAcc                                                                                              >  2:210732/1‑147 (MQ=255)
                                                         cgGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCG‑TTAc                                                                                        >  4:186665/1‑149 (MQ=255)
                                                                cTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCG‑TTACCCAGcct                                                                                 >  6:355841/1‑149 (MQ=255)
                                                                   cACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCCTCGCCCG‑TTACCCAGcctcct                                                                              >  5:356000/1‑149 (MQ=255)
                                                                     cTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCG‑TTACCCAGCCTCCtt                                                                             >  2:167397/1‑148 (MQ=255)
                                                                      tGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCG‑TTACCCAGCCTCCTTc                                                                            >  4:137529/1‑148 (MQ=255)
                                                                        tCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTTCTGCTCAACCTTGAACTGCGGCTGCTGGACCACCACGCCATCGACCG‑TTACCCAGCCTCCTTCAAt                                                                         <  6:190038/149‑1 (MQ=255)
                                                                        tCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCG‑TTAc                                                                                        <  8:114767/134‑1 (MQ=255)
                                                                        tCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCG‑TTAc                                                                                        >  7:114767/1‑134 (MQ=255)
                                                                        tCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCG‑TTACCCAGCCTCCTTCAAt                                                                         >  7:1795/1‑149 (MQ=255)
                                                                            tGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCG‑TTACCCAGCCTCCTTCAATGTAc                                                                     <  1:2808/149‑1 (MQ=255)
                                                                              ccGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCG‑TTACCCAGCCTCCTTCAATGTACTg                                                                   >  7:230925/1‑149 (MQ=255)
                                                                               cGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCG‑TTACCCAGCCTCCTTCAATGTACTGc                                                                  <  4:224326/149‑1 (MQ=255)
                                                                                ggctggctGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCG‑TTACCCAGCCTCCTTCAATGTACTGCt                                                                 >  4:154048/1‑149 (MQ=255)
                                                                                 gctggctgGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCG‑TTACCCAGCCTCCTTCAATGTACTGCTc                                                                <  6:377127/149‑1 (MQ=255)
                                                                                  ctggctggTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCG‑TTACCCAGCCTCCTTCAATGTACTGCTCg                                                               <  3:23944/149‑1 (MQ=255)
                                                                                    ggctggTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCG‑TTACCCAGCCTCCTTCAATGTACTGCTCGGc                                                             <  8:136345/149‑1 (MQ=255)
                                                                                     gctggTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCG‑TTACCCAGCCTCCTTCAATGTACTGCTCGGc                                                             <  1:330491/148‑1 (MQ=255)
                                                                                       tggTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCG‑TTACCCAGCCTCCTTCAATGTACTGCTCGGCCt                                                           >  3:89264/1‑148 (MQ=255)
                                                                                       tggTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCG‑TTACCCAGCCTCCTTCAATGTACTGCTCGGCCTc                                                          >  4:87124/1‑149 (MQ=255)
                                                                                          ttCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCG‑TTACCCAGCCTCCTTCAATGTACTGCTCGGCCTcgcg                                                       <  2:288294/149‑1 (MQ=255)
                                                                                              gcagcagGGTTACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCG‑TTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGcggcgg                                                   <  1:319216/149‑1 (MQ=255)
                                                                                              gcagcagGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCG‑TTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGcggcgg                                                   <  1:167397/149‑1 (MQ=255)
                                                                                              gcagcagGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCG‑TTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGcggcgg                                                   <  2:146483/149‑1 (MQ=255)
                                                                                                  cagGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCG‑TTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGcggcgg                                                   >  3:252954/1‑145 (MQ=255)
                                                                                                  cagGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCG‑TTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGcggcgg                                                   <  4:252954/145‑1 (MQ=255)
                                                                                                    gggTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCG‑TTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGagcag                                             <  4:163897/149‑1 (MQ=255)
                                                                                                         aCCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCG‑TTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAg                                        >  5:111685/1‑149 (MQ=255)
                                                                                                               tCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCG‑TTACCGAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTc                                  >  5:119236/1‑149 (MQ=255)
                                                                                                                   gggTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGACTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGAt                               <  1:297636/149‑1 (MQ=255)
                                                                                                                   gggTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCG‑TTACCAAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGAt                               >  8:316066/1‑148 (MQ=255)
                                                                                                                    ggTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCG‑TTACCTAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGa                             <  6:119236/149‑1 (MQ=255)
                                                                                                                    ggTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCG‑TTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGa                             <  1:377903/149‑1 (MQ=255)
                                                                                                                    ggTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCG‑TTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGa                             <  8:1795/149‑1 (MQ=255)
                                                                                                                     gTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCG‑TTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAg                            <  3:154048/149‑1 (MQ=255)
                                                                                                                           gCAAGGGCGCCCGGCAGCAGCTCGATGCGCTGCAGCTCAACCTTGAACTGCGGCTGCTCGACCACCCCGCCATGGACCG‑TTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTAGATTAGGCGgtt                      >  4:50025/1‑149 (MQ=255)
                                                                                                                            cACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCG‑TTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTg                     >  2:235208/1‑149 (MQ=255)
                                                                                                                               gggCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCG‑TTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGAc                  <  6:111685/149‑1 (MQ=255)
                                                                                                                                ggCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCG‑TTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACa                 <  5:355841/149‑1 (MQ=255)
                                                                                                                                 gcgcCCGGCAGCAGCTCGACGCGCTGCTGCTCAACATTGAACTGCGGCTGCTCGCCCACCACGCCATCGACCG‑TTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGAAAg                >  1:130459/1‑149 (MQ=255)
                                                                                                                                   gcCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCG‑TTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGc              <  7:101612/149‑1 (MQ=255)
                                                                                                                                       ggCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCG‑TTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGAc          >  6:178584/1‑149 (MQ=255)
                                                                                                                                        gcagcaGCGCGTCGAGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCG‑TTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACg         >  2:39850/1‑149 (MQ=255)
                                                                                                                                        gcagcaGCGCGTCGAGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCG‑TTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACg         >  3:266602/1‑149 (MQ=255)
                                                                                                                                         cagcagCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCG‑TTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACgg        <  3:186665/149‑1 (MQ=255)
                                                                                                                                          agcagcTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCG‑TTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACggg       <  2:96061/149‑1 (MQ=255)
                                                                                                                                           gcagcGCGTCGAGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCG‑TTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGcc                                           <  8:207329/112‑1 (MQ=255)
                                                                                                                                           gcagcGCGTCGAGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCG‑TTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGcc                                           >  7:207329/1‑112 (MQ=255)
                                                                                                                                                tCGACGCGCTGCTGCTAAACCTTGAACTGCGGCTGCTCGACCAACACGCCATCGACCG‑TTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCATCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCgt  >  8:186275/1‑147 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
AGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCG‑TTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCGG  >  NC_002947/1366404‑1366695

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