Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 274,154 A→T 5.9% L139Q (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,1540AT 5.9% 199.9 / 2.7 68L139Q (CTG→CAG) PP_0221methionine ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (30/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 = 3.36e-01

TCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAGGCACAAGGTTGAAGCCTGCAGCA  >  NC_002947/274007‑274300
                                                                                                                                                     |                                                                                                                                                       
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‑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                                                                                                                                    <  1:332645/125‑1 (MQ=255)
                                         ccgAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATc                                                                                                                                    >  2:332645/3‑127 (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)
                                                                                                                                                   gcagGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAGGCACAAGGTTGAAGCCTgcagca  <  1:25797/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
TCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGG‑TGATCAAC‑CCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTG‑TTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAGGCACAAGGTTGAAGCCTGCAGCA  >  NC_002947/274007‑274300

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