Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,366,551 A→T 9.1% V50D (GTC→GAC)  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,5510AT 9.1% 171.1 / 5.3 63V50D (GTC→GAC) PP_1191S4 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (26/31);  new base T (3/3);  total (29/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 = 9.70e-01

AGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCG  >  NC_002947/1366404‑1366694
                                                                                                                                                    |                                                                                                                                               
aGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCAACGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑Ac                                                                                                                                                >  8:17221/1‑149 (MQ=255)
 gggCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTCGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑Acg                                                                                                                                               >  1:135850/1‑149 (MQ=255)
  ggCGCGCACGCCTTCGCGGTGGGCCTCTCTGAGGATGCCCATGTTCATGCTGGCGCGGGCATCTTCAATGTCAATGACGGCTGTCTGTTTCAGCAGCAGTGTCACCGGATACAGGGTTTCGGCACCGGCGCCCGGCAGCAGCTCG‑Acgc                                                                                                                                              >  2:257443/1‑149 (MQ=255)
  ggCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑Acg                                                                                                                                               <  4:46639/148‑1 (MQ=255)
  ggCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGCGTTTCGGCCCGGGCGCCCGGCAGCCGCTCG‑Acgc                                                                                                                                              >  6:245641/1‑149 (MQ=255)
      cgcACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCgctg                                                                                                                                          <  1:90919/149‑1 (MQ=255)
      cgcACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCgctg                                                                                                                                          <  5:181923/149‑1 (MQ=255)
                  cggcggGCCGCGCTGAGGCGGCCCATCTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCtt                                                                                                                              <  8:203246/149‑1 (MQ=255)
                  cgGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCtt                                                                                                                              <  8:194034/149‑1 (MQ=255)
                   ggTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTg                                                                                                                             <  6:108724/149‑1 (MQ=255)
                     tGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGaa                                                                                                                           <  7:26375/149‑1 (MQ=255)
                     tGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGaa                                                                                                                           <  2:164434/149‑1 (MQ=255)
                         ccTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGc                                                                                                                       >  4:264290/1‑149 (MQ=255)
                         ccTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGc                                                                                                                       >  6:329897/1‑149 (MQ=255)
                         ccTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGc                                                                                                                       >  8:235529/1‑149 (MQ=255)
                           tCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCgg                                                                                                                     >  3:228276/1‑149 (MQ=255)
                             gCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCgg                                                                                                                     <  4:8284/147‑1 (MQ=255)
                               tGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGgctgc                                                                                                                 >  2:208960/1‑149 (MQ=255)
                                  ggCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGAACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTc                                                                                                               <  1:208960/149‑1 (MQ=255)
                                     ttcccATTTTCATGCTGTCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGacc                                                                                                           <  7:198883/147‑1 (MQ=255)
                                                aTGCTGGCGCGGGCTTCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATc                                                                                                <  2:105194/149‑1 (MQ=255)
                                                      gcgcGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCATGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCg                                                                                           <  6:162591/148‑1 (MQ=255)
                                                      gcgcGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGt                                                                                          <  3:264290/149‑1 (MQ=255)
                                                        gcgGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTAG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTa                                                                                        >  4:81137/1‑149 (MQ=255)
                                                                   cACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGcctc                                                                               <  7:235529/147‑1 (MQ=255)
                                                                           aTGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTa                                                                     >  1:170897/1‑149 (MQ=255)
                                                                           aTGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTa                                                                     >  2:135052/1‑149 (MQ=255)
                                                                            tGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTAc                                                                    <  4:228276/149‑1 (MQ=255)
                                                                                 gctggctgGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTc                                                               >  8:292398/1‑149 (MQ=255)
                                                                                   tggctggTTCAGCAGAAGGGTCACCGGCTCCAGGGTTTCGGCACTGGCTCCAGGCAGCAGCTCT‑AAGCGCTGATGATCAACCTTGAACTGTGTCTGCTCGACCACAACGCCATCTACCGTTACCCAGCCTCCTTCAATGTACTGCTctg                                                             <  1:257443/149‑3 (MQ=255)
                                                                                    ggctggTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGCGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCCCCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGc                                                            >  8:6872/1‑149 (MQ=255)
                                                                                     gctggTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGcc                                                           >  6:273383/1‑149 (MQ=255)
                                                                                     gctggTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGc                                                            >  8:143035/1‑148 (MQ=255)
                                                                                         gTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTcgc                                                       <  2:106733/149‑1 (MQ=255)
                                                                                             agcagcagGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCgcggcg                                                   <  5:245641/149‑1 (MQ=255)
                                                                                              gcagcagGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGcggcgg                                                  <  2:170897/149‑1 (MQ=255)
                                                                                              gcagcagGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACACAGCCTCCTTCACTGTACTGCTCGGCCTCGcggcgg                                                  >  8:70475/1‑149 (MQ=255)
                                                                                              gcagcagGGTCACCGGCTCCAGGGTTTCGGCACGGACGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACACCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGcggcgg                                                  <  7:298235/149‑1 (MQ=255)
                                                                                                agcagGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGa                                                <  5:329897/149‑1 (MQ=255)
                                                                                                     ggTCACCGGCTCCAGGGTTTAGGCACGGGCCCACGGCAGCAGCTAG‑ACGCGCTGGTGCTCAACCATGAACTGCGGCTGCTCGACCCCCACGCCATCGACCGTTAACCAGCCTCCTTCAATTTACTGCTCGGACTCGCGGCGCGagcacc                                           >  2:194465/1‑147 (MQ=255)
                                                                                                         aCCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTCAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAg                                       <  7:6872/149‑1 (MQ=255)
                                                                                                           cGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAgct                                     >  1:94965/1‑149 (MQ=255)
                                                                                                                 cAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTAGCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTc                                 <  7:292398/147‑1 (MQ=255)
                                                                                                                 cAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGa                                                <  4:130867/132‑1 (MQ=255)
                                                                                                                 cAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGa                                                >  3:130867/1‑132 (MQ=255)
                                                                                                                 cAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCTCGGCTGCTcc                                <  8:20115/148‑2 (MQ=255)
                                                                                                                 cAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGa                               >  1:193632/1‑149 (MQ=255)
                                                                                                                  aGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGGCCTGCGCGAt                              >  3:1369/1‑149 (MQ=255)
                                                                                                                     gTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCTCGGCTGCTcccgccg                           <  3:81137/149‑7 (MQ=255)
                                                                                                                             aCGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTg                    >  8:92609/1‑148 (MQ=255)
                                                                                                                                   gcCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGc             <  5:273383/149‑1 (MQ=255)
                                                                                                                                    cccGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCg            <  8:273677/149‑1 (MQ=255)
                                                                                                                                      cGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGa          <  2:193632/149‑1 (MQ=255)
                                                                                                                                       ggCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGAc         >  5:123555/1‑149 (MQ=255)
                                                                                                                                        gcagcaGCGCG‑TCGAGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAAt                                                                        >  7:36448/1‑85 (MQ=39)
                                                                                                                                        gcagcaGCGCG‑TCGAGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAAt                                                                        <  8:36448/85‑1 (MQ=39)
                                                                                                                                        gcagcaGCGCG‑TCGAGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACg        >  2:287336/1‑149 (MQ=255)
                                                                                                                                          agcagcGCG‑TCGAGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACggg      >  5:150433/1‑149 (MQ=255)
                                                                                                                                           gcagcTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGt     >  1:233835/1‑149 (MQ=255)
                                                                                                                                           gcagcGCG‑TCGAGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGt     <  1:287336/149‑1 (MQ=255)
                                                                                                                                             agcGCG‑TCGAGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTc   <  6:150433/149‑1 (MQ=255)
                                                                                                                                              gcTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCg  >  5:233314/1‑149 (MQ=255)
                                                                                                                                              gcTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCg  <  1:135052/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
AGGGCGCGCACGCCTTCGCGGTGGGCCTCGCTGAGGCTGCCCATGTTCATGCTGGCGCGGGCATCTTCACTGTCCATGCCGGCTGGCTGGTTCAGCAGCAGGGTCACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCG‑ACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCG  >  NC_002947/1366404‑1366694

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