Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,096,109 G→A 5.7% R102H (CGC→CAC)  PP_2704 → M20/M25/M40 family peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,096,1090GA 5.7% 260.6 / 4.4 87R102H (CGC→CAC) PP_2704M20/M25/M40 family peptidase
Reads supporting (aligned to +/- strand):  ref base G (46/36);  new base A (2/3);  total (48/39)
Fisher's exact test for biased strand distribution p-value = 6.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

TGCATCGCAATGTGGGCGGGCATGGTGTGGTCGGCGTGCTGAAAGTGGGGGACAGCCCGCGCAGCATCGGTTTGCGTGCAGACATGGACGCACTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTACACGCCAGTTCAACGGCACGTTGAACGTGATATTCCAGCCCGACGAGGAAGGCCTGGCGGGGGCCAAGGCGATGATCGAC  >  NC_002947/3095965‑3096257
                                                                                                                                                |                                                                                                                                                    
tGCATCGCAATGTGGGCGGGCATGGTGTGGTCGGCGGGCTGAAAGTGGGGGACAGCCCGCGCAGCATCGGTTTGCGTGCAGACATGGACGCACTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCgcatgc                                                                                                                                                 >  2:108430/1‑150 (MQ=255)
  cATCGCAATGTGGGCGGGCATGGTGTGGTCGGCGTGCTGAAAGTGGGGGACAGCCCGCGCAGCATCGGTTTGCGTGCAGACATGGACGCACTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCgcatgcat                                                                                                                                               >  1:595575/1‑150 (MQ=255)
     cGCAATGTGGGCGGGCATGGTGTGGACTGCGTGCTCAACGTGGGGGACAGCCCGCCCAGCATTGGTTTGCGTGCAGACATGGACGCACTGCCCAGGGGCGAAAAGAATGCGATTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGcc                                                                                                                                            >  1:1655759/1‑150 (MQ=255)
      gCAATGTGGGCGGGCATGGTGTGGTCGGCGTGCTGAAAGTGGGGGACAGCCCGCGCAGCATCGGTTTGCGTGCAGACATGGACGCACTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCt                                                                                                                                           <  2:2056733/150‑1 (MQ=255)
       cAATGTGGGCGGGCATGGTGTGGTCGGCGTGCTGAAAGTGGGGGACAGCCCGCGCAGCATCGGTTTGCGTGCAGACATGGACGCACTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCACaggcatgcatg                                                                                                                                          >  1:197665/1‑147 (MQ=255)
          tgtgGGCGGGCATGGTGTGGTCGGCGTGCTGAAAGTGGGGGACAGCACGCGCAGCATCGGTTTGCGTGCAGACATGGACGCACTGCCCATGGGCGAAAAGAATGCGTTTGCGTATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTgg                                                                                                                                       <  1:1951988/150‑1 (MQ=255)
            tgGGCGGGCATGGTGTGGTCGGCGTGCTGAAAGTGGGGGACAGCCCGCGCAGCATCGGTTTGCGTGCAGACATGGACGCACTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCACAGGCATGCATGCGGcc                                                                                                                                     >  2:2030550/1‑150 (MQ=255)
              ggcgggcATGGTGTGGTCGGCGTGCTGAAAGTGGGGGACAGCCCGCGCAGCATCGGTTTGCGTGCAGACATGGACGCACTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCAt                                                                                                                                   <  2:536618/150‑1 (MQ=255)
              ggcgggcATGGTGTGGTCGGCGTGCTGAAAGTGGGGGACAGCCCGCGCAGCATCGGTTTGCGTGCAGACATGGACGCACTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCACAGGCATGCATGCGGccgg                                                                                                                                   <  1:1068122/150‑3 (MQ=255)
              ggcgggcATGGTGTGGTCGGCGTGCTGAAAGTGGGGGACAGCCCGCGCAGCATCGGTTTGCGTGCAGACATGGACGCACTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCACAGGCATGCATGCGGccgg                                                                                                                                   <  2:197665/150‑3 (MQ=255)
                 gggcATGGTGTGGTCGGCGTGCTGAAAGTGGGGGACAGCCCGCGCAGCATCGGTTTGCGTGCAGACATGGACGCACTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCACAGGCATGCATGCGGccgggga                                                                                                                                <  1:2030550/150‑6 (MQ=255)
                    cATGGTGTGGTCGGCGTGCTGAAAGTGGGGGACAGCCCGCGCAGCATCGGTTTGCGTGCAGACATGGACGCACTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACggg                                                                                                                             >  1:965940/1‑150 (MQ=255)
                     aTGGTGTGGTCGGCGTGCTGAAAGTGGGGGACAGCCCGCGCAGCATCGGTTTGCGTGCAGACATGGACGCACTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGc                                                                                                                            >  1:1827254/1‑150 (MQ=255)
                          gtgGTCGGCGTGCTGAAAGTGGGGGACAGCCCGCGCAGCATCGGTTTGCGTGCAGACATGGACGCACTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATAcc                                                                                                                       >  1:1590720/1‑150 (MQ=255)
                                     gCTGAAAGTGGGGGACAGCCCGCGCAGCATCGGTTTGCGTGCAGACATGGACGCACTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCctgct                                                                                                            <  2:1827254/150‑1 (MQ=255)
                                     gCTGAAAGTGGGGGACAGCCCGCGCAGCATCGGTTTGCGTGCAGACATGGACGCACTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCACGCACTTATGACGGGCATACCGCGATCctgct                                                                                                            <  1:1269213/150‑1 (MQ=255)
                                      cTGAAAGTGGGGGACAGCCCGCGCAGCATCGGTTTGCGTGCAGACATGGACGCACTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACggg                                                                                                                             >  2:625352/1‑132 (MQ=255)
                                      cTGAAAGTGGGGGACAGCCCGCGCAGCATCGGTTTGCGTGCAGACATGGACGCACTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACggg                                                                                                                             <  1:625352/132‑1 (MQ=255)
                                       tGAAAGTGGGGGACAGCTCGCGCAGCATCGGTTTGCGTGCAGACATGGACGCACTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTgg                                                                                                          >  2:1169355/1‑150 (MQ=255)
                                       tGAAAGTGGGGGACAGCCCGCGCAGCATCGGTTTGCGTGCAGACATGGACGCACTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTgg                                                                                                          >  2:1577732/1‑150 (MQ=255)
                                       tGAAAGTGGGGGACAGCCCGCGCAGCATCGGTTTGCGTGCAGACATGGACGCACTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTgg                                                                                                          >  1:842633/1‑150 (MQ=255)
                                        gAAAGTGGGGGACAGCCCGCGCAGCATCGGTTTGCGTGCAGACATGGACGCACTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGc                                                                                                         >  2:60671/1‑150 (MQ=255)
                                                gggACAGCCCGCGCAGCATCGGTTTGCGTGCAGACATGGACGCACTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTg                                                                                                 >  1:341083/1‑150 (MQ=255)
                                                 ggACAGCCCGCGCAGCATCGGTTTGCGTGCAGACATGGACGCACTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTgc                                                                                                >  1:1238787/1‑150 (MQ=255)
                                                   aCAGCCCGCGCAGCATCGGTTTGCGTGCAGACATGGACGCACTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTgcgc                                                                                              >  2:771286/1‑150 (MQ=255)
                                                     aGCCCGCGCAGCATCGGTTTGCGTGCAGACATGGACGCACTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCaa                                                                                            >  2:312921/1‑150 (MQ=255)
                                                      gCCCGCGCAGCATCGGTTTGCGTGCAGACATGGACGCACTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCtgtg                                                                                                                                        >  1:2038117/1‑105 (MQ=255)
                                                      gCCCGCGCAGCATCGGTTTGCGTGCAGACATGGACGCACTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCtgtg                                                                                                                                        <  2:2038117/105‑1 (MQ=255)
                                                      gCCCGCGCAGCATCGGTTTGCGTGCAGACATGGACGCACTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCAACTAAGATCCCCGGCCGCATGCATGCGTGTGGCCGTGTGGGGCATACTGCGATCCTGCTGGCAGCCGCTGCGCAAc                                                                                           >  1:541233/1‑150 (MQ=255)
                                                       cccGCGCAGCATCGGTTTGCGTGCAGACATGGACGCACTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACt                                                                                          >  2:1336387/1‑150 (MQ=255)
                                                                 aTCGGTTTGCGTGCAGACATGGACGCACTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTaca                                                                                >  1:2192758/1‑150 (MQ=255)
                                                                  tCGGTTTGCGTGCAGACATGGACGCACTGCCCATGGGCGAAAAGCATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGACTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTacac                                                                               >  1:2129448/1‑150 (MQ=255)
                                                                  tCGGTTTGCGTGCAGACATGGACGCACTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTacac                                                                               >  2:163155/1‑150 (MQ=255)
                                                                  tCGGTTTGCGTGCAGACATGGACGCACTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTacac                                                                               >  2:2376908/1‑150 (MQ=255)
                                                                    ggTTTGCGTGCAGACATGGACGCACTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTACACGc                                                                             <  2:341083/150‑1 (MQ=255)
                                                                     gTTTGCGTGCAGACATGGACGCACTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTACACGcc                                                                            >  2:1373169/1‑150 (MQ=255)
                                                                           gTGCAGACATGGACGCACTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTACACGCCAGTTCa                                                                      >  1:1335648/1‑150 (MQ=255)
                                                                             gCAGACATGGACGCACTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTACACGCCAGTTCAAc                                                                    <  1:1336387/150‑1 (MQ=255)
                                                                                 aCATGGACGCACTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTACACGCCAGTTCAACGGCa                                                                <  2:1188335/150‑1 (MQ=255)
                                                                                   aTGGACGCACTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGGTCCCCGGCCGCATGCAATCCTGTGGCCATGACTGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTACACGCAAGTGCAACGGCACg                                                              >  1:1375058/1‑150 (MQ=255)
                                                                                   aTGGACGCACTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTACACGCCAGTTCAACGGCACg                                                              >  1:1716024/1‑150 (MQ=255)
                                                                                   aTGGACGCACTGCCCATGGGCGAAAAGAATGCGGTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTACACGCCAGTTCAACGGCACg                                                              >  2:1933942/1‑150 (MQ=255)
                                                                                    tGGACGCACTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTACACGCCAGTTCAACGGCACGt                                                             <  2:842633/150‑1 (MQ=255)
                                                                                     ggACGCACTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTACACGCCAGTTCAACGGCACGtt                                                            >  1:219940/1‑150 (MQ=255)
                                                                                      gACGCACTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTACACGc                                                                             >  1:2189613/1‑132 (MQ=255)
                                                                                      gACGCACTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTACACGc                                                                             <  2:2189613/132‑1 (MQ=255)
                                                                                       aCGCACTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTACACGCCAGTTCAACGGCACGTTGa                                                          >  2:2032356/1‑150 (MQ=255)
                                                                                         gCACTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCatgcatgc                                                                                                                                             >  1:627605/1‑65 (MQ=255)
                                                                                         gCACTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCatgcatgc                                                                                                                                             <  2:627605/65‑1 (MQ=255)
                                                                                         gCACTGCCCATGGCCGAAAAGAATGTGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATTCATGCCTGTGTCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTACACGCCAGTTCAACGTGCCGTTGAAc                                                        <  2:721074/150‑1 (MQ=255)
                                                                                           aCTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTACACGCCAGTTCAAc                                                                    >  1:262028/1‑136 (MQ=255)
                                                                                           aCTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTACACGCCAGTTCAAc                                                                    <  2:262028/136‑1 (MQ=255)
                                                                                           aCTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTACACGCCAGTTCAACGTGCCGTTGAActg                                                      <  1:1455661/150‑3 (MQ=255)
                                                                                           aCTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTACACGCCAGTTCAACGGCACGTTGAACGt                                                      >  2:1170090/1‑150 (MQ=255)
                                                                                           aCTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTACACGCCAGTTCAACGGCACGTTGAACGt                                                      >  2:1294833/1‑150 (MQ=255)
                                                                                            cTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTACACGCCAGTTCAACGGCACGTTGAACGTg                                                     <  1:1170090/150‑1 (MQ=255)
                                                                                                ccATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTACACGCCAGTTCAACGGCACGTTGAACGTGatat                                                 >  2:1660211/1‑150 (MQ=255)
                                                                                                 cATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTggcgg                                                                                     >  2:117794/1‑113 (MQ=255)
                                                                                                 cATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTggcgg                                                                                     <  1:117794/113‑1 (MQ=255)
                                                                                                  aTGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTACACGCCAGTTCAACGGCACGTTGAACGTGATATTc                                               >  2:2239703/1‑150 (MQ=255)
                                                                                                    gggCGAAAAGAATGCGTTTGCTCATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTACACGCCAGTTCAACGGCACTTTGAACGTGATATTCCa                                             <  2:1238787/150‑1 (MQ=255)
                                                                                                      gCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTACACGCCAGTTCAACGGCACGTTGAACGTGATATTCCAGc                                           <  1:2032356/150‑1 (MQ=255)
                                                                                                        gAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCctgctg                                                                                                           >  1:1676112/1‑84 (MQ=255)
                                                                                                        gAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCctgctg                                                                                                           <  2:1676112/84‑1 (MQ=255)
                                                                                                         aaaaGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTACACGCCAGTTCAACGGCACGTTGAACGTGATATTCCAGCCcg                                        >  2:575736/1‑150 (MQ=255)
                                                                                                          aaaGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACTGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTACACGCCAGTTCAACGGCACGTTGAAc                                                        <  2:2144315/133‑1 (MQ=255)
                                                                                                          aaaGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACTGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTACACGCCAGTTCAACGGCACGTTGAAc                                                        >  1:2144315/1‑133 (MQ=255)
                                                                                                               aTGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTACACGCCAGTTCAACGGCACGTTGAACGTGATATTCCAGCCCGACGAgg                                  <  1:312921/150‑1 (MQ=255)
                                                                                                                         gcATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTACACGCCAGTTCAACGGCACGTTGAACGTGata                                                  >  1:167159/1‑124 (MQ=255)
                                                                                                                         gcATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTACACGCCAGTTCAACGGCACGTTGAACGTGata                                                  <  2:167159/124‑1 (MQ=255)
                                                                                                                         gcATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTACACGCCAGTTCAACGGCACGTTGAACGTGATATTCCAGCCCGACGAGGAAGGCCTGGc                        >  2:1439039/1‑150 (MQ=255)
                                                                                                                            tGCCAGCAAGACCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTACACGCCAGTTCAACGGCACGTTGAACGTGATATTCCAGCCCGACGAGGAAGGCCTGGCggg                     <  1:436093/150‑1 (MQ=255)
                                                                                                                             gCCAGCAAGGTCCCGGGCCGCATGCATCCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTACACGCCAGTTCAACGGCACGTTGAACGTGAGTTTCCAGCCCGACGCGGAAGGCCTGGCgggg                    <  1:2266881/150‑1 (MQ=255)
                                                                                                                              ccAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTACACGCCAGTTCAACGGCACGTTGAACGTGATATTCCAGCCCGACGAGGAAGGCCTGGCggggg                   <  2:1335648/150‑1 (MQ=255)
                                                                                                                               cAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTACACGCCAGATCAACGGCACGTTGAACGGGATATTCCAGCCCGCCGAGGAAGGCCTGGCGGGGGc                  >  1:827890/1‑150 (MQ=255)
                                                                                                                                aGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTACACGCCAGTTCAACGGCACGTTGAACGTGATATTCCAGCCCGACGAGGAAGGCCTGGCGGGGGcc                 >  2:1608414/1‑150 (MQ=255)
                                                                                                                                  cAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTACACGCCAGTTCAACGGCACGTTGAACGTGATATTCCAGCCCGACGAGGAAGGCCTGGCGGGGGCCaa               >  2:98961/1‑150 (MQ=255)
                                                                                                                                    aGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTACACGCCAGTTCAACGGCACGTTGAACGTGATATTCCAGCCCGACGAGGAAGGCCTGGCGGGGGCCAAgg             <  1:1169355/150‑1 (MQ=255)
                                                                                                                                     gATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTACACGCCAGTTCAACGGCACGTTGAACGTGATATTCCAGCCCGACGAGGAAGGCCTGGCGGGGGCCAAGGc            >  2:1464146/1‑150 (MQ=255)
                                                                                                                                       tCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTACACGCCAGTTCAACGGCACGTTGAACGTGATATTCCAGCCCGACGAGGAAGGCCTGGCGGGGGCCAAGGCga          <  2:1590720/150‑1 (MQ=255)
                                                                                                                                        ccccGGCCGCATGCATGCCTGTGGCCATGACGGGCATAcc                                                                                                                       <  2:414235/40‑1 (MQ=255)
                                                                                                                                        ccccGGCCGCATGCATGCCTGTGGCCATGACGGGCATAcc                                                                                                                       >  1:414235/1‑40 (MQ=255)
                                                                                                                                         cccGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTACACGCCAGTTCAACGGCACGTTGAACGTGATATTCCAGCCCGACGAGGAAGGCCTGGCGGGGGCCAAGGCgatg        <  1:2376908/150‑1 (MQ=255)
                                                                                                                                          ccGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTACACGCCAGTTCAACGGCACGTTGAACGTGATATTCCAGCCCGACGAGGAAGGCCTGGCGGGGGCCAAGGCga          >  2:1672244/1‑147 (MQ=255)
                                                                                                                                          ccGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTACACGCCAGTTCAACGGCACGTTGAACGTGATATTCCAGCCCGACGAGGAAGGCCTGGCGGGGGCCAAGGCga          <  1:1672244/147‑1 (MQ=255)
                                                                                                                                            ggCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTACACGCCAGTTCAACGGCACGTTGAACGTGATATTCCAGCCCGACGAGGAAGGCCTGGCGGGGGCCAAGGCGATGATc     <  2:1334943/150‑1 (MQ=255)
                                                                                                                                            ggCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTACACGCCAGTTCAACGGCACGTTGAACGTGATATTCCAGCCCGACGAGGAAGGCCTGGCGGGGGCCAAGGCGATGATc     <  1:163155/150‑1 (MQ=255)
                                                                                                                                             gCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTACACGCCAGTTCAACGGCACGTTGAACGTGATATTCCAGCCCGACGAGGAAGGCCTGGCGGGGGCCAAGGCGATGATcg    <  1:1439039/150‑1 (MQ=255)
                                                                                                                                               cGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTACACGCCAGTTCAACGGCACGTTGAACGTGATATTCCAGCCCGACGAGGAAGGCCTGGCGGGGGCCAAGGCGATGATcgac  <  2:2092864/150‑1 (MQ=255)
                                                                                                                                                gcatgcatGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTACACGCCAGTTCAACGGCACGTTGAACGTGa                                                    <  2:1646605/99‑1 (MQ=255)
                                                                                                                                                gcatgcatGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTACACGCCAGTTCAACGGCACGTTGAACGTGa                                                    >  1:1646605/1‑99 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
TGCATCGCAATGTGGGCGGGCATGGTGTGGTCGGCGTGCTGAAAGTGGGGGACAGCCCGCGCAGCATCGGTTTGCGTGCAGACATGGACGCACTGCCCATGGGCGAAAAGAATGCGTTTGCGCATGCCAGCAAGATCCCCGGCCGCATGCATGCCTGTGGCCATGACGGGCATACCGCGATCCTGCTGGCAGCCGCTGCGCAACTGGCGGCTACACGCCAGTTCAACGGCACGTTGAACGTGATATTCCAGCCCGACGAGGAAGGCCTGGCGGGGGCCAAGGCGATGATCGAC  >  NC_002947/3095965‑3096257

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