Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 274,151 A→T 5.9% L140Q (CTG→CAG)  PP_0221 ← methionine ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947274,1510AT 5.9% 206.2 / 2.7 68L140Q (CTG→CAG) PP_0221methionine ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (31/33);  new base T (2/2);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.79e-01

GGTCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAGGCACAAGGTTGAAGCCTGC  >  NC_002947/274005‑274296
                                                                                                                                                    |                                                                                                                                                      
ggTCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCTCCTCGACCTGAAGGATCTTAAGCTTTTTAGGGTTGGCGAGTATGTCAGCCTAGGTGTACTTGTAGCACACGCCGGGATTGATGG‑TGATCAAC‑CCGGCCTTGGCcagcagcag                                                                                                                                                       >  8:203036/1‑148 (MQ=255)
  tCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCcagcagcagcag                                                                                                                                                    <  3:68953/149‑1 (MQ=255)
    ggCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGc                                                                                                                                                  >  8:329053/1‑149 (MQ=255)
       gtcgtcgtAGGCGCGTACCAGTTGCACCGCCTCGACCTGATGTATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCaa                                                                                                                                               <  2:314053/149‑2 (MQ=255)
          gtcgtAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCcagcagca                                                                                                                                                        <  3:219749/137‑1 (MQ=255)
          gtcgtAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCcagcagca                                                                                                                                                        >  4:219749/1‑137 (MQ=255)
          gtcgtAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCACGGc                                                                                                                                            >  4:413188/1‑149 (MQ=255)
                ggCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAACAGCGGCCGTGGCCTGCTGCTGCTGGCCAAGGCCggg                                                                                                                                        <  6:182733/148‑3 (MQ=37)
                ggCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑t                                                                                                                                     >  2:270691/1‑149 (MQ=255)
                 gcgcGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑tt                                                                                                                                    >  3:267145/1‑149 (MQ=255)
                      tACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAACAGCGGCCGTGGCCTGCTGCTGCTGGCCAAGGCCG‑GG‑TTGATca                                                                                                                               <  8:242476/149‑2 (MQ=25)
                        ccAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCggg                                                                                                                             >  2:270644/1‑149 (MQ=255)
                               cACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTgg                                                                                                                      <  3:77184/149‑1 (MQ=255)
                               cACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCCGGTCGTTgg                                                                                                                      >  7:228909/1‑149 (MQ=255)
                                   gCCTCGACCTGAAGGATCATCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGAt                                                                                                                  <  5:12758/149‑1 (MQ=255)
                                      tCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGc                                                                                                               >  5:241528/1‑149 (MQ=255)
                                       cGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGc                                                                                                               >  2:25797/1‑148 (MQ=255)
                                        gACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGc                                                                                                               <  1:270691/147‑1 (MQ=255)
                                           ccgAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATc                                                                                                                                >  2:332645/3‑127 (MQ=255)
                                           ccgAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATc                                                                                                                                <  1:332645/125‑1 (MQ=255)
                                              aaGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTgg                                                                                                       >  7:391253/1‑149 (MQ=255)
                                              aaGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTgg                                                                                                       >  1:48334/1‑149 (MQ=255)
                                               aGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGgc                                                                                                      <  2:135035/149‑1 (MQ=255)
                                                   tCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCt                                                                                                  >  6:139130/1‑149 (MQ=255)
                                                     ttCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCAGGGCTTGAGGG‑TGATCAAC‑ACGGCCTTGGCCAGCATCAGCAGGCCACGGCAGATG‑TTGATCTGGTCGTTGGCGATGGCCACCCTGGCGCCtt                                                                                                 >  5:338347/1‑148 (MQ=255)
                                                      tCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCtt                                                                                                 >  4:240069/1‑147 (MQ=255)
                                                       cAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGTGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGTCGCCtt                                                                                                 <  7:203036/146‑1 (MQ=255)
                                                       cAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGATTAAAAAA‑CCCGCCTTTGGCAAGAGGAGCCCGCCACGGGCGGGGTTTGAACGGTTCTTTGGCGATGGCCACGATgccccat                                                                                                  >  1:131063/1‑141 (MQ=255)
                                                       cAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGcc                                                                                                   <  4:267145/144‑1 (MQ=255)
                                                       cAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGcc                                                                                                   <  6:199104/144‑1 (MQ=255)
                                                       cAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCtt                                                                                                 <  4:205913/146‑1 (MQ=255)
                                                       cAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCt                                                                                                  <  7:163102/145‑1 (MQ=255)
                                                       cAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTc                                                                                                <  3:426178/147‑1 (MQ=255)
                                                       cAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCa                                                                                               <  3:240069/148‑1 (MQ=255)
                                                        aGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAAGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAg                                                                                              <  2:423271/148‑1 (MQ=255)
                                                         gCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCATTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAAA‑CCGGCCTTGGCCAGCAGAAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAgg                                                                                             >  7:293551/1‑148 (MQ=255)
                                                                 gggTTGGCGAGGATGTCAGCCTCGGTGGGCTTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTCGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCgg                                                                                 <  7:80512/148‑1 (MQ=255)
                                                                  ggTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGc                                                                                >  2:348727/1‑148 (MQ=255)
                                                                    ttGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCa                                                                               >  4:62431/1‑147 (MQ=255)
                                                                          aGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCtttt                                                                       <  5:127518/149‑1 (MQ=255)
                                                                           ggATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTg                                                                      <  3:399424/149‑1 (MQ=255)
                                                                              tggcAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTAc                                                                   <  3:246977/146‑1 (MQ=255)
                                                                                  aGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACtttt                                                               <  7:329053/149‑1 (MQ=255)
                                                                                    ccTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTc                                                             >  7:290557/1‑149 (MQ=255)
                                                                                         gTGGCCTTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTa                                                        >  6:249400/1‑149 (MQ=255)
                                                                                          tGGCCTTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTa                                                        >  7:212331/1‑148 (MQ=255)
                                                                                          tGGCCTTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTa                                                        >  6:151034/1‑148 (MQ=255)
                                                                                           ggCCTTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTa                                                        >  5:192648/1‑147 (MQ=255)
                                                                                            gCCTTGTAGCCGACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAg                                                     <  5:160233/149‑1 (MQ=255)
                                                                                            gCCTTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCCGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAg                                                     >  6:266241/1‑149 (MQ=255)
                                                                                             ccTTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAgg                                                    >  5:178949/1‑149 (MQ=255)
                                                                                              cTTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGc                                                   <  3:413188/149‑1 (MQ=255)
                                                                                                 gTAGCCCACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGcc                                                  <  5:247121/147‑1 (MQ=255)
                                                                                                      ccACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACgttg                                           >  7:29753/1‑149 (MQ=255)
                                                                                                       gaCGCCTGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCAGCAGCATCAGGCCACGGACGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACgttgt                                          <  8:293551/148‑1 (MQ=255)
                                                                                                           ccGGGCTTGAGGG‑GGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTtgat                                      <  6:241528/149‑1 (MQ=255)
                                                                                                            cGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTtgatg                                     <  4:9791/149‑1 (MQ=255)
                                                                                                               gCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGatgatg                                  >  2:43160/1‑149 (MQ=255)
                                                                                                                     ggg‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCggg                            <  6:192648/149‑1 (MQ=255)
                                                                                                                       g‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCTCA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCggg                            <  1:343246/149‑1 (MQ=255)
                                                                                                                           aTCAAC‑CCGGGCTTTGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTTATGATGCCCGGGGCGt                        <  6:338347/148‑1 (MQ=255)
                                                                                                                            tCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGt                        >  8:228068/1‑147 (MQ=255)
                                                                                                                              aaC‑CCGGCCTTGGCCAGCCGCAGCAGGCCACGGCCGCTT‑GTGATCGGGTCGTTGGCGTTGGCCACGCTGGCGCCTTCAGG‑‑CACTGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTaa                      <  1:348727/149‑1 (MQ=255)
                                                                                                                                c‑ccGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAGGca                  >  3:37774/1‑149 (MQ=255)
                                                                                                                                  gcGGCCGTGGCCTGCTGCTGCTGGCCAAGGCCG‑GG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAGGcac                 >  8:343721/2‑148 (MQ=38)
                                                                                                                                   cGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAGGcaca                <  6:178949/149‑1 (MQ=255)
                                                                                                                                   cGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAGGcac                 >  1:212053/1‑148 (MQ=255)
                                                                                                                                       cgtGGCCTGCTGCTGCTGGCCAAGGCCG‑GG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAGGCACAAGGtt           >  7:288013/3‑149 (MQ=38)
                                                                                                                                          ggCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAGGCACAAGGTTg          <  7:228068/148‑1 (MQ=255)
                                                                                                                                           gCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAGGCACAAGGTTGaa        <  8:29753/149‑1 (MQ=255)
                                                                                                                                            ccAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAgg                    >  1:128553/1‑136 (MQ=255)
                                                                                                                                            ccAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAgg                    <  2:128553/136‑1 (MQ=255)
                                                                                                                                             cagcagcagcagGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAGGCACAAGGTTGAAGc      <  8:212331/149‑1 (MQ=255)
                                                                                                                                             cagcagcagcagGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAGGCACAAGGTTGAAGc      >  8:53940/1‑149 (MQ=255)
                                                                                                                                              agcagcagcagGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAGGCACAAGGTTGAAGcc     <  2:131063/149‑1 (MQ=255)
                                                                                                                                                cagcagcagGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTATTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAGGCACAAGGTTGAAGCCTg   >  1:101046/1‑149 (MQ=255)
                                                                                                                                                 agcagcagGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAGGCACAAGGTTGAAGCCTgc  <  3:30364/149‑1 (MQ=255)
                                                                                                                                                   cagcagGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACTTTGTTGATGATGCCCGGGACGTAAGGCACAAGGTTGAAGCCTg   >  3:80342/1‑146 (MQ=255)
                                                                                                                                                    agcagGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCggg                            >  4:315385/1‑120 (MQ=255)
                                                                                                                                                    agcagGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCggg                            <  3:315385/120‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
GGTCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAGGCACAAGGTTGAAGCCTGC  >  NC_002947/274005‑274296

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