Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,939,498 G→A 6.6% A89T (GCA→ACA)  PP_1739 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,939,4980GA 6.6% 174.5 / 2.6 60A89T (GCA→ACA) PP_1739membrane protein
Reads supporting (aligned to +/- strand):  ref base G (28/28);  new base A (2/2);  total (30/30)
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.10e-01

GGGCTGTGCGCAGTGCCGACCGCCTATTGGCAGTGCGAGGT‑GCGTTCGCTGCTGGGGTTGGGCATTCACCATCAAGTCATGGCCTGGGGCGGATTGGGCCTTGCGCTACTGGCCCAGTTCGTTTCCGCGCGTGCCGGCTGGTGGCTGGCAGTGTTGGCACTGGTGTTCG‑CGGTGCC‑GGCGCTGGTGCTGTACACCGCCAGCATTGCGGTGTTTGCGCTGGTGCTGGCGGGGTTGCGCCTGGTCAGGCAACCCCCATTGTAAGCAGCTGCGGGATCGCGCCAGCTTTCAGAAATTGA  >  NC_002947/1939351‑1939646
                                                                                                                                                    |                                                                                                                                                      
gggCTGTGCGCAGTGCCGACCGCCTATTGGCAGTGCGAGGT‑GCGTTCGCTGCTGGGGTTGGGCATTCACCATCAAGTCATGGCCTGGGGCGGATTGGGCCTTGCGCTACTGGCCCAGTTCGTTTCCGCGCGTGCCGGCTGGtggctgg                                                                                                                                                        <  5:136275/148‑1 (MQ=255)
      tgCGCAGTGCCGACCGCCTATTGGCAGTGCGAGGT‑GCGTTCGCTGCTGGGGTTGGGCATTCACCATCAAGTCATGGCCTGGGGCGGATTGGGCCTTGCGCTACTGGCCCAGTTCGTTTCCGCGCGTGCCGGCTGGTGGCTGGCAGTGtt                                                                                                                                                 <  8:230612/149‑1 (MQ=255)
      tgCGCAGTGCCGACCCCCTATTGGCAGTGCGAGGT‑GCGTTCGCTGCTGGGGTTGGGCATTCACCATCAAGTCATGGCCTGGGCCGGATTGGGCCTTGCGCTACTGCCCCAGGACGTTTCCGCGCGTGCCGGCTCGTGGCTGGCAGTGtt                                                                                                                                                 >  4:101500/1‑149 (MQ=255)
         gcAGTGCCGACCGCCTATTGGCAGTGCGAGGT‑GCGTTCGCTGCTGGGGTTGGGCATTCACCATCAAGTCATGGCCTGGGGCGGATTGGGCCTTGCGCTACTGGCCCAGTTCGTTTCCGCGCGTGCCGGCTGGTGGCTGGCAGTGTTgg                                                                                                                                               >  1:188127/1‑148 (MQ=255)
             tGCCGACCGCCTATTGGCAGTGCGAGGTGGCGTTCGCTGCTGGGGTTGGGCATTCACCATCAAGTCATGGCCTGGGGCGGATTGGGCCTTGCGCTACTGGCCCAGTTCGTTTCCGCGCGTGCCGGCTGGTGGCTGGCa                                                                                                                                                      >  1:221455/1‑138 (MQ=255)
             tGCCGACCGCCTATTGGCAGTGCGAGGT‑GCGTTCGCTGCTGGGGTTGGGCATTCACCATCAAGTCATGGCCTGGGGCGGATTGGGCCTTGCGCTACTGGCCCAGTTCGTTTCCGCGCGTGCCGGCTGGTGGCTGGCa                                                                                                                                                      <  2:221455/137‑1 (MQ=255)
             tGCCGACCGCCTATTGGCAGTGCGAGGT‑GCGTTCGCTGCTGGGGTTGGGCATTCACCATCAAGTCATGGCCTGGGGCGGATTGGGCCTTGCGCTACTGGCCCAGTTCGTTTCCGCGCGTGCCGGCTGGTGGCTGGCAGTGTTGGCACTg                                                                                                                                          <  7:17374/149‑1 (MQ=255)
               ccGACCGCCTATTGGCAGTGCGAGGT‑GCGTTCGCTGCTGGGGTTGGGCATTCACCATCAAGTCATGGCCTGGGGCGGATTGGGCCTTGCGCTACTGGCCCAGTTCGTTTCCGCGCGTGCCGGCTGGTGGCTGGCAGTGTTGGCACTg                                                                                                                                          >  6:168112/1‑147 (MQ=255)
               ccGACCGCCTATTGGCAGTGCGAGGT‑GCGTTCGCTGCTGGGGTTGGGCATTCACCATCAAGTCATGGCCTGGGGCGGATTGGGCCTTGCGCTACTGGCCCAGTTCGTTTCCGCGCGTGCCGGCTGGTGGCTGGCAGTGTTGGCACTg                                                                                                                                          <  5:168112/147‑1 (MQ=255)
                    cGCCTATTGGCAGTGCGAGGT‑GCGTTCGCTGCTGGGGTTGGGCATTCACCATCAAGTCATGGCCTGGGGCGGATTGGGCCTTGCGCTACTGGCCCAGTTCGTTTCCGCGCGTGCCGGCTGGTGGCTGGCAGTGTTGGCACTGGTGTTcg                                                                                                                                   >  3:341222/1‑149 (MQ=255)
                      ccTATTGGCAGTGCGAGGT‑GCGTTCGCTGCTGGGGTTGGGCATTCACCATCAAGTCATGGCCTGGGGCGGATTGGGCCTTGCGCTACTGGCCCAGTTCGTTTCCGCGCGTGCCGGCTGGTGGCTGGCAGTGTTGGCACTGGTGTTcg‑cg                                                                                                                                >  3:205872/1‑149 (MQ=255)
                      ccTATTGGCAGTGCGAGGT‑GCGTTCGCTGCTGGGGTTGGGCATTCACCATCAAGTCATGGCCTGGGGCGGATTGGGCCTTGCGCTACTGGCCCAGTTCGTTTCCGCGCGTGCCGGCTGGTGGCTGACAGTGTTGGCACTGGTGTTcg‑cg                                                                                                                                >  5:20705/1‑149 (MQ=255)
                       cTATTGGCAGTGCGAGGT‑GCGTTCGCTGCTGGGGTTGGGCATTCACCATCAAGTCATGGCCTGGGGCGGATTGGGCCTTGCGCTACTGGCCCAGTTCGTTTCCGCGCGTGCCGGCTGGTGGCTGGCAGTGTTGGCACTGGTGTTCG‑Cgg                                                                                                                               <  2:20874/149‑1 (MQ=255)
                          ttGGCAGTGCGAGGT‑GCGTTCGCTGCTGGGGTTGGGCATTCACCATCAAGTCATGGCCTGGGGCGGATTGGGCCTTGCGCTACTGGCCCAGTTCGTTTCCGCGCGTGCCGGCTGGTGGCTGGCAGTGTTGGCACTGGTGTTCG‑CGGt                                                                                                                              <  8:179452/147‑1 (MQ=255)
                          ttGGCAGTGCGAGGT‑GCGTTCGCTGCTGGGGTTGGGCATTCACCATCAAGTCATGGCCTGGGGCGGATTGGGCCTTGCGCTACTGGCCCAGTTCGTTTCCGCGCGTGCCGGCTGGTGGCTGGCAGTGTTGGCACTGGTGTTCG‑CGGTg                                                                                                                             <  4:341222/148‑1 (MQ=255)
                          ttGGCAGTGCGAGGT‑GCGTTCGCTGCTGGGGTTGGGCATTCACCATCAAGTCATGGCCTGGGGCGGATTGGGCCTTGCGCTACTGGCCCAGTTCGTTTCCGCGCGTGCCGGCTGGTGGCTGGCAGTGTTGGCACTGGTGTTCG‑CGGTg                                                                                                                             <  7:41682/148‑1 (MQ=255)
                          ttGGCAGTGCGAGGT‑GCGTTCGCTGCTGGGGTTGGGCATTCACCATCAAGTCATGGCCTGGGGCGGATTGGGCCTTGCGCTACTGGCCCAGTTCGTTTCCGCGCGTGCCGGCTGGTGGCTGGCAGTGTTGGCACTGGTGTTCG‑CGGTGc                                                                                                                            <  5:251917/149‑1 (MQ=255)
                            ggCAGTGGGAGGT‑GCGTTCGCTGCTGGGGTTGGGCATTCAGCATCAAGTCATGGGCTGGGGCGGATTGGGCCTTGCGCTACTGGCCCTGTTCGTTTCCGCGCTTGCCGGCTGGTGGCTGGCAGTGtt                                                                                                                                                 <  2:186306/127‑1 (MQ=255)
                            ggCAGTGCGAGGT‑GCGTTCGCTGCTGGGGATGGGCATTCACCATCAAGTCATGGCCTGGGGCGGATTGGGCCTTGCGCTACTGGCCCAGTTCGTTTCCGCGCGTGCCGGCTGGTGGCTGGCAGTGtt                                                                                                                                                 >  1:186306/1‑127 (MQ=255)
                             gCAGTGCGAGGT‑GCGTTCGCTGCTGGGGTTGGGCATTCACCATCAAGTCATGGCCTGGGGCGGATTGGGCCTTGCGCTACTGGCCCAGTTCGTTTCCGCGCGTGCCGGCTGGTGGCTGGCAGTGTTg                                                                                                                                                >  6:136408/1‑127 (MQ=255)
                             gCAGTGCGAGGT‑GCGTTCGCTGCTGGGGTTGGGCATTCACCATCAAGTCATGGCCTGGGGCGGATTGGGCCTTGCGCTACTGGCCCAGTTCGTTTCCGCGCGTGCCGGCTGGTGGCTGGCAGTGTTg                                                                                                                                                <  5:136408/127‑1 (MQ=255)
                                 tGCGAGGT‑GCGTTCGCTGCTGGGGTTGGGCATTCACCATCAAGTCATGGCCTGGGGCGGATTGGGCCTTGCGCTACTGGCCCAGTTCGTTTCCGCGCGTGCCGGCTGGTGGCTGACa                                                                                                                                                      <  5:163189/117‑1 (MQ=255)
                                 tGCGAGGT‑GCGTTCGCTGCTGGGGTTGGGCATTCACCATCAAGTCATGGCCTGGGGCGGATTGGGCCTTGCGCTACTGGCCCAGTTCGTTTCCGCGCGTGCCGGCTGGTGGCTGACa                                                                                                                                                      >  6:163189/1‑117 (MQ=255)
                                     aGGT‑GCGTTCGCTGCTGGGGTTGGGCATTCACCATCAAGTCATGGCCTGGGGCGGATTGGGCCTTGCGCTACTGGCCCAGTTCGTTTCCGCGCGTGCCGGCTGGTGGCTGGCAGTGTTGGCACTGGTGTTCG‑CGGTGCC‑GGCGCTGGTg                                                                                                                <  5:177222/149‑1 (MQ=255)
                                        t‑GCGTTCGCTGCTGGGGTTGGGCATTCACCATCAAGTCATGGCCTGGGGCGGATTGGGCCTTGCGCTACTGGCCCAGTTCGTTTCCGCGCGTGCCGGCTGGTGGGTGGGAGTGTTGGCACTGGTGTTCG‑CGGTGCC‑GGCGATGGTGCTg                                                                                                             >  7:163500/1‑149 (MQ=255)
                                          gCGTTCGCTGCTGGGGTTGGGCATTCACCATCAAGTCATGGCCTGGGGCGGATTGGGCCTTGCGCTACTGGCCCAGTTCGTTTCCGCGCGTGCCGGCTGGTGGCTGGCAGTGTTGGCACTGGTGTTCG‑CGGTGCC‑GGCGCTGGTGCTGt                                                                                                            <  3:266973/149‑1 (MQ=255)
                                             ttCGCTGCTGGGGTTGGGCATTCACCATCAAGTCATGGCCTGGGGCGGATTGGGCCTTGCGCTACTGGCCCAGTTCGTTTCCGCGCGTGCCGGCTGGTGGCTGGCAGTGTTGGCACTGGTGTTCG‑CGGTGCC‑GGCGCTGGTGCTGTaca                                                                                                         >  2:260506/1‑149 (MQ=255)
                                                  tgctgGGGTTGGGCATTCACCATCAAGTCATGGCCTGGGGCGGATTGGGCCTTGCGCTACTGGCCCAGTTCGTTTCCGCGCGTGCCGGCTGGTGGCTGGCAGTGTTGGCACTGGTGTTCG‑CGGTGc                                                                                                                            <  4:200913/126‑1 (MQ=255)
                                                  tgctgGGGTTGGGCATTCACCATCAAGTCATGGCCTGGGGCGGATTGGGCCTTGCGCTACTGGCCCAGTTCGTTTCCGCGCGTGCCGGCTGGTGGCTGGCAGTGTTGGCACTGGTGTTCG‑CGGTGc                                                                                                                            >  3:200913/1‑126 (MQ=255)
                                                         gTTGGGCATTCACCATCAAGTCATGGCCGGGGGCGGATTGGGCCTTGCGCTACTGGCCCAGTTCGTTTCCGCGCGTGCCGGCTGGTGGCTGGCAGTGTTGGCACTGGTGTTCG‑CGGTGCC‑GGCGCTGGTGCTGTACACCGc                                                                                                     <  7:263029/141‑1 (MQ=255)
                                                         gTTGGGCATTCACCATCAAGTCATGGCCGGGGGCGGATTGGGCCTTGCGCTACTGGCCCAGTTCGTTTCCGCGCGTGCCGGCTGGTGGCTGGCAGTGTTGGCACTGGTGTTCG‑CGGTGCC‑GGCGCTGGTGCTGTACACCGc                                                                                                     >  8:263029/1‑141 (MQ=255)
                                                             ggCATTCACCATCAAGTCATGGCCTGGGGCGGATTGGGCCTTGCGCTACTGGCCCAGTTCGTTTCCGCGCGTGCCGGCTGGTGGCTGGCAGTGTTGGCACTGGTGTTCG‑CGGTGCC‑GGCGCTGGTGCTGTACACCGCCa                                                                                                   <  8:326815/139‑1 (MQ=255)
                                                             ggCATTCACCATCAAGTCATGGCCTGGGGCGGATTGGGCCTTGCGCTACTGGCCCAGTTCGTTTCCGCGCGTGCCGGCTGGTGGCTGACAGTGTTGGCACTGGTGTTCG‑CGGTGCC‑GGCGCTGGTGCTGTACACCGCCAGCATTGCGgt                                                                                         <  6:20705/149‑1 (MQ=255)
                                                             ggCATTCACCATCAAGTCATGGCCTGGGGCGGATTGGGCCTTGCGCTACTGGCCCAGTTCGTTTCCGCGCGTGCCGGCTGGTGGATGGCAGTGTTGGCACTGGTGTTCG‑CGGTGCC‑GGCGCTGGTGCTGTACACCGCCa                                                                                                   >  7:326815/1‑139 (MQ=255)
                                                                aTTCACCATCAAGTCATGGCCTGGGGCGGATTGGGCCTTGCGCTACTGGCCCAGTTCGTTTCCGCGCGTGCCGGCTGGTGGCTGGCAGTGTTGGCACTGGTGTTCG‑CGGTGCC‑GGCGCTGGTGCTGTACACCGCCAGCATTGCGGTGtt                                                                                      <  6:65166/149‑1 (MQ=255)
                                                                aTTAACCATCAAGTCATGGCCTGGGGCGGATTGGGCCTTGCGCTACTGGCCCAGTTCGTTTCCGCGCGTGCCGGCTGGTGGCTGGCAGTGTTGGCACTGGTGTTCG‑CGGTGCC‑GGCGCTGGTGCTGTACACCGCCAGCATTGCGGTGtt                                                                                      >  3:150524/1‑149 (MQ=255)
                                                                    aCCATCAAGTCATGGCCTGGGGCGGATTGGGCCTTGCGCTACTGGCCCAGTTCGTTTCCGCGCGTGCCGGCTGGTGGCTGGCAGTGTTGGCACTGGTGTTCG‑CGGTGCC‑GGCGCTGGTGCTGTACACCGACAGCATTGCGGTGTTTgcg                                                                                  >  3:81239/1‑149 (MQ=255)
                                                                     ccATCAAGTCATGGCCTGGGGCGGATTGGGCCTTGCGCTACTGGCCCAGTTCGTTTCCGCGCGTGCCGGCTGGTGGCTGGCAGTGTTGGCACTGGTGTTCG‑CGGTGCC‑GGCGCTGGTGCTGTACACCGCCAGCATTGCGGTGTTTgcgc                                                                                 <  2:111540/149‑1 (MQ=255)
                                                                     ccATCAAGTCATGGCCTGGGGCGGATTGGGCCTTGCGCTACTGGCCCAGTTCGTTTCCGCGCGTGCCGGCTGGTGGCTGGCAGTGTTGGCACTGGTGTTCG‑CGGTGCC‑GGCGCTGGTGCTGTACACCGCCAGCATTGCGGTGTTTgcgc                                                                                 <  2:188127/149‑1 (MQ=255)
                                                                      cATCAAGTCATGGCCTGGGGCGGATTGGGCCTTGCGCTACTGGCCCAGTTCGTTTCCGCGCGTGCCGGCTGGTGGCTGGCAGTGTTg                                                                                                                                                <  1:34373/87‑1 (MQ=255)
                                                                      cATCAAGTCATGGCCTGGGGCGGATTGGGCCTTGCGCTACTGGCCCAGTTCGTTTCCGCGCGTGCCGGCTGGTGGCTGGCAGTGTTg                                                                                                                                                >  2:34373/1‑87 (MQ=255)
                                                                                    cTGGGGCGGATTGGGCCTTGCGCTACTGGCCCAGTTCGTTTCCGCGCGTGCCGGCTGGTGGCTGGCAGTGTTGGCACTGGTGTTCG‑CGGTGCC‑GGCGCTGGTGCTGTACACCGCCAGCATTGCGGTGTTTGCGCTGGTGCTGGCGGGGt                                                                  >  4:252123/1‑149 (MQ=255)
                                                                                         gCGGATTGGGCCTTGCGCTACTGGCCCAGTTCGTTTCCGCGCGTGCCGGCTGGTGGCTGGCAGTGTTGGCACTGGTGTTCG‑CGGTGCC‑GGCGCTGGTGCTGTACACCGCCAGCATTGCGGTGTTTGCGCTGGTGCTGGCGGGGTTgcgc                                                             >  4:127875/1‑149 (MQ=255)
                                                                                                 ggCCTTGCGCTACTGGCCCAGTTCGTTTCCGCGCGTGCCGGCTGGTGGATGGCAGTGTTGGCACTGGTGTTCG‑CGGTGCC‑GGCGCTGGTGCTGTACACCGCAAGCATTGCGGTGTTTGCGCTGGTGCTGGCGGGGTTGCGCCTGGTCAg                                                     >  5:239931/1‑149 (MQ=255)
                                                                                                   ccTTGCGCTACTGGCCCAGTTCGTTTCCGCGCGTGCCGGCTGGTGGCTGGCAGTGTTGGCACTGGTGTTCG‑CGGTGCC‑GGCGCTGGTGCTGTACACCGCCAGCATTGCGGTGTTTGCGCTGGTGCTGGCGGGGTTGCGCCTGGTCAgg                                                    >  8:359321/1‑148 (MQ=255)
                                                                                                   ccTTGCGCTACTGGCCCAGTTCGTTTCCGCGCGTGCCGGCTGGTGGCTGGCAGTGTTGGCACTGGTGTTCG‑CGGTGCC‑GGCGCTGGTGCTGTACACCGCCAGCATTGCGGTGTTTGCGCTGGTGCTGGCGGGGTTGCGCCTGGTCAgg                                                    >  6:315827/1‑148 (MQ=255)
                                                                                                                  ccAGTTCGTTTCCGCGCGTGGCGGCTGGTGGCTGGCAGTGTTGGCACTGGTGTTCG‑CGGTGCC‑GGCGCTGGTGCTGTACACCGCCAGCATTGCGGTGTTTGCGCTGGTGCTGGCGGGGTTGCGCCTGGTCAGGCAACCCCCATTGTAag                                    <  2:122091/149‑1 (MQ=255)
                                                                                                                  ccAGTTCGTTTCCGCGCGTGCCGGCTGGTGGCTGGCAGTGTTGGCACTGGTGTTCG‑CGGTGCC‑GGCGCTGGTGCTGTACACCGCCAGCATTGCGGTGTTTGCGCTGGTGCTGGCGGGGTTGCGCCTGGTCAGGCAACCCCCATTGTAag                                    >  2:141743/1‑149 (MQ=255)
                                                                                                                  ccAGTTCGTTTCCGCGCGTGCCGGCTGGTGGCTGGCAGTGTTGGCACTGGTGTTCG‑CGGTGCC‑GGCGCTGGTGCTGTACACCGCCAGCATTGCGGTGTTTGCGCTGGTGCTGGCGGGGTTGCGCCTGGTCAGGCAACCCCCATTGTAag                                    <  3:252123/149‑1 (MQ=255)
                                                                                                                  ccAGTTCGTTTCCGCGCGTGCCGGCTGGTGGCTGGCAGTGTTGGCACTGGTGTTCG‑CGGTGCC‑GGCGCTGGTGCTGTACACCGCCAGCATTACGGTGTTTGCGCTGGTGCTGGCGGGGTTGCGCCTGGTCAGGCAACCCCCATTGTAag                                    >  1:323064/1‑149 (MQ=255)
                                                                                                                  ccAGTTCGTTTCCGCGCGTGCCGGCTGGTGGCTGGCAGTGTTGGCACTGGTGTTCG‑CGGTGCC‑GGCGCTGGTGCTGCACACCGCCAGCATTGCGGTGTTTGCGCTGGTGCTGGCGGGGTTGCGCCTGGTCAGGCAACCCCCATTGTAag                                    >  2:104624/1‑149 (MQ=255)
                                                                                                                     gTTCGTTTCCGCGCGTGCCGGCTGGTGGCTGGCAGTGTTGGCACTGGTGTTCG‑CGGTGCC‑GGCGCTGGTGCTGTACACCGCCAGCATTGCGGTGTTTGCGCTGGTGCTGGCGGGGTTGCGCCTGGt                                                        >  3:122271/1‑126 (MQ=255)
                                                                                                                     gTTCGTTTCCGCGCGTGCCGGCTGGTGGCTGGCAGTGTTGGCACTGGTGTTCG‑CGGTGCC‑GGCGCTGGTGCTGTACACCGCCAGCATTGCGGTGTTTGCGCTGGTGCTGGCGGGGTTGCGCCTGGt                                                        <  4:122271/126‑1 (MQ=255)
                                                                                                                         gTTTCCGCGCGTGCCGGCTGGTGGCTGGCAGTGTTGGCACTGGTGTTCGACGTTGCC‑GGCGATGTTGATGTACACCGCCAGCATTGCGGTGTTTGCGGTGGTCATGGCGGGGTTCCGCATGGTCAGGCAACCCCCATTTTAGGCAGATg                              >  1:88504/1‑149 (MQ=255)
                                                                                                                            tCCGCGCGTGCCGGCTGGTGGCTGGCAGTGTTGGCACTGGTGTTCG‑CGGTGCC‑GGCGCTGGTGCTGTACACCGCCAGCATTGCGGTGTTTGCGCTGGTGCTGGCGGGGTTGCGCCTGGTCAGGCAACCCCCATTGTAAGCAGCTGCgg                           <  1:260506/148‑1 (MQ=255)
                                                                                                                               gcgcgTGCCGGCTGGTGGCTGGCAGTGTTGGCACTGGTGTTCG‑CGGTGCC‑GGCGCTGGTGCTGTACACCGCCAGCATTGCGGTGTTTGCGCTGGTGCTGGCGGGGTTGCGCCTGGTCAGGCAACCCCCATTGTAAGCAGCTGCGGGATc                       >  3:297960/1‑149 (MQ=255)
                                                                                                                                cgcgTGCCGGCTGGTGGCTGGCAGTGTTGGCACTGGTGTTCG‑CGGTGCC‑GGCGCTGGTGCTGTACACCGCCAGCATTGCGGTGTTTGCGCTGGTGCTGGCGGGGTTGCGCCTGGTCAGGCAACCCCCATTGTAAGCAGCTGCGGGATcg                      >  6:378311/1‑149 (MQ=255)
                                                                                                                                 gcgTGCCGGCTGGTGGCTGGCAGTGTTGGCACTGGTGTTCG‑CGGTGCC‑GGCGCTGGTGCTGTACACCGCCAGCATTGCGGTGTTTGCGCTGGTGCTGGCGGGGTTGCGCCTGGTCAGGCAACCCCCATTGTAAGCAGCTGCGGGATcgc                     <  4:81239/149‑1 (MQ=255)
                                                                                                                                    tGCCGGCTGGTGGCTGGCAGTGTTGGCACTGGTGTTCG‑CGGTGCC‑GGCGCTGGTGCTGTACACCGCCAGCATTGCGGTGTTTGCGCTGGTGCTGGCGGGGTTGCGCCTGGTCAGGCAACCCCCATTGTAAGCAGCTGCGGGATCGCGcc                  <  7:359321/149‑1 (MQ=255)
                                                                                                                                      ccGGCTGGTGGCTGGCAGTGTTGGCACTGGTGTTCG‑CGGTGCC‑GGCGCTGGTGCTGTACACCGCCAGCATTGCGGTGTTTGCGCTGGTGCTGGCGGGGTTGCGCCTGGTCAGGCAACCCCCATTGTAAGCAGCTGCGGGATCGCGCCAg                <  5:378311/149‑1 (MQ=255)
                                                                                                                                          cTGGTGGCTGGCAGTGTTGGCACTGGTGTTCG‑CGGTGCC‑GGCGCTGGTGCTGTACACCGCCAGCATTGCGGTGTTTGCGCTGGTGCTGGCGGGGTTGCGCCTGGTCAGGCAACCCCCATTGTAAGCAGCTGCGGGATCGCGCCAGCttt            <  5:315827/149‑1 (MQ=255)
                                                                                                                                                gctggcAGTGTTGGCACTGGTGTTCG‑CGGTGCC‑GGCGCTGGTGCTGTACACCGCCAGCATTGCGGTGTTTGCGCTGGTGCTGGCGGGGTTGCGCCTGGTCAGGCAACCCCCATTGTAAGCAGCTGCGGGATCGCGCCAGCTTTCAGaaa      >  6:370955/1‑149 (MQ=255)
                                                                                                                                                  tggcAGTGTTGGCACTGGTGTTCG‑CGGTGCCTGGCGCTGGTGCTGTACACCGCCAGCATTGCGGTGTTTGCGCTGGTGCTGGCGGGGTTGCGCCTGGTCAGGCAACCCCCATTGTAAGCAGCTGCGGGATCGCGCCAGCTTTCAGaaa      >  1:40336/1‑148 (MQ=255)
                                                                                                                                                  tggcAGTGTTGGCACTGGTGTTCG‑CGGTGCC‑GGCGCTGGTGCTGTACACCGCCAGCATTGCGGTGTTTGCGCTGGTGCTGGCGGGGTTGCGCCTGGTCAGGCAACCCCCATTGTAAGCAGCTGCGGGATCGCGCCAGCTTTCAGaaa      <  2:40336/147‑1 (MQ=255)
                                                                                                                                                  tggcAGTGTTGGCACTGGTGTTCG‑CGGTGCC‑GGCGCTGGTGCTGTACACCGCCAGCATTGCGGTGTTTGCGCTGGTGCTGGCGGGGTTGCGCCTGGTCAGGCAACCCCCATTGTAAGCAGCTGCGGGATCGCGCCAGCTTTCAGAAAtt    >  2:171791/1‑149 (MQ=255)
                                                                                                                                                    gcAGTGTTGGCACTGGTGTTCG‑CGGTGCC‑GGCGCTGGTGCTGTACACCGCCAGCATTGCGGTGTTTGCGCTGGTGCTGGCGGGGTTGCGCCTGGTCAGGCAACCCCCATTGTAAGCAGCTGCGGGATCGCGCCAGCTTTCAGAAAtt    >  8:287872/1‑147 (MQ=255)
                                                                                                                                                    gcAGTGTTGGCACTGGTGTTCG‑CGGTGCC‑GGCGCTGGTGCTGTACACCGCCAGCATTGCGGTGTTTGCGCTGGTGCTGGCGGGGTTGCGCCTGGTCAGGCAACCCCCATTGTAAGCAGCTGCGGGATCGCGCCAGCTTTCAGAAAtt    >  5:250140/1‑147 (MQ=255)
                                                                                                                                                    gcAGTGTTGGCACTGGTGTTCG‑CGGTGCC‑GGCGCTGGTGCTGTACACCGCCAGCATTGCGGTGTTTGCGCTGGTGCTGGCGGGGTTGCGCCTGGTCAGGCAACCCCCATTGTAAGCAGCTGCGGGATCGCGCCAGCTTTCAGAAATTGa  <  6:239931/149‑1 (MQ=255)
                                                                                                                                                    gcAGTGTTGGCACTGGTGTTCG‑CGGTGCC‑GGCGCTGGTGCTGTACACCGCCAGCATTGCGGTGTTTGCGCTGGTGCTGGCGGGGTTGCGCCTGGTCAGGCAACCCCCATTGTAAGCAGCTGCGGGATCGCGCCAGCTTTCAGAAATTGa  <  8:91029/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
GGGCTGTGCGCAGTGCCGACCGCCTATTGGCAGTGCGAGGT‑GCGTTCGCTGCTGGGGTTGGGCATTCACCATCAAGTCATGGCCTGGGGCGGATTGGGCCTTGCGCTACTGGCCCAGTTCGTTTCCGCGCGTGCCGGCTGGTGGCTGGCAGTGTTGGCACTGGTGTTCG‑CGGTGCC‑GGCGCTGGTGCTGTACACCGCCAGCATTGCGGTGTTTGCGCTGGTGCTGGCGGGGTTGCGCCTGGTCAGGCAACCCCCATTGTAAGCAGCTGCGGGATCGCGCCAGCTTTCAGAAATTGA  >  NC_002947/1939351‑1939646

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