Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,366,560 C→A 7.0% R49L (CGC→CTC)  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,5600CA 7.0% 179.3 / 2.8 58R49L (CGC→CTC) PP_1191S4 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (24/29);  new base A (2/2);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

CTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCGGACATGGGG  >  NC_002947/1366423‑1366709
                                                                                                                                         |                                                                                                                                                     
cTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCCGGACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCaa                                                                                                                                           >  1:608436/1‑150 (MQ=255)
 ttCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAAc                                                                                                                                          <  2:939583/150‑1 (MQ=255)
 ttCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAAc                                                                                                                                          >  2:923723/1‑150 (MQ=255)
    actttGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCCGCAGGGTCACCGGCTCCAGGGTTTCTGCACGGTCGCCCGGCGTCAGGTCGACGCGCGGCGGCTCAGCCtt                                                                                                                                       >  1:1007320/5‑150 (MQ=255)
      ggTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGa                                                                                                                                     >  1:1018705/1‑150 (MQ=255)
       gTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGaa                                                                                                                                    <  1:780122/150‑1 (MQ=255)
          ggCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTg                                                                                                                                 <  1:1005041/150‑1 (MQ=255)
           gCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGc                                                                                                                                <  1:278490/150‑1 (MQ=255)
            ccTCGCTGAGGCTGCCCATGTTCATACTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCg                                                                                                                               >  1:30816/1‑150 (MQ=255)
             cTCGCTGAGGCTGCCCCTGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGTCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCgg                                                                                                                              <  2:566135/150‑1 (MQ=255)
             cTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCgg                                                                                                                              >  1:157625/1‑150 (MQ=255)
                           ccATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGAccacca                                                                                                                >  2:261782/1‑150 (MQ=255)
                                ttCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCa                                                                                                           >  2:1056620/1‑150 (MQ=255)
                                  cATGCTGGCTCGGGCATCTTCTCTGTCCATGCCGGTTGTCTGGTTTAGCAGCAGGGTCACCTGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCAGGCCATc                                                                                                         <  2:228676/150‑1 (MQ=255)
                                      cTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGAcc                                                                                                     <  2:1018705/150‑1 (MQ=255)
                                       tGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCGGTCGACCg                                                                                                    >  1:47320/1‑150 (MQ=255)
                                           gcgGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTAc                                                                                                <  1:615917/150‑1 (MQ=255)
                                               gCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTc                                                                                                                        <  2:903839/122‑1 (MQ=255)
                                               gCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTc                                                                                                                        >  1:903839/1‑122 (MQ=255)
                                               gCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAg                                                                                            >  1:530810/1‑150 (MQ=255)
                                                    ttCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGac                                                                                                                     >  2:751500/1‑120 (MQ=255)
                                                    ttCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGac                                                                                                                     <  1:751500/120‑1 (MQ=255)
                                                    ttCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGcctcc                                                                                       >  1:305571/1‑150 (MQ=255)
                                                    ttCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGcctcc                                                                                       >  1:839113/1‑150 (MQ=255)
                                                      cACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCtt                                                                                     <  2:305571/150‑1 (MQ=255)
                                                        cTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCGCGTCGAgctgctgcc                                                                                                                                              >  1:886792/1‑90 (MQ=39)
                                                        cTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCGCGTCGAgctgctgcc                                                                                                                                              <  2:886792/91‑2 (MQ=39)
                                                           tCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATg                                                                                <  2:30816/150‑1 (MQ=255)
                                                                    gctggctgGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTCCAATGTACTGCTCg                                                                       <  1:2075/150‑1 (MQ=255)
                                                                            gTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGAccacca                                                                                                                <  2:722933/101‑1 (MQ=255)
                                                                            gTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGAccacca                                                                                                                >  1:722933/1‑101 (MQ=255)
                                                                            gTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTcgcg                                                               <  1:1128399/150‑1 (MQ=255)
                                                                             ttCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCgcgg                                                              >  2:163368/1‑150 (MQ=255)
                                                                              tCAGCAGCAGGTTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCCGCTCGACGCGCTGCCGATTAAAAATGCAGTGCGTCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCgcggc                                                             <  2:587493/150‑1 (MQ=255)
                                                                              tCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCgcggc                                                             <  1:774927/150‑1 (MQ=255)
                                                                                  cagcagGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGa                                                         >  2:332338/1‑150 (MQ=255)
                                                                                       gggTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGagcagc                                                    >  2:119164/1‑150 (MQ=255)
                                                                                          tCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGa                                                 <  1:261782/150‑1 (MQ=255)
                                                                                           cACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAg                                                <  2:852949/150‑1 (MQ=255)
                                                                                            aCCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTTGCAGTGACGGGGTTCAACCTTGAACTGCGGCGGCGCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAgc                                               >  1:193109/1‑150 (MQ=255)
                                                                                                 cTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCATCCGAGATGCTc                                          <  2:1007320/150‑1 (MQ=255)
                                                                                                        gTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAgg                                   >  1:1106216/1‑150 (MQ=255)
                                                                                                        gTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAgg                                   <  1:923723/150‑1 (MQ=255)
                                                                                                         tttCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGc                                  <  1:133429/150‑1 (MQ=255)
                                                                                                              gCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTg                             <  1:332338/150‑1 (MQ=255)
                                                                                                                aCGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGCGCAGCCCATCTTCTCGATGAGGCGTTTGGa                           >  2:528323/1‑150 (MQ=255)
                                                                                                                   ggCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAg                        <  1:528323/150‑1 (MQ=255)
                                                                                                                        ccGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCgcggc                                                             >  1:775229/1‑108 (MQ=255)
                                                                                                                        ccGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCgcggc                                                             <  2:775229/108‑1 (MQ=255)
                                                                                                                        ccGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGAt                                       <  2:198616/130‑1 (MQ=255)
                                                                                                                        ccGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGAt                                       >  1:198616/1‑130 (MQ=255)
                                                                                                                           gcagcaGCGCGGCGAGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACgg                >  2:776870/1‑150 (MQ=255)
                                                                                                                             agcagcGCGTCGAGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGt              <  1:776870/150‑1 (MQ=255)
                                                                                                                               cagcTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTc            <  1:1056620/150‑1 (MQ=255)
                                                                                                                                   tCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGa                                                 >  2:130235/1‑109 (MQ=255)
                                                                                                                                   tCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGa                                                 <  1:130235/109‑1 (MQ=255)
                                                                                                                                     gACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCGGACAt      >  1:395447/1‑150 (MQ=255)
                                                                                                                                      aCGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCGGACATg     <  1:1114580/150‑1 (MQ=255)
                                                                                                                                        gcgcTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCGGACATggg   <  1:606298/150‑1 (MQ=255)
                                                                                                                                         cgcTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCGGACATgggg  <  1:768042/150‑1 (MQ=255)
                                                                                                                                         cgcTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCGGACATgggg  <  1:119164/150‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                     
CTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCGGACATGGGG  >  NC_002947/1366423‑1366709

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