Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,864,704 T→C 5.9% I227V (ATC→GTC)  gllP ← porin‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,864,7040TC 5.9% 197.9 / 2.3 68I227V (ATC→GTC) gllPporin‑like protein
Reads supporting (aligned to +/- strand):  ref base T (33/31);  new base C (2/2);  total (35/33)
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.01e-01

GAGAACAGGCCCGAAGCGGTGCGGTTGTTCAGTTCGCCTGCCTTGGCGCCACCGTCGTCCTTGCCCCAGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCC‑AGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAA‑GTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTCCTGCATGCTGGCGTCGTTGCGCGGGCTG  >  NC_002947/2864556‑2864852
                                                                                                                                                     |                                                                                                                                                     
gagaACAGGCCCGAAGCGGTGCGGTTGTTCAGTTCGCCTGCCTTGGCGCCACCGTCGTCCTTGCCCCAGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCC‑AGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGAt                                                                                                                                                       <  4:130546/149‑1 (MQ=255)
     cAGGCCCGAAGCGGTGCGGTTGTTCAGTTCGCCTGCCTTGGCGCCACCGTCGTCCTTGCCCCAGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCC‑AGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCtt                                                                                                                                                  <  8:12226/149‑1 (MQ=255)
      aGGCCCGAAGCGGTGCGGTTGTTCAGTTCGCCTGCCTTGGCGCCACCGTCGTCCTTGCCCCAGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCC‑AGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTc                                                                                                                                                 >  4:162223/1‑149 (MQ=255)
             aaGCGGTGCGGTTGTTCAGTTCGCCTGCCTTGGCGCCACCGTCGTCCTTGCCCCAGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCC‑AGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTgg                                                                                                                                          <  5:367822/149‑1 (MQ=255)
              aGCGGTGCGGTTGTTCAGTTCGCCTGCCTTGGCGCCACCGTCGTCCTTGCCCCAGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGACC‑AGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTggg                                                                                                                                         <  8:55720/149‑1 (MQ=255)
              aGCGGTGCGGTTGTTAAGTTCGCCTGCCTTGGGGGCACCGTCGTCCTTGCCCCACTAGTATCCCAGGTTTGCGCCCAGTACCCAGTTGCCCTTGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGAAGACATCTACCAGcagc                                                                                                                                          <  1:106744/149‑4 (MQ=255)
                   tgcggtTGTTCAGTTCGCCTGCCTTGGCGCCACCGTCGTCCTTGCCCCAGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCC‑AGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATAc                                                                                                                                    <  2:313360/149‑1 (MQ=255)
                       gttgttCAGTTCGCCTGCCTTGGCTCCACCGTCGTCCTTGCCCCAGAAGTAGCCCAGGTGGGCGCCCAGTACCCAGTTGCAC‑AGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATAc                                                                                                                                    <  7:181160/145‑1 (MQ=255)
                       gttgttCAGTTCGCCTGCCTTGGCGCCACAGTCGTCCTTGCCCCAGAAGTAGACCAGGTTGGCGCCCAGTACCCAGTTGCCC‑AGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCAGGTAGATGTATTCCAGCTGGGCATaa                                                                                                                                    >  8:181160/1‑144 (MQ=255)
                           ttCAGTTCGCCTGCCTTGGCGCCACCGTCGTCCTTGCCCCAGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCC‑AGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGcc                                                                                                                            >  6:319093/1‑149 (MQ=255)
                           ttCAGTTCGCCTGCCTTGGCGCCACCGTCGTCCTTGCCCCAGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCC‑AGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGGAGATGTCTTCCAGCTGGGCATACCAGGCGcc                                                                                                                            >  6:313831/1‑149 (MQ=255)
                            tCAGTTCGCCTGCCTTGGCGCCACCGTCGTCCTTGCCCCAGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCC‑AGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGAAGACATCTACCAGcagcgct                                                                                                                                       >  3:351781/1‑131 (MQ=39)
                            tCAGTTCGCCTGCCTTGGCGCCACCGTCGTCCTTGCCCCAGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCC‑AGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGAAGACATCTACCAGcagcgct                                                                                                                                       <  4:351781/137‑7 (MQ=39)
                              aGTTCGCCTGCCTTGGCGCCACCGTCGTCCTTGCCCCAGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCC‑AGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGAc                                                                                                                         >  7:35678/1‑149 (MQ=255)
                               gTTCGCCTGCCTTGGCGCCACCGTCGTCCTTGCCCCAGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTTCCC‑AGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGAcc                                                                                                                        <  4:209803/149‑1 (MQ=255)
                                ttCGCCTGCCTTGGCGCCACCGTCGTCCTTGCCCCAGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCC‑AGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGAt                                                                                                                                                       >  5:33604/1‑117 (MQ=255)
                                ttCGCCTGCCTTGGCGCCACCGTCGTCCTTGCCCCAGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCC‑AGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGAt                                                                                                                                                       <  6:33604/117‑1 (MQ=255)
                                  cGCCTGCCTTGGCGCCACCGTCGTCCTTGCCCCAGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCC‑AGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCt                                                                                                                                                   >  4:30506/1‑119 (MQ=255)
                                  cGCCTGCCTTGGCGCCACCGTCGTCCTTGCCCCAGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCC‑AGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCt                                                                                                                                                   <  3:30506/119‑1 (MQ=255)
                                  cGCCTGCCTTGGCGCCACCGTCGTCCTTGCCCCAGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCC‑AGCGGCTGTTGATCCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCCGGCGCCGACCAgg                                                                                                                     >  6:384982/1‑149 (MQ=255)
                                   gcctgcctTGGCGCCACCGTCGTCCTTGCCCCAGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCC‑AGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAgg                                                                                                                     <  5:282816/148‑1 (MQ=255)
                                     ctgcctTGGCGCCACCGTCGTCCTTGCCCCAGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCC‑AGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGAt                                                                                                                                                       >  5:309388/1‑112 (MQ=255)
                                     ctgcctTGGCGCCACCGTCGTCCTTGCCCCAGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCC‑AGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGAt                                                                                                                                                       <  6:309388/112‑1 (MQ=255)
                                           tGGCGCCACCGTCGTCCTTGCCCCAGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCC‑AGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTcc                                                                                                            >  4:360964/1‑149 (MQ=255)
                                                 cACCGTCGTCCTTGCCCCAGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCC‑AGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGtt                                                                                                       >  8:14722/1‑148 (MQ=255)
                                                  aCCGTCGTCCTTGCCCCAGAAGTAGCCCATGTTGGCGCCCAGTACCCAGTTTCCC‑AGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGa                                                                                                     <  4:37995/149‑1 (MQ=255)
                                                        gtcCTTGCCCCAGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCC‑AGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCCGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCgg                                                                                                <  6:56710/148‑1 (MQ=255)
                                                                cccAGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCC‑AGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTcgccg                                                                                       >  5:44767/1‑149 (MQ=255)
                                                                 ccAGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCC‑AGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTcgccgc                                                                                      >  7:272076/1‑149 (MQ=255)
                                                                 ccAGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCC‑AGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTcgccgc                                                                                      >  7:362635/1‑149 (MQ=255)
                                                                 ccAGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCC‑AGCGGCTGTTGATGCTGCACATGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTcgccgc                                                                                      >  5:200490/1‑149 (MQ=255)
                                                                   aGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCC‑AGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGccgccg                                                                                    <  1:108046/149‑1 (MQ=255)
                                                                      aGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCC‑AGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGc                                                                                  >  1:204698/1‑148 (MQ=255)
                                                                      aGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCC‑AGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCg                                                                                 <  8:35678/149‑1 (MQ=255)
                                                                         aGCCCAGGTTGGCGCCCAGTACCCAGTTGCCC‑AGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAA‑g                                                                             >  4:301212/1‑149 (MQ=255)
                                                                               ggTTGGCGCCCAGTACCCAGTTGCCC‑AGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAA‑GTTGaaa                                                                       <  6:200490/149‑1 (MQ=255)
                                                                                    gcgcCCAGTACCCAGTTTCCC‑AGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAA‑GTTGAAACGTTc                                                                  <  6:217565/149‑1 (MQ=255)
                                                                                     cgcCCAGTACCCAGTTGCCC‑AGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAA‑GTTGAAACGGTcc                                                                 >  8:82594/1‑149 (MQ=255)
                                                                                         cAGTACCCAGTTGCCC‑AGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAA‑GTTGAAACGGTCCGAgg                                                             >  2:300672/1‑149 (MQ=255)
                                                                                           gTACCCAGTTGCCC‑AGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAA‑GTTGAAACGGTCCGAGGTg                                                           <  8:32968/149‑1 (MQ=255)
                                                                                               ccAGTTGCCC‑AGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAA‑GTTGAAACGGTCCGAGGTGAAgg                                                       >  8:173372/1‑149 (MQ=255)
                                                                                                   ttGCCC‑AGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAA‑GTTGAAACGGTCCGAGGTGAAGGCtgt                                                   <  8:272076/149‑1 (MQ=255)
                                                                                                     gCCC‑AGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAA‑GTTGAAACGCTCCGAGGTGAAGGCTGTgc                                                 >  8:271650/1‑149 (MQ=255)
                                                                                                      ccc‑AGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAA‑GTTGAAACGGTCCGAg                                                              >  2:69511/1‑135 (MQ=255)
                                                                                                      ccc‑AGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAA‑GTTGAAACGGTCCGAg                                                              <  1:69511/135‑1 (MQ=255)
                                                                                                       cc‑AGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAA‑GTTGAAACGGTCCGAGGTGAAGGCTGTgcgc                                               <  1:236968/149‑1 (MQ=255)
                                                                                                        c‑AGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAA‑GTTGAAACGGTCCGAGGTGAAGGCTGTGcgcg                                              <  7:173372/149‑1 (MQ=255)
                                                                                                          aGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAA‑GTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGt                                             >  2:211884/1‑149 (MQ=255)
                                                                                                                    gATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAA‑GTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCt                                    >  4:276738/1‑148 (MQ=255)
                                                                                                                      tgctgcACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAA‑GTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCAt                                 >  3:82567/1‑149 (MQ=255)
                                                                                                                       gctgcACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAA‑GTTGAAACGGTcc                                                                 <  4:111059/116‑1 (MQ=255)
                                                                                                                       gctgcACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAA‑GTTGAAACGGTcc                                                                 >  3:111059/1‑116 (MQ=255)
                                                                                                                       gctgcACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAA‑GTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATg                                >  2:204346/1‑149 (MQ=255)
                                                                                                                         tgcACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAA‑GTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTc                              >  4:336996/1‑149 (MQ=255)
                                                                                                                          gcACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGTTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAA‑GTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTcc                             >  1:265645/1‑149 (MQ=255)
                                                                                                                          gcACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAA‑GTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTcc                             >  6:172493/1‑149 (MQ=255)
                                                                                                                           cACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAA‑GTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTCCt                            <  2:204698/149‑1 (MQ=255)
                                                                                                                            aCCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAA‑GTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTCCTg                           <  2:222992/149‑1 (MQ=255)
                                                                                                                              cTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAA‑GTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTCCTGCa                         >  3:9118/1‑149 (MQ=255)
                                                                                                                               tGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAA‑GTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTCCTGCAt                        >  5:277180/1‑149 (MQ=255)
                                                                                                                                ggAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAA‑GTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTCCTGCATg                       >  3:6399/1‑149 (MQ=255)
                                                                                                                                  aaGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAA‑GTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTCCTGCATGCt                     <  3:276738/149‑1 (MQ=255)
                                                                                                                                    gTAGCGCGGCTGGTAGTTGTCTTCCTGGTGGGCATACCAGGGGCCGACCAGGGTGGGCTCCTGGTTGAAGGGGTATTCGCCGCCGGCGAATCTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTCCTGCATGCTg                    <  1:300672/149‑1 (MQ=255)
                                                                                                                                     tAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAA‑GTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTCCTGCATGCTGGc                  <  7:14722/149‑1 (MQ=255)
                                                                                                                                     tAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAA‑GTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTCCTGCATGCTGGc                  <  4:6399/149‑1 (MQ=255)
                                                                                                                                        cgcTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAA‑GTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTCCTGCATGCTGGcgtc               <  6:44767/149‑1 (MQ=255)
                                                                                                                                         gctgctGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAA‑GTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTCCTGCATGCTGGcgtcg              >  6:87729/1‑149 (MQ=255)
                                                                                                                                          ctgctgGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAA‑GTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTCCTGCATGCTGGcgtcgt             <  3:360964/149‑1 (MQ=255)
                                                                                                                                          cagctgGAAGACATCTACCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAA‑GTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTCCTGCATGCTGGcgtcgt             >  3:208543/3‑149 (MQ=255)
                                                                                                                                               ggTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAA‑GTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTCCTGCATGCTGGCGTCGTTgcgc        >  3:314846/1‑149 (MQ=255)
                                                                                                                                                gTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAA‑GTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTCCTGCATGCTGGCGTCGTTGcgcg       <  2:265645/149‑1 (MQ=255)
                                                                                                                                                   gATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGATGCGGTATTCGCCGCCGGCGAA‑GTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTCCTGCATGCTGGCGTCGTTGCGCGGGc    <  5:384982/149‑1 (MQ=255)
                                                                                                                                                   gATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAA‑GTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTCCTGCATGCTGGCGTCGTTGCGCGGGc    <  1:211884/149‑1 (MQ=255)
                                                                                                                                                     tGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAA‑GTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTCCTGCATGCTGGCGTCGTTGCGCGGGCTg  <  6:277180/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAGAACAGGCCCGAAGCGGTGCGGTTGTTCAGTTCGCCTGCCTTGGCGCCACCGTCGTCCTTGCCCCAGAAGTAGCCCAGGTTGGCGCCCAGTACCCAGTTGCCC‑AGCGGCTGTTGATGCTGCACCTGGAAGTAGCGCTGCTGGTAGATGTCTTCCAGCTGGGCATACCAGGCGCCGACCAGGGTGCGCTCCTGGTTGAAGCGGTATTCGCCGCCGGCGAA‑GTTGAAACGGTCCGAGGTGAAGGCTGTGCGCGTGAACAGGCTCATGTCCTGCATGCTGGCGTCGTTGCGCGGGCTG  >  NC_002947/2864556‑2864852

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