Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 657,266 G→A 5.5% S233N (AGC→AAC)  PP_0564 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947657,2660GA 5.5% 278.4 / 2.4 90S233N (AGC→AAC) PP_0564hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (43/42);  new base A (3/2);  total (46/44)
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.02e-01

ATCAAGTCGCTGACCCGCTACCAGACCCGTCTGGCTGCCGGCGACAAGGGCGCTATCTTCTATACCCGCGCCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATGGCTGCAGATCGACAACGTGTTCTATGAAGCCCGTGGCCAAGCCTGGG  >  NC_002947/657124‑657409
                                                                                                                                              |                                                                                                                                                  
aTCAAGTCGCTGACCCGCTACCAGACCCGTCTGGCTGCCGGCGACAAGGGCGCTATCTTCTATACCCGCGCCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTgg                                                                                                                                           >  3:300112/1‑149 (MQ=255)
    aGTCGCTGACCCGCTACCAGACCCGTCTGGCTGCCGGCGACAAGGGCGCTATCTTCTATACCCGCGCCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCgg                                                                                                                                       >  6:210239/1‑149 (MQ=255)
    aGTCGCTGACCCGCTACCAGACCCGTCTGGCTGCCGGCGACAAGGGCGCTATCTTCTATACCCGCGCCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCgg                                                                                                                                       <  6:226561/149‑1 (MQ=255)
    aGTCGCTGACCCGCTACCAGACCCGTCTGGCTGCCGGCGACAAGGGCGCTATCTTCTATACCCGCGCCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAACG‑‑‑CTGGCCga                                                                                                                                       >  5:327317/1‑148 (MQ=255)
     gTCGCTGACCCGCTACCAGACCCGTCTGGCTGCCGGCGACAAGGGCGCTATCTTCTATACCCGCGCCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCggg                                                                                                                                      <  5:354501/149‑1 (MQ=255)
     gTCGCTGACCCGCTACCAGACCCGTCTGGCTGCCGGCGACAAGGGCGCTATCTTCTATACCCGCGCCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAACG‑‑‑CTGGCCgac                                                                                                                                      <  6:327317/149‑3 (MQ=255)
      tCGCTGACCCGCTACCAGACCCGTCTGGCTGCCGGCGACAAGGGCGCTATCTTCTATACCCGCGCCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGt                                                                                                                                     >  7:369554/1‑149 (MQ=255)
         cTGACCCGCTACCAGACCCGTCTGGCTGCCGGCGACAAGGGCGCTATCTTCTATACCCGCGCCGACAACCTGAACAACTGGCGGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCaa                                                                                                                                  <  3:116307/149‑1 (MQ=255)
            aCCCGCTACCAGACCCGTCTGGCTGCCGGCGACAAGGGCGCTATCTTCTATACCCGCGCCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAg                                                                                                                                 <  1:39833/147‑1 (MQ=255)
            aCCCGCTACCAGACCCGTCTGGCTGCCGGCGACAAGGGCGCTATCTTCTATACCCGCGCCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGc                                                                                                                                >  4:254742/1‑148 (MQ=255)
                    ccAGACCCGTCTGGCTGCCGGCGACAAGGGCGCTATCTTCTATACCCGCGCCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAAcacc                                                                                                                        <  2:328007/148‑1 (MQ=255)
                    ccAGACCCGTCTGGCTGCCGGCGACAAGGGCGCTATCTTCTATACCCGCGCCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TCGGCCGGGTCAAGCTCAAcacca                                                                                                                       >  5:133902/1‑149 (MQ=255)
                    ccAGACCCGTCTGGCTGCCGGCGACAAGGGCGCTATCTTCTATACCCGCGCCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGACCAGCG‑‑‑TTGGCCGGGTCAAGCTCAAcacca                                                                                                                       >  1:382270/1‑149 (MQ=255)
                       gACCCGTCTGGCTGCCGGCGACAAGGGCGCTATCTTCTATACCCGCGCCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGc                                                                                                                    >  2:21332/1‑149 (MQ=255)
                          ccGTCTGGCTGCCGGCGACAAGGGCGCTATCTTCTATACCCGCGCCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTg                                                                                                                  >  8:312789/1‑148 (MQ=255)
                          ccGTCTGGCTGCCGGCGACAAGGGCGCTATCTTCTATACCCGCGCCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGACCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGa                                                                                                                 <  2:382270/149‑1 (MQ=255)
                            gTCTGGCTGCCGGCGACAAGGGCGCTATCTTCTATACCCGCGCCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTg                                                                                                                  >  8:73355/1‑146 (MQ=255)
                                 gCTGCCGGCGACAAGGGCGCTATCTTCTATACCCGCGCCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGa                                                                                                          <  8:224654/149‑1 (MQ=255)
                                   tGCCGGCGACAAGGGCGCTATCTTCTATACCCGCGCCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCg                                                                                                                                        >  5:172388/1‑117 (MQ=255)
                                   tGCCGGCGACAAGGGCGCTATCTTCTATACCCGCGCCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCg                                                                                                                                        <  6:172388/117‑1 (MQ=255)
                                   tGCCGGCGACAAGGGCGCTATCTTCTATACCCGCGCCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAgt                                                                                                        >  2:117936/1‑149 (MQ=255)
                                      cGGCGACAAGGGCGCTATCTTCTATACCCGCGCCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTg                                                                                                                  >  2:212731/1‑136 (MQ=255)
                                      cGGCGACAAGGGCGCTATCTTCTATACCCGCGCCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTg                                                                                                                  <  1:212731/136‑1 (MQ=255)
                                        gCGACCAGGGCGCTATCGTCTCTACCCGCGCCGACAACCGGAACAACTGGCTGGGTGACGGCGCCACCCGGCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAgtcggtc                                                                                                   <  6:220308/149‑1 (MQ=255)
                                                 gcgcTATCTTCTATACCCGCGCCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGc                                                                                          <  5:210239/149‑1 (MQ=255)
                                                 gcgcTAGCTTCTATACCCGCGCCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGc                                                                                          <  7:151150/149‑1 (MQ=255)
                                                   gcTATCTTCTATACCCGCGCCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGc                                                                                          >  5:242687/1‑147 (MQ=255)
                                                         ttctATACCCGCGCCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGcgcc                                                                                  >  6:353070/1‑149 (MQ=255)
                                                         ttctATACCCGCGCCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGcgcc                                                                                  >  6:7041/1‑149 (MQ=255)
                                                         ttctATACCCGCGCCGACAACCTGAACAACTGGCTAGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGcgcc                                                                                  >  1:258185/1‑149 (MQ=255)
                                                            tataCCCGCGCCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGcgccg                                                                                 >  5:323422/1‑147 (MQ=255)
                                                               aCCCGCGCCGACAACCTGAACACCTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGt                                                                            >  6:198194/1‑149 (MQ=255)
                                                               aCCCGCGCCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGt                                                                            >  6:307970/1‑149 (MQ=255)
                                                                 ccGCGCCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAgt                                                                                                        >  5:296634/1‑119 (MQ=255)
                                                                 ccGCGCCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCACAACACCACGCTGAAGACCGAgt                                                                                                        <  6:296634/119‑1 (MQ=255)
                                                                  cgcgCCGACAACCTTAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGa                                                                         <  2:187522/149‑1 (MQ=255)
                                                                  cgcgCCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCGCCCTTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGt                                                                            <  1:117936/149‑1 (MQ=255)
                                                                   gcgcCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGAt                                                                        <  1:391434/149‑1 (MQ=255)
                                                                        gACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGTCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATgag                                                                     >  6:100748/1‑147 (MQ=255)
                                                                        gACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGGCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATgagga                                                                   <  8:369554/149‑1 (MQ=255)
                                                                             ccTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCtgg                                                              >  8:264334/1‑149 (MQ=255)
                                                                               tGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCtggtg                                                            >  2:201058/1‑149 (MQ=255)
                                                                                    aacTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAAc                                                       <  1:21332/149‑1 (MQ=255)
                                                                                     acTGGCTGGGTGACGGCGCCACCCGTCTCGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAAcc                                                      <  2:167845/149‑1 (MQ=255)
                                                                                      ctggctggGGGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAAccc                                                     <  4:300112/149‑1 (MQ=255)
                                                                                       tggctggGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAAccc                                                     >  4:318893/1‑148 (MQ=255)
                                                                                         gctggGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCAt                                                  >  3:133914/1‑149 (MQ=255)
                                                                                         gctggGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCAt                                                  <  3:254742/149‑1 (MQ=255)
                                                                                           tggGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATgg                                                >  8:88596/1‑149 (MQ=255)
                                                                                            gggTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATGGc                                               >  2:83735/1‑149 (MQ=255)
                                                                                             ggTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATgg                                                <  6:70615/147‑1 (MQ=255)
                                                                                             ggTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATGGCt                                              >  3:244202/1‑149 (MQ=255)
                                                                                             ggTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATGGCt                                              >  8:202355/1‑149 (MQ=255)
                                                                                               tGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATGGCTGc                                            <  6:323422/149‑1 (MQ=255)
                                                                                                 aCGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATGGCTGCa                                           <  5:7041/148‑1 (MQ=255)
                                                                                                   gTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATGGCTg                                             >  8:80618/1‑144 (MQ=255)
                                                                                                   gTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATGGCTg                                             <  6:135423/144‑1 (MQ=255)
                                                                                                   gTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATGGCTg                                             >  5:135423/1‑144 (MQ=255)
                                                                                                   gTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATGGCTGCAGAt                                        <  1:201058/149‑1 (MQ=255)
                                                                                                      gCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAAccc                                                     <  6:139839/133‑1 (MQ=255)
                                                                                                      gCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAAccc                                                     >  5:139839/1‑133 (MQ=255)
                                                                                                         aCCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATGGCTGCAGATCGACaa                                  >  7:257198/1‑149 (MQ=255)
                                                                                                         aCCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATGGCTGCAGATCGACaa                                  <  5:100748/149‑1 (MQ=255)
                                                                                                          cccGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATGGCTGCAGATCGACAAc                                 <  7:312789/149‑1 (MQ=255)
                                                                                                          cccGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGGCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATGGCTGCAGATCGACAAc                                 <  7:88596/149‑1 (MQ=255)
                                                                                                           ccGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAgtcgg                                                                                                     <  4:305592/80‑1 (MQ=255)
                                                                                                           ccGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAgtcgg                                                                                                     >  3:305592/1‑80 (MQ=255)
                                                                                                             gTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATGGCTGCAGATCGACAACgtg                              >  7:231811/1‑149 (MQ=255)
                                                                                                              tCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCACCACCCCGCGGAAGACCGAGTCGGTCGTTGCCGGCCCGGCGCCGCAGGGGGATGAGGCGCTGGTGGCGACCCCATGGCTGCAGATCGACAACgtgt                             >  4:9570/1‑149 (MQ=255)
                                                                                                                   gCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATGGCTGCAGATCGACAACGTGTTCTAt                        <  7:80618/149‑1 (MQ=255)
                                                                                                                   gCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATGGCTGCAGATCGACAACGTGTTCTAt                        <  7:73355/149‑1 (MQ=255)
                                                                                                                    cTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATGGCTGCAGATCGACAACGTGTTCTATg                       <  1:114523/149‑1 (MQ=255)
                                                                                                                    cTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATGGCTGCAGATCGACAACGTGTTCTATg                       <  4:244202/149‑1 (MQ=255)
                                                                                                                       gCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATGGCTGCAGATCGACAACGTGTTCTATGaa                     <  2:258185/148‑1 (MQ=255)
                                                                                                                       gCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATGGCTGCAGATCGACAACGTGTTCTATGAAg                    >  8:48030/1‑149 (MQ=255)
                                                                                                                        cTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATGGCTGCAGATCGACAAc                                 >  2:76890/1‑135 (MQ=255)
                                                                                                                        cTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATGGCTGCAGATCGACAAc                                 <  1:76890/135‑1 (MQ=255)
                                                                                                                        cTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATGGCTGCAGATCGACAACGTGTTCTATGAAGc                   >  4:121382/1‑149 (MQ=255)
                                                                                                                           tCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATGGCCGCAGATCGACAACGTGTTCTATGAAGCCCg                >  3:364157/1‑149 (MQ=255)
                                                                                                                             gCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATGGCTGCAGATCGACAACGTGTTCTATGAAGCCCGTg              >  1:184553/1‑149 (MQ=255)
                                                                                                                              cAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATGGCTGCAGATCGACAACGTGTTCTATGAAGCCCGTga             >  7:296938/1‑148 (MQ=255)
                                                                                                                               aGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATGGCTGCAGATCGACAACGTGTTCTATGAAGCCCGTGGc            <  5:353070/149‑1 (MQ=255)
                                                                                                                                    cTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATGGCTGCAGATCGACAACGTGTTCTATGAAGCCCGTGGCCAAGc       <  4:133914/149‑1 (MQ=255)
                                                                                                                                    cTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATGGCTGCAGATCGACAACGTGTTCTATGAAGCCCGTGGCCAAGc       <  6:242687/149‑1 (MQ=255)
                                                                                                                                    cTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATGGCTGCAGATCGACAACGTGTTCTATGAAGCCCGTGGCCAAGc       >  8:63373/1‑149 (MQ=255)
                                                                                                                                      cccGGCCAACG‑‑‑CTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATGGCTGCAGATCGACAAc                                 >  8:339631/3‑121 (MQ=255)
                                                                                                                                      cccGGCCAACG‑‑‑CTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATGGCTGCAGATCGACAAc                                 <  7:339631/119‑1 (MQ=255)
                                                                                                                                       tCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATGGCTGCAGATCGACAACGTGTTCTATGAAGCCCGTGGCCAAGcc      <  8:231811/147‑1 (MQ=255)
                                                                                                                                       ccGGCCAACG‑‑‑CTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATGGCTGCAGATCGACAACGTGTTCTATGAAGCCCGTGGCCAAGCCTg    >  1:286591/2‑149 (MQ=255)
                                                                                                                                         ggCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATGGCTGCAGATCGACAACGTGTTCTATGAAGCCCGTGGCCAAGCCTggg  <  7:264334/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
ATCAAGTCGCTGACCCGCTACCAGACCCGTCTGGCTGCCGGCGACAAGGGCGCTATCTTCTATACCCGCGCCGACAACCTGAACAACTGGCTGGGTGACGTCGCCACCCGTCTGGGCTCGCTGTCGCAGCGCCTGTCGGCCAGCG‑‑‑TTGGCCGGGTCAAGCTCAACACCACGCTGAAGACCGAGTCGGTCGTTGCCGGCCAGGCGCCGCAGGTGGATGAGGAGCTGGTGGAAACCCCATGGCTGCAGATCGACAACGTGTTCTATGAAGCCCGTGGCCAAGCCTGGG  >  NC_002947/657124‑657409

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