Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,000,461 G→C 5.3% G120A (GGC→GCC)  PP_1786 → glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,000,4610GC 5.3% 233.8 / 4.4 76G120A (GGC→GCC) PP_1786glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base G (35/36);  new base C (2/2);  total (38/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.13e-01

TGCCTCGGCTTCGCACGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGAACGGGTTCAGCCTTTCAGCA  >  NC_002947/2000315‑2000603
                                                                                                                                                    |                                                                                                                                               
tGCCTCGGCCTCGCACGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑Ccggcgg                                                                                                                                              <  1:763320/150‑1 (MQ=255)
  ccTCGGCTTCGCACGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGgcggcg                                                                                                                                            >  2:1570366/1‑150 (MQ=255)
   cTCGGCTTCGCACGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGt                                                                                                                                           <  2:500057/150‑1 (MQ=255)
   cTCGGCTTCGCACGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGt                                                                                                                                           <  2:826369/150‑1 (MQ=255)
            gCACGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGCAAGTTGACGTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCa                                                                                                                                  >  2:1698416/1‑150 (MQ=255)
            gCACGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCa                                                                                                                                  >  2:389117/1‑150 (MQ=255)
            gCACGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCGTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCa                                                                                                                                  >  1:1258612/1‑150 (MQ=255)
             cACGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTGTATGCGGCCGGCAGCGGCGTTATGGTCGAGCCGGATGTTGCGACAGTTGACCTGGATTTTTGCTTGTACCTTTAGCCTGCCGTCAACG‑‑CTGGCGGCGTTGAAACTCag                                                                                                                                 <  2:202932/150‑2 (MQ=255)
             cACGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCaa                                                                                                                                 >  2:250758/1‑150 (MQ=255)
              aCGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAACG‑‑CCGCCGGCTTTGaggtcag                                                                                                                                  >  1:1082375/1‑142 (MQ=255)
              aCGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAACG‑‑CCGCCGGCTTTGaggtcag                                                                                                                                  <  2:1082375/148‑7 (MQ=255)
              aCGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTc                                                                                                                                   >  2:766755/1‑147 (MQ=255)
              aCGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTc                                                                                                                                   <  1:766755/147‑1 (MQ=255)
              aCGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAg                                                                                                                                >  1:1543641/1‑150 (MQ=255)
                gAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCt                                                                                                                              >  2:804844/1‑150 (MQ=255)
                gAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTGTACCCTGACCTCAAAT‑‑CCGGCTGCGTTGAAACTCAAGCt                                                                                                                              >  2:1753365/1‑150 (MQ=255)
                     tggtgATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGc                                                                                                                         <  1:1301535/150‑1 (MQ=255)
                         gATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGaa                                                                                                                     >  1:199415/1‑150 (MQ=255)
                         gATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGaa                                                                                                                     <  1:1538126/150‑1 (MQ=255)
                          aTGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGTAAg                                                                                                                    <  1:804844/150‑1 (MQ=255)
                          aTGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCCACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAg                                                                                                                    <  1:577765/150‑1 (MQ=255)
                              ggAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGGTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGc‑aca                                                                                                               >  2:896432/1‑150 (MQ=255)
                                  cTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTc                                                                                                           >  1:1063748/1‑150 (MQ=255)
                                     aTCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCaa                                                                                                        <  2:630573/150‑1 (MQ=255)
                                     aTCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCaa                                                                                                        <  1:948283/150‑1 (MQ=255)
                                          tGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAACCTTTGGCGGAAGC‑ACATTTCCAACGCtt                                                                                                   <  1:1698416/150‑1 (MQ=255)
                                              gCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGa                                                                                               >  1:303458/1‑150 (MQ=255)
                                               cAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAACG‑‑CCGCCGGCTTTGag                                                                                                                                       <  2:1095362/110‑2 (MQ=255)
                                               cAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAACG‑‑CCGCCGGCTTTGag                                                                                                                                       >  1:1095362/1‑109 (MQ=255)
                                                  ttCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACggg                                                                                           <  2:424492/150‑1 (MQ=255)
                                                     tcAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCgg                                                                                        >  2:1447537/1‑150 (MQ=255)
                                                              ccGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGAc                                                                               >  1:1053103/1‑150 (MQ=255)
                                                                gACAGCGCCGTTATGGTCAAGCCGGATGTTGCGACAGGTTAGCTGGAGTGTTACCTGACCCTTTACCCTGAGATAAAAGTACCTTTACCGTT‑‑ATCTGAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGAAGGGCGGGGGGGCGACtt                                                                             >  1:896375/1‑150 (MQ=255)
                                                                       ccGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTActc                                                                      >  1:383277/1‑150 (MQ=255)
                                                                        cGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTActct                                                                     <  1:1406723/150‑1 (MQ=255)
                                                                        cGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTActct                                                                     <  1:1447537/150‑1 (MQ=255)
                                                                            aTGGTCAAGCCGGATGTTGCTAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTAt                                                                                                                            <  1:1510617/92‑1 (MQ=255)
                                                                            aTGGTCAAGCCGGATGTTGCTAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTAt                                                                                                                            >  2:1510617/1‑92 (MQ=255)
                                                                                tCAAGCCGGATGTTGCGACAGTTGGCCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCACGCTATTGCGGACGC‑ACACTTCCAACGCTTTGGACGGGCGGAGGGGCGAGTGCCTACTCTCGTTGCCt                                                             >  2:1648582/1‑150 (MQ=255)
                                                                                     ccGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCt                                                        <  2:896375/150‑1 (MQ=255)
                                                                                         aTGTTGCGAAAGTTGACGTGGATTTTTACCTGACCTTTTACCCTGTCCTCAAAG‑‑CCGGCGGAGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGTTaaga                                                    <  2:759681/150‑1 (MQ=255)
                                                                                          tGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTaagaa                                                   <  1:250758/150‑1 (MQ=255)
                                                                                              gCGAAAGTTGACCTGGATTTTTACCTGACCCTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCAc                                               <  2:1258612/150‑1 (MQ=255)
                                                                                               cGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCg                                                                                                 >  1:39661/1‑99 (MQ=255)
                                                                                               cGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCg                                                                                                 <  2:39661/99‑1 (MQ=255)
                                                                                                gAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTCGTGTAGC‑ACATTGCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGa                                             >  1:1670465/1‑150 (MQ=255)
                                                                                                gAAAGTTAACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGa                                             <  2:300867/150‑1 (MQ=255)
                                                                                                     ttGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTc                                        <  2:1063748/150‑1 (MQ=255)
                                                                                                          cTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGaaa                                                                                                                                      >  2:1214106/1‑52 (MQ=255)
                                                                                                          cTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGaaa                                                                                                                                      <  1:1214106/52‑1 (MQ=255)
                                                                                                           tGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGc                                           <  1:344216/141‑1 (MQ=255)
                                                                                                           tGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGc                                           >  2:344216/1‑141 (MQ=255)
                                                                                                            ggATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTaga                                 <  1:389117/150‑1 (MQ=255)
                                                                                                            ggATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTaga                                 <  2:1663441/150‑1 (MQ=255)
                                                                                                               tttttCCCTGACCTTCTGCCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAACCTCGT‑CGGCTTCGCCATCTGCAACTCGTCGGACGGTCGGAGCGGCGCCTCCCTACTATCGGTTGATGGGTTAAGAAGCACGAGCTTCCACTAGAGCt                              <  1:1876933/150‑1 (MQ=255)
                                                                                                                     ccTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAGACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTAcc                        >  1:1256951/1‑150 (MQ=255)
                                                                                                                      cTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCg                       >  1:1913495/1‑150 (MQ=255)
                                                                                                                       tGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGa                      >  1:1559330/1‑150 (MQ=255)
                                                                                                                       tGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGa                      >  1:215827/1‑150 (MQ=255)
                                                                                                                        gACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCAc                                     >  1:1293995/1‑134 (MQ=255)
                                                                                                                        gACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCAc                                     <  2:1293995/134‑1 (MQ=255)
                                                                                                                          ccTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATCTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGAACg                   >  2:299367/1‑150 (MQ=255)
                                                                                                                           cgtttACCCGGACCTCAAAT‑‑CCGGCGGCGGGGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGAACgg                  >  2:504209/3‑150 (MQ=255)
                                                                                                                            tgttACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGAACggg                 >  2:1268431/3‑150 (MQ=255)
                                                                                                                                 cccTAACCTACAAG‑‑CCGGGGGCGGTGAAACGCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGAACGGGTTCAg            <  2:1543641/150‑1 (MQ=255)
                                                                                                                                  ccTGCCCTCATAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGAGCGTTCCTACTCTCGATGCCTGGGCTCAGACGCACTAGCTTCCACTAGAGCTGTTACCGAACGGGTTCAGc           <  1:1570366/150‑1 (MQ=255)
                                                                                                                                  ccTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACTTTTCGAACGCGTCGGACGTTCTCAGGGGCAAAATCCTGCTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGAACTGGTTCAGc           >  1:1875333/1‑150 (MQ=255)
                                                                                                                                  ccTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGAACGGGtgcggc           >  1:1765051/1‑145 (MQ=255)
                                                                                                                                   cTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGAACGGGTTCAGcc          >  2:1474199/1‑150 (MQ=255)
                                                                                                                                    tGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGAACGGGTTCAGCCt         <  2:1913495/150‑1 (MQ=255)
                                                                                                                                        cTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTa                          >  2:317317/1‑129 (MQ=255)
                                                                                                                                        cTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTa                          <  1:317317/129‑1 (MQ=255)
                                                                                                                                         tCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGAACgg                  >  2:1168505/1‑136 (MQ=255)
                                                                                                                                         tCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGAACgg                  <  1:1168505/136‑1 (MQ=255)
                                                                                                                                         tCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGAACGGGTTCAGc           >  1:1009591/1‑143 (MQ=255)
                                                                                                                                         tCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGAACGGGTTCAGc           <  2:1009591/143‑1 (MQ=255)
                                                                                                                                          cAAAG‑‑CCGGCGGTGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTTGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGAACGGGTTCAGCCTTTcagc   <  1:827683/150‑1 (MQ=255)
                                                                                                                                           aaaG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGAACGGGTTCAGCCTTTcagca  <  1:1753365/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
TGCCTCGGCTTCGCACGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGAACGGGTTCAGCCTTTCAGCA  >  NC_002947/2000315‑2000603

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑