Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,031,919 T→A 6.2% L5M (TTG→ATG) ‡ PP_0889 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,031,9190TA 6.2% 187.6 / 2.7 65L5M (TTG→ATG) ‡PP_0889hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (27/34);  new base A (2/2);  total (29/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

CAATATCTCGGCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACGACCCTGCCACCCTCGACTGG  >  NC_002947/1031771‑1032065
                                                                                                                                                    |                                                                                                                                                   
cAATATCTCGGCCACCGCACGGGCGCTGGGCATGTACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCt                                                                                                                                                     >  8:61893/1‑149 (MQ=255)
 aaTATCTCGGCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCtt                                                                                                                                                    <  8:256405/149‑1 (MQ=255)
 aaTATCTCGGCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCtt                                                                                                                                                    >  8:136115/1‑149 (MQ=255)
  atatCTCGGCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCATCTGTGAGCGCGGGCTTg                                                                                                                                                   >  3:310505/1‑149 (MQ=255)
     tctcGGCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCAAGccc                                                                                                                                                >  4:308916/1‑149 (MQ=255)
       tcGGCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCcg‑c                                                                                                                                             >  6:185528/1‑149 (MQ=255)
         ggCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCAAGCCcg‑cgc                                                                                                                                           <  3:308916/149‑2 (MQ=255)
           ccACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCcg‑c                                                                                                                                             <  1:312932/145‑1 (MQ=255)
           ccACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCcg‑c                                                                                                                                             >  2:312932/1‑145 (MQ=255)
           ccACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAg                                                                                                                                         <  5:185528/149‑1 (MQ=255)
              ccGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGc                                                                                                                                      >  5:250448/1‑149 (MQ=255)
                 cACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTg                                                                                                                                    >  7:331706/1‑148 (MQ=255)
                  aCGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCa                                                                                                                                  <  8:302267/149‑1 (MQ=255)
                                  cACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCaa                                                                                                                  <  7:136115/149‑1 (MQ=255)
                                   aCCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAg                                                                                                                 >  6:126663/1‑149 (MQ=255)
                                    ccGACGGACGTTGCAGCGGAAGCTGCAGAAGCGGCCGGTTCGGCGCTGTTTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGGGCGCGCGGGCTTGGCCGACGAAGGGCTGCAAAGCAGTCCCAATCAAg                                                                                                                 <  1:26453/149‑1 (MQ=255)
                                     cGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCAAGCCcg‑cgc                                                                                                                                           <  6:22524/121‑2 (MQ=255)
                                     cGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCAAGCCcg‑cgc                                                                                                                                           >  5:22524/1‑120 (MQ=255)
                                         ggacgTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGc                                                                                                           <  7:61893/149‑1 (MQ=255)
                                             gTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGa                                                                                                       <  8:331706/149‑1 (MQ=255)
                                                gCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATg                                                                                                     <  6:250448/148‑1 (MQ=255)
                                                  aGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCATCTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAg                                                                                                  <  4:310505/149‑1 (MQ=255)
                                                           tGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTCTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTg                                                                                         >  1:110550/1‑149 (MQ=255)
                                                             cAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTg                                                                                         <  7:168815/147‑1 (MQ=255)
                                                             cAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTg                                                                                         >  8:168815/1‑147 (MQ=255)
                                                                    gtccGGTTCGGCGCTTATTCCTTATCGTTGATGTTGGCCTGTTCGCTGGGGTATCCGCTCCTGCAACTGTTAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTc                                                                                                         <  3:134296/122‑1 (MQ=38)
                                                                     gCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAg                                                                                                                               <  2:36890/101‑1 (MQ=255)
                                                                     gCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAg                                                                                                                               >  1:36890/1‑101 (MQ=255)
                                                                     gCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGa                                                                                                             <  5:57627/119‑1 (MQ=255)
                                                                     gCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGa                                                                                                             >  6:57627/1‑119 (MQ=255)
                                                                     gCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGa                                                                                                    <  1:347953/128‑1 (MQ=255)
                                                                     gCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGa                                                                                                    >  2:347953/1‑128 (MQ=255)
                                                                      ccGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTgc                                                                              <  5:126663/149‑1 (MQ=255)
                                                                             ggCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGtt                                                                       >  6:345355/1‑149 (MQ=255)
                                                                               cgcTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGc                                                                     <  6:155632/149‑1 (MQ=255)
                                                                                 cTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAg                                                                   <  7:320334/149‑1 (MQ=255)
                                                                                 cTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAg                                                                   <  6:171185/149‑1 (MQ=255)
                                                                                 cTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAg                                                                   <  7:324980/149‑1 (MQ=255)
                                                                                         ttATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATc                                                                                                                    >  5:215048/1‑92 (MQ=255)
                                                                                         ttATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATc                                                                                                                    <  6:215048/92‑1 (MQ=255)
                                                                                         ttATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTg                                                                                         <  3:146655/119‑1 (MQ=255)
                                                                                         ttATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTg                                                                                         >  4:146655/1‑119 (MQ=255)
                                                                                              gatgatGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTg                                                                                         <  2:98463/114‑1 (MQ=255)
                                                                                              gatgatGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTg                                                                                         >  1:98463/1‑114 (MQ=255)
                                                                                              gatgatGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGc                                                      >  7:295907/1‑149 (MQ=255)
                                                                                              gatgatGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGc                                                      >  8:165694/1‑149 (MQ=255)
                                                                                               atgatgTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTgc                                                                              <  5:48998/124‑1 (MQ=255)
                                                                                               atgatgTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTgc                                                                              >  6:48998/1‑124 (MQ=255)
                                                                                               atgatgTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCt                                                     <  8:295907/149‑1 (MQ=255)
                                                                                                     ttGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCACCAAAGCTGAAGAc                                               >  7:95822/1‑149 (MQ=255)
                                                                                                        gCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTg                                             <  2:110550/148‑1 (MQ=255)
                                                                                                         ccTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAgcgc                                                            >  3:103885/1‑132 (MQ=255)
                                                                                                         ccTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAgcgc                                                            <  4:103885/132‑1 (MQ=255)
                                                                                                               cgcgGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGAc                                      <  5:10215/148‑1 (MQ=255)
                                                                                                                   ggggAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGaa                                                  >  4:129321/1‑132 (MQ=255)
                                                                                                                   ggggAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGaa                                                  <  3:129321/132‑1 (MQ=255)
                                                                                                                       aaTCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTgg                             <  7:165694/149‑1 (MQ=255)
                                                                                                                            gCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGa                                                >  2:120154/1‑125 (MQ=255)
                                                                                                                            gCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGa                                                <  1:120154/125‑1 (MQ=255)
                                                                                                                                cTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCaaa                                                        >  3:279314/1‑113 (MQ=255)
                                                                                                                                cTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCaaa                                                        <  4:279314/113‑1 (MQ=255)
                                                                                                                                cTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACGAc                    <  5:345355/149‑1 (MQ=255)
                                                                                                                                       tgtgAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGc                                           <  6:194288/119‑1 (MQ=255)
                                                                                                                                       tgtgAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGc                                           >  1:136663/1‑119 (MQ=255)
                                                                                                                                       tgtgAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGc                                           <  2:136663/119‑1 (MQ=255)
                                                                                                                                       tgtgAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGc                                           >  5:194288/1‑119 (MQ=255)
                                                                                                                                            gcgcGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGAAGTGGCAGAACGACCCTGCCACCCTc        >  6:363431/1‑149 (MQ=255)
                                                                                                                                                  gCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACGACCCTGCCACCCTCGACTgg  >  6:358241/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CAATATCTCGGCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCG‑CGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACGACCCTGCCACCCTCGACTGG  >  NC_002947/1031771‑1032065

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