Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,000,466 G→T 5.3% V122F (GTT→TTT)  PP_1786 → glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,000,4660GT 5.3% 240.9 / 2.5 76V122F (GTT→TTT) PP_1786glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base G (36/36);  new base T (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 = 6.88e-01

CCTCGGCTTCGCACGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGAACGGGTTCAGCCTTTCAGCA  >  NC_002947/2000317‑2000603
                                                                                                                                                       |                                                                                                                                          
ccTCGGCTTCGCACGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGgcggcg                                                                                                                                            >  2:1570366/1‑150 (MQ=255)
 cTCGGCTTCGCACGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGt                                                                                                                                           <  2:826369/150‑1 (MQ=255)
 cTCGGCTTCGCACGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGt                                                                                                                                           <  2:500057/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                                                                                                                                   <  1:766755/147‑1 (MQ=255)
            aCGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTc                                                                                                                                   >  2:766755/1‑147 (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:1538126/150‑1 (MQ=255)
                       gATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGaa                                                                                                                     >  1:199415/1‑150 (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:1447537/150‑1 (MQ=255)
                                                                      cGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTActct                                                                     <  1:1406723/150‑1 (MQ=255)
                                                                          aTGGTCAAGCCGGATGTTGCTAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTAt                                                                                                                            >  2:1510617/1‑92 (MQ=255)
                                                                          aTGGTCAAGCCGGATGTTGCTAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTAt                                                                                                                            <  1:1510617/92‑1 (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                                                                                                                                      <  1:1214106/52‑1 (MQ=255)
                                                                                                        cTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGaaa                                                                                                                                      >  2:1214106/1‑52 (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:215827/1‑150 (MQ=255)
                                                                                                                     tGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGa                      >  1:1559330/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                  <  1:1168505/136‑1 (MQ=255)
                                                                                                                                       tCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGAACgg                  >  2:1168505/1‑136 (MQ=255)
                                                                                                                                       tCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGAACGGGTTCAGc           <  2:1009591/143‑1 (MQ=255)
                                                                                                                                       tCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGAACGGGTTCAGc           >  1:1009591/1‑143 (MQ=255)
                                                                                                                                        cAAAG‑‑CCGGCGGTGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTTGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGAACGGGTTCAGCCTTTcagc   <  1:827683/150‑1 (MQ=255)
                                                                                                                                         aaaG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGAACGGGTTCAGCCTTTcagca  <  1:1753365/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                          
CCTCGGCTTCGCACGAAAGTGGTGATGCGGAGCTTATCGCTGAAGCAGTTCTCAATGCGGCCGACAGCGCCGTTATGGTCAAGCCGGATGTTGCGAAAGTTGACCTGGATTTTTACCTGACCTTTTACCCTGACCTCAAAG‑‑CCGGCGGCGTTGAAACTCAAGCTATTGCGGAAGC‑ACATTTCCAACGCTTCGGACGGGCGGAGGGGCGACTTCCTACTCTCGATGCCTGGGCTAAGAAGCACGAGCTTCCACTAGAGCTGTTACCGAACGGGTTCAGCCTTTCAGCA  >  NC_002947/2000317‑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: ‑