Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,791,876 G→C 7.3% G162R (GGT→CGT)  rseP → regulatory intramembrane protein RIP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,791,8760GC 7.3% 302.3 / 12.7 98G162R (GGT→CGT) rsePregulatory intramembrane protein RIP
Reads supporting (aligned to +/- strand):  ref base G (44/45);  new base C (3/4);  total (47/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.52e-01

TGGCTATCCTGTTTTTCTGGGTGCTGGCGATGCTGGGTACCCAGCAAATCCGCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAGGTCAAGCTGGACAGTTGGCTCAAGGGCGCCGACGAGCCGGAC  >  NC_002947/1791727‑1792022
                                                                                                                                                     |                                                                                                                                                  
tgGCTATCCTGTTTTTCTGGGTGCTGGCGATGCTGGGTACCCAGCAAATCCGCCCGGTTATCGGAGCGGTCGACAGTGGCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACg                                                                                                                                                    >  1:982337/1‑150 (MQ=255)
 ggCTATCCTGTTTTTCTGGGTGCTGGCGATGCTGGGCACCCAGCAAATCCGCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACgg                                                                                                                                                   >  2:2184344/1‑150 (MQ=255)
  gCTATCCTGTTTTTCTGGGTGCTGGCGATGCTGGGTACCCAGCAAATCCGCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGt                                                                                                                                                  <  2:1894167/150‑1 (MQ=255)
      tCCTGTTTTTCTGGGTGCTGGCGATGCTGGGTACCCAGCAAATCCGCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGt                                                                                                                                              <  2:2219281/150‑1 (MQ=255)
         tGTTTTTCTGGGTGCTGGCGATGCTGGGTACCCAGCAAATCCGCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTcgg                                                                                                                                           >  1:1923751/1‑150 (MQ=255)
         tGTTTTTCTGGGTGCTGGCGATGCTGGGTACCCAGCAAATCCGCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGACGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCCAGCCGACCAACGGTTGGTcgg                                                                                                                                           >  1:1135186/1‑150 (MQ=255)
          gTTTTTCTGGGTGCTGGCGATGCTGGGTACCCAGCAAATCCGCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTcggc                                                                                                                                          <  1:1060766/150‑1 (MQ=255)
           tttttCTGGGTGCTGGCGATGCTGGGTACCCAGCAAATCCGCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTcggcg                                                                                                                                         <  1:692437/150‑1 (MQ=255)
           tttttCTGGGTGCTGGCGATGCTGGGTACCCAGCAAATCCGCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTcggcg                                                                                                                                         >  1:950979/1‑150 (MQ=255)
           tttttCTGGGTGCTGGCGATGCTGGGTACCCAGCAAATCCGCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTcggcg                                                                                                                                         <  2:1506533/150‑1 (MQ=255)
            ttttCTGGGTGCTGGCGATGCTGGGTACCCAGCAAATCCGCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAGCATCTGCAGGACTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACCGTTGGTcggctt                                                                                                                                        <  1:1626985/150‑3 (MQ=255)
             tttCTGGGTGCTGGCGATGCTGGGTACCCAGCAAATCCGCCCGGTTATCGGCGCGCTCGACAGTGGCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACCGTTGGTcggcttg                                                                                                                                       <  1:728669/150‑4 (MQ=255)
                  gggTGCTGGCGATGCTGGGTACCCAGCAAATCCGCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAAtt                                                                                                                                  >  1:689606/1‑150 (MQ=255)
                      gCTGGCGATGCTGGGTACCCAGCAAATCCGCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGTTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCa                                                                                                                              <  2:400997/150‑1 (MQ=255)
                          gCGATGCTGGGTACCCAGCAAATCCGCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTcggc                                                                                                                                          <  1:118636/134‑1 (MQ=255)
                          gCGATGCTGGGTACCCAGCAAATCCGCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTcggc                                                                                                                                          >  2:118636/1‑134 (MQ=255)
                           cGATGCTGGGTACCCAGCAAATCCGCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTgg                                                                                                                         >  2:1320021/1‑150 (MQ=255)
                           atttGCTGGGTACCCAGCAAATCCGCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGtt                                                                                                                                                 >  2:939670/4‑126 (MQ=255)
                           atttGCTGGGTACCCAGCAAATCCGCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGtt                                                                                                                                                 <  1:939670/123‑1 (MQ=255)
                             aTGCTGGGTACCCAGCAAATCCGCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGtt                                                                                                                       >  2:1682350/1‑150 (MQ=255)
                                 tGGGTACCCAGCAAATCCGCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTtcgtc                                                                                                                   >  1:1760042/1‑150 (MQ=255)
                                       cccAGCAAATCCGCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTgg                                                                                                             <  1:2184344/150‑1 (MQ=255)
                                       cccAGCAAATCCGCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTgg                                                                                                             <  1:2356880/150‑1 (MQ=255)
                                       cccAGCAAATCCGCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTgg                                                                                                             <  2:950979/150‑1 (MQ=255)
                                         cAGCAAATCCGCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGc                                                                                                           >  2:721014/1‑150 (MQ=255)
                                         cAGCAAATCCGCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGc                                                                                                           >  2:1808348/1‑150 (MQ=255)
                                         cAGCAAATCCGCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGc                                                                                                           >  2:1992037/1‑150 (MQ=255)
                                          aGCAAATCCGCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCg                                                                                                          >  2:1520970/1‑150 (MQ=255)
                                           gCAAATCCGCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCga                                                                                                         <  2:1785161/150‑1 (MQ=255)
                                             aaaTCCGCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCgaga                                                                                                       >  1:890643/1‑150 (MQ=255)
                                                tCCGCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCg                                                                                                    >  2:768566/1‑150 (MQ=255)
                                                       ggTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCtt                                                                                             >  1:2395305/1‑150 (MQ=255)
                                                        gTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTg                                                                                            <  2:988833/150‑1 (MQ=255)
                                                            tCGGCGCGGTCGACAGTGGCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGa                                                                                        <  1:1682350/150‑1 (MQ=255)
                                                              ggCGCGGTCGACAGTGGCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGAtt                                                                                      <  2:689606/150‑1 (MQ=255)
                                                              ggCGCGGTCGACAGTGGCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACCGTTGGTcggc                                                                                                                                          <  2:202332/98‑1 (MQ=255)
                                                              ggCGCGGTCGACAGTGGCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACCGTTGGTcggc                                                                                                                                          >  1:202332/1‑98 (MQ=255)
                                                               gcgcGGTCGACAGTGGCAGCCTGGCAGCATCTGCAGGGGTGACGGCGGGTTAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACTGTTGGTTGGCGGTCAATTTGCATCTGGTTCGTCGTCTCTGCGAGAGCGGTACCTTGCAGATTg                                                                                     <  2:568242/150‑1 (MQ=255)
                                                                 gcgGTCGACAGTGGCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGc                                                                                   <  1:344789/150‑1 (MQ=255)
                                                                 gcgGTCGACAGTGGCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGc                                                                                   <  2:2291835/150‑1 (MQ=255)
                                                                 gcgGTCGACAGGGGCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGc                                                                                   <  1:1647440/150‑1 (MQ=255)
                                                                    gTCGACAGTGGCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAg                                                                                                                             >  1:944616/1‑105 (MQ=255)
                                                                    gTCGACAGTGGCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAg                                                                                                                             <  2:944616/105‑1 (MQ=255)
                                                                     gtacttgTGGCAGCCTGGCAGCATCTGCAGGGCTGACCACGGGTGAGGAAATCGCTTCTGGGGATGGCAATCCGACCAACGGTTGGTCGGCGGTCAACTTGCAGCTGGTTCGACGTCTGGGCGAGAGCGGTACATTCCAGATTGGCGTgg                                                                               >  2:170316/7‑149 (MQ=255)
                                                                      cGACAGTGGCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTgcg                                                                              <  2:2208251/150‑1 (MQ=255)
                                                                        aCAGTGGCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACCGTTGGTcgg                                                                                                                                           <  1:1420997/87‑1 (MQ=255)
                                                                        aCAGTGGCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACCGTTGGTcgg                                                                                                                                           >  2:1420997/1‑87 (MQ=255)
                                                                          aGTGGCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGAt                                                                          <  1:1520970/150‑1 (MQ=255)
                                                                            tGGCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGa                                                                        <  2:2247058/150‑1 (MQ=255)
                                                                             ggCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAg                                                                       <  1:1320021/150‑1 (MQ=255)
                                                                              gcaagcTGGGAGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGTCAGCCGACCAATTGTTCGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAgg                                                                      <  1:170316/145‑1 (MQ=255)
                                                                              gCAGCCTGGCAGCATCTGCAGTGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGc                                                                                   >  1:804531/1‑137 (MQ=255)
                                                                              gCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGTTCAGGTAATTGTTTCTGTAGATGGCAAGGCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCTGATTGGCGTGCGCGATgggg                                                                      <  2:400848/150‑1 (MQ=255)
                                                                              gCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGc                                                                                   <  2:804531/137‑1 (MQ=255)
                                                                              gCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGTgg                                                                      <  2:287471/150‑1 (MQ=255)
                                                                              gCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAgg                                                                      >  2:846805/1‑150 (MQ=255)
                                                                                      gcagcaTCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAgc                                                              <  2:1760042/150‑1 (MQ=255)
                                                                                         gcaTCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGcc                                                           <  1:721014/150‑1 (MQ=255)
                                                                                            tCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAg                                                        <  2:890643/150‑1 (MQ=255)
                                                                                              tGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGAAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAgcg                                                      >  1:1322523/1‑150 (MQ=255)
                                                                                               gCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAgcgc                                                     >  2:542934/1‑150 (MQ=255)
                                                                                                cAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGcc                                                    >  1:974548/1‑150 (MQ=255)
                                                                                                  gggCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATg                                                                         >  2:2396022/1‑127 (MQ=255)
                                                                                                  gggCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATg                                                                         <  1:2396022/127‑1 (MQ=255)
                                                                                                      tGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGc                                              >  2:1313104/1‑150 (MQ=255)
                                                                                                        aCCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAg                                            >  1:709538/1‑150 (MQ=255)
                                                                                                        aCCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAg                                            >  1:1080941/1‑150 (MQ=255)
                                                                                                           gcgGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGGGGGCGCCAGCGCCGAGCGCCAGCTGCAGGTc                                         >  1:13673/1‑150 (MQ=255)
                                                                                                               gTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAGGTCAAGc                                     <  2:1923751/150‑1 (MQ=255)
                                                                                                                   ggAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAGGTCAAGCTGGa                                 <  2:1080941/150‑1 (MQ=255)
                                                                                                                    gAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAGGTCAAGCTGGAc                                >  1:2146684/1‑150 (MQ=255)
                                                                                                                       aTTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCaat                                                              <  1:1049270/117‑3 (MQ=255)
                                                                                                                       aTTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCaat                                                              >  2:1049270/1‑115 (MQ=255)
                                                                                                                        ttGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAGGTCAAGCTGGACAGtt                            >  1:131373/1‑150 (MQ=255)
                                                                                                                              cTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAGGTCAAGCTGGACAGTTGGCTCa                      >  2:1617671/1‑150 (MQ=255)
                                                                                                                              cTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAGGTCAAGCTGGACAGTTGGCTCa                      >  1:1685273/1‑150 (MQ=255)
                                                                                                                                 tgGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTg                                                                                     >  1:1304252/1‑84 (MQ=255)
                                                                                                                                 tgGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTg                                                                                     <  2:1304252/84‑1 (MQ=255)
                                                                                                                                   gATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAGGTCAAGCTGGa                                 >  2:1970649/1‑134 (MQ=255)
                                                                                                                                   gATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAGGTCAAGCTGGa                                 <  1:1970649/134‑1 (MQ=255)
                                                                                                                                   gATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAGGTCAAGCTGGACAGTTGGCTCAAGGgc                 >  1:1383210/1‑150 (MQ=255)
                                                                                                                                   gATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAGGTCAAGCTGGACAGTTGGCTCAAGGgc                 <  1:846805/150‑1 (MQ=255)
                                                                                                                                   gATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAGGTCAAGCTGGACAGTTGGCTCAAGGgc                 >  2:847220/1‑150 (MQ=255)
                                                                                                                                   gATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAGGTCAAGCTGGACAGTTGGCTCAAGGgc                 <  1:907713/150‑1 (MQ=255)
                                                                                                                                      ggCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAGGTCAAGCTGGACAGTTg                           <  1:1721009/137‑1 (MQ=255)
                                                                                                                                      ggCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCAAGCTGCAGGTCAAGCTGGACAGTTg                           >  2:1721009/1‑137 (MQ=255)
                                                                                                                                       gCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCa                                                                                          >  1:2235371/1‑73 (MQ=255)
                                                                                                                                       gCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCa                                                                                          <  2:2235371/73‑1 (MQ=255)
                                                                                                                                        cAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCCGCAGGTCAAGCTGGACAGTTGGCTCAAGGGCGCcga            >  1:287077/1‑150 (MQ=255)
                                                                                                                                          aGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGATGGCGCCAGCGCCGAGCGCCAGCTGCAGGTCAAGCTGGACAGTTGGCTCAAGGGCGCcgacg          >  1:898007/1‑150 (MQ=255)
                                                                                                                                          aGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAGGTCAAGCTGGACAGTTGGCTCAAGGGCGccttcg          <  2:2146684/150‑5 (MQ=255)
                                                                                                                                           gCCGACCAACCGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAGGTCAAGCTGGACAGTTGGCTCAAGGGCGCcgacga         >  1:877454/1‑150 (MQ=255)
                                                                                                                                             cGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAGGTCAAGCTGGACAGTTGGCTCAAGGGCGCCGACGAGc       <  2:974548/150‑1 (MQ=255)
                                                                                                                                               acaaaCGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCCGCAGGTCAAGCTGGACAGTTGGCTCAAGGGCGCCGACGAGCCg     <  2:287077/147‑1 (MQ=255)
                                                                                                                                                 cAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAGGTCAAGCTGGACAGTTGGCTCAAGGGCGCCGACGAGCCGGa   >  1:866948/1‑150 (MQ=255)
                                                                                                                                                  aaCGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGc                                                 >  2:243256/1‑103 (MQ=255)
                                                                                                                                                  aaCGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGc                                                 <  1:243256/103‑1 (MQ=255)
                                                                                                                                                  aaCGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAGGTCAAGCTGGACAGTTGGCTCAAGGGCGCCGACGAGcc      <  2:1185839/146‑1 (MQ=255)
                                                                                                                                                  aaCGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAGGTCAAGCTGGACAGTTGGCTCAAGGGCGCCGACGAGcc      >  1:1185839/1‑146 (MQ=255)
                                                                                                                                                  aaCGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAGGTCAAGCTGGACAGTTGGCTCAAGGGCGCCGACGAGCCGGAc  <  2:709538/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TGGCTATCCTGTTTTTCTGGGTGCTGGCGATGCTGGGTACCCAGCAAATCCGCCCGGTTATCGGCGCGGTCGACAGTGGCAGCCTGGCAGCATCTGCAGGGCTGACCGCGGGTCAGGAAATTGTTTCTGTGGATGGCAAGCCGACCAACGGTTGGTCGGCGGTCAATTTGCAGCTGGTTCGTCGTCTGGGCGAGAGCGGTACCTTGCAGATTGGCGTGCGCGATGAGGGCGCCAGCGCCGAGCGCCAGCTGCAGGTCAAGCTGGACAGTTGGCTCAAGGGCGCCGACGAGCCGGAC  >  NC_002947/1791727‑1792022

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: