Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,509,458 A→T 9.0% S355C (AGC→TGC)  urtA → urea ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,509,4580AT 9.0% 186.2 / 5.8 66S355C (AGC→TGC) urtAurea ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (23/37);  new base T (2/4);  total (25/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAGGCCAAGGGCCTGCCAGGTGCCGACAAGGCCGTGACCAACGACCCCATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGACGGCCAGTTCAGCGTGGTGTGGGAAACCGAGCA  >  NC_002947/5509310‑5509582
                                                                                                                                                     |                                                                                                                            
caaggccaaggGCCTGCCAGGTGCCGACAAGGCCGTGACCAACGACCCCATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCaga                                                                                                                              <  5:15737/149‑1 (MQ=255)
caaggccaaggGCCTGCCAGGTGCCGACAAGGCCGTGACCAACGACCCCATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCaga                                                                                                                              <  1:240243/149‑1 (MQ=255)
caaggccaaggGCCTGCCAGGTGCCGACAAGGCCGTGACCAACGACCCCATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCaga                                                                                                                              <  3:97874/149‑1 (MQ=255)
      caaggGCCTGCCAGGTGCCGACAAGGCCGTGACCAACGACCCCATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCg                                                                                                                        <  2:297470/149‑1 (MQ=255)
       aaggGCCTGCCAGGTGCCGACAAGGCCGTGACCAACGACCCCATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCTCTGCCCGGCCAGTGCTTCgc                                                                                                                       <  2:188223/149‑2 (MQ=255)
       aaggGCCTGCCAGGTGCCGACAAGGCCGTGACCAACGACCCCATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGa                                                                                                                       <  7:30681/149‑1 (MQ=255)
        aggGCCTGCCAGGTGCCGACAAGGCCGTGACCAACGACCCCATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCTCTGCCCGGCCAGTGCTTCgcg                                                                                                                      <  5:97583/149‑3 (MQ=255)
        aggGCCTGCCAGGAGCCGACAAGGCCGTGACCAACGACCCCATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACTTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAg                                                                                                                      <  3:32935/149‑1 (MQ=255)
            ccTGCCAGGTGCCGACAAGGCCGTGACCAACGACCCCATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCtt                                                                                                                          >  8:104339/1‑141 (MQ=255)
            ccTGCCAGGTGCCGACAAGGCCGTGACCAACGACCCCATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCtt                                                                                                                          <  7:104339/141‑1 (MQ=255)
                     tGCCGACAAGGCCGTGACCAACGACCCCATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCt                                                                                                         <  3:52244/149‑1 (MQ=255)
                          aCAAGGCCGTGACCAACGACCCCATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAGGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCTCTGCCCGGCCAGTGCTTCgc                                                                                                                       <  6:263828/130‑2 (MQ=255)
                          aCAAGGCCGTGACAAACGACCCCATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAGGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCTCTGCCCGGCCAGTGCTTCgc                                                                                                                       >  5:263828/1‑129 (MQ=255)
                            aaGGCCGTGACCAACGACCCCATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCgg                                                                                                            <  8:214086/139‑1 (MQ=255)
                            aaGGCCGTGACCAACGACCCCATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCgg                                                                                                            >  7:214086/1‑139 (MQ=255)
                              ggCCGTGACCAACGACCCCATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGa                                                                                                <  2:181268/149‑1 (MQ=255)
                               gCCGTGACCAACGACCCCATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGAc                                                                                               >  7:134785/1‑149 (MQ=255)
                                     aCCAACGACCCCATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGa                                                                                         <  4:183097/149‑1 (MQ=255)
                                      ccAACGACCCCATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCTAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGCCTTCACCCTGACCATGGa                                                                                         >  5:310089/1‑148 (MQ=255)
                                         aCGACCCCATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGc                                                                                                          >  2:245366/1‑128 (MQ=255)
                                         aCGACCCCATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGc                                                                                                          <  1:245366/128‑1 (MQ=255)
                                          cGACCCCATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGGTTCACCCTGACCATGGACAAGa                                                                                    >  7:98827/1‑149 (MQ=255)
                                             ccccATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGacc                                                                                  >  4:301748/1‑148 (MQ=255)
                                              cccATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGaccaa                                                                                <  2:58441/149‑1 (MQ=255)
                                                cATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGaccaacc                                                                              <  4:106066/149‑1 (MQ=255)
                                                    gAAGCCACCTACGTGTGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAaccacca                                                                          <  3:272959/149‑1 (MQ=255)
                                                      aGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAaccaccacc                                                                        >  7:13900/1‑149 (MQ=255)
                                                        ccACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTg                                                                      <  7:376/149‑1 (MQ=255)
                                                          aCCTTCGTGGGCATCCACATGTGGGCCCAGGCTGTCGATAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGGTTCACCCTGACCATGGACAAGACCAACCACCACCTGca                                                                    <  7:88786/149‑1 (MQ=255)
                                                          aCCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGca                                                                    >  6:150605/1‑149 (MQ=255)
                                                               cGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGc                                                               <  7:4909/149‑1 (MQ=255)
                                                                gTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGcc                                                              >  5:238725/1‑149 (MQ=255)
                                                                          acaTGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGGTTCACCCTGACCATGGACAAGa                                                                                    <  7:46982/117‑1 (MQ=255)
                                                                          acaTGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGa                                                                                    >  8:46982/1‑117 (MQ=255)
                                                                            aTGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCgg                                                  >  6:221404/1‑149 (MQ=255)
                                                                             tgtgGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGc                                                 <  8:98827/149‑1 (MQ=255)
                                                                               tgGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCggg                                                                                                           <  4:209677/89‑1 (MQ=255)
                                                                               tgGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCggg                                                                                                           >  3:209677/1‑89 (MQ=255)
                                                                                gggCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGc                                                 >  5:125413/1‑146 (MQ=255)
                                                                                 ggCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCgaga                                             <  1:194432/149‑1 (MQ=255)
                                                                                  gCCCAGGCTGTCGAGAAAGCCAAGTCGCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGc                                                 >  1:160970/1‑145 (MQ=255)
                                                                                  gCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGtg                                                          >  7:168458/1‑135 (MQ=255)
                                                                                  gCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGtg                                                          <  8:168458/135‑1 (MQ=255)
                                                                                  gCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGAt                                            >  6:49061/1‑149 (MQ=255)
                                                                                   cccAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGacc                                                                                  <  8:232435/110‑1 (MQ=255)
                                                                                   cccAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGacc                                                                                  >  7:232435/1‑110 (MQ=255)
                                                                                       ggCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAgg                                       <  5:150605/149‑1 (MQ=255)
                                                                                        gCTGGCGAGAAAGCGAAGTCCACTGACGTTGACAAGGTCCGGGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGa                                      <  6:310089/149‑1 (MQ=255)
                                                                                         cTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGAt                                     <  6:238725/149‑1 (MQ=255)
                                                                                            tCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGAc                                  <  5:144076/149‑1 (MQ=255)
                                                                                            tCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGAc                                  >  4:247511/1‑149 (MQ=255)
                                                                                              gagaAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGACgg                                >  5:277057/1‑149 (MQ=255)
                                                                                               agaAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGACGGc                               <  5:221404/149‑1 (MQ=255)
                                                                                               agaAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGACGGc                               <  6:277057/149‑1 (MQ=255)
                                                                                                 aaaGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGCATGACGGc                               >  4:8970/1‑147 (MQ=255)
                                                                                                        aGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGACGGCCAGTTCAGc                      <  5:49061/149‑1 (MQ=255)
                                                                                                          tCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGACGGCCAGTTCAGCgt                    <  2:138952/149‑1 (MQ=255)
                                                                                                              aCTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGACGGCCAGTTCAGCgtggt                 >  6:57211/1‑149 (MQ=255)
                                                                                                                  aCGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGACGGCCAGTTCAGCGTGgtgt               <  8:292595/147‑1 (MQ=255)
                                                                                                                  aCGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGACGGCCAGTTCAGCGTGGtgtg              <  3:247511/148‑1 (MQ=255)
                                                                                                                    gTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGtga                                                         <  6:67385/103‑1 (MQ=255)
                                                                                                                        aCAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGACGGCCAGTTCAGCGTGGTGTGGGAAAc        <  8:13900/148‑1 (MQ=255)
                                                                                                                         cAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGACGGCCAGTTCAGCGTGGTGTGGGAAACCg      <  2:160970/149‑1 (MQ=255)
                                                                                                                          aaGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGACGGCCAGTTCAGCGTGGTGTGGGAAACCGa     <  3:8970/149‑1 (MQ=255)
                                                                                                                             gTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGACGGCCAGTTCAGCGTGGTGTGGGAAACCGagca  <  3:301748/149‑1 (MQ=255)
                                                                                                                                ctcgAAGCTCTGCCCGGCCAGTGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGa                                      <  4:56909/108‑1 (MQ=39)
                                                                                                                                ctcgAAGCTCTGCCCGGCCAGTGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGa                                      >  3:56909/3‑110 (MQ=39)
                                                                                                                                             ccGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCa                                         >  2:157514/1‑94 (MQ=255)
                                                                                                                                             ccGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCa                                         <  1:157514/94‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                            
CAAGGCCAAGGGCCTGCCAGGTGCCGACAAGGCCGTGACCAACGACCCCATGGAAGCCACCTACGTGGGCATCCACATGTGGGCCCAGGCTGTCGAGAAAGCCAAGTCCACTGACGTTGACAAGGTCCGCGAAGCACTGGCCGGGCAGAGCTTCGAGGCGCCGTCGGGCTTCACCCTGACCATGGACAAGACCAACCACCACCTGCACAAGCCGGTGATGATCGGCGAGATTCAGGATGACGGCCAGTTCAGCGTGGTGTGGGAAACCGAGCA  >  NC_002947/5509310‑5509582

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