Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,101,813 A→T 5.5% H130L (CAT→CTT)  mlaD → phospholipid ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,101,8130AT 5.5% 611.8 / 17.1 180H130L (CAT→CTT) mlaDphospholipid ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (80/90);  new base T (2/8);  total (82/98)
Fisher's exact test for biased strand distribution p-value = 1.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

Only 100 of 188 total aligned reads displayed.

TCGGGTGACGTTGCAGCTGGACAAGTCGGTCGATAACCTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAGCGCAACCGGCTAATTAAGGAGTACCCATGATTTCGATCCTGCGACGTGGCCTGCTGGTGCTGCTGGCGG  >  NC_002947/1101664‑1101962
                                                                                                                                                     |                                                                                                                                                     
   ggTGACGTTGCAGCTGGACAAGTCGGTCGATAACCTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCATGGATGGTTGCACCATCCTTGa                                                                                                                                                    <  2:869024/150‑1 (MQ=255)
   ggTGACGTTGCAGCTGGACAAGTCGGTCGATAACCTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCATGGATGGTTGCACCATCCTTGa                                                                                                                                                    <  1:2855276/150‑1 (MQ=255)
   ggTGACGTTGCAGCTGGACAAGTCGGTCGATAACCTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTTGCACCATCCTTGa                                                                                                                                                    <  2:4174938/150‑1 (MQ=255)
   ggTGACGTTGCAGCTGGACAAGTCGGTCGATAACCTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGa                                                                                                                                                    <  2:1781700/150‑1 (MQ=255)
   ggTGACGTTGCAGCTGGACAAGTCGGTCGATAACCTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGa                                                                                                                                                    <  1:2032344/150‑1 (MQ=255)
   ggTGACGTTGCAGCTGGACAAGTCGGTCGATAACCTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCATTGGCGGTGAAGACGAAGTGCTCAAGGATGGTTGCACCATCCTTGa                                                                                                                                                    <  2:1186721/150‑1 (MQ=255)
         gTTGCAGCTGGACAAGTCGGTCGATAACCTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGAcacaca                                                                                                                                              >  2:3548114/1‑150 (MQ=255)
          ttGCAGCTGGACAAGTCGGTCGATAACCTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAg                                                                                                                                             >  2:3588556/1‑150 (MQ=255)
            gCAGCTGGACAAGTCGGTCGATAACCTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTc                                                                                                                                           >  2:4282607/1‑150 (MQ=255)
             cAGCTGGACAAGTCGGTCGATAACCTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCg                                                                                                                                          <  1:5218729/150‑1 (MQ=255)
              aGCTGGACAAGTCGGTCGATAACCTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCgg                                                                                                                                         <  2:2733422/150‑1 (MQ=255)
               gCTGGACAAGTCGGTCGATAACCTGCCGACCGACTCCACCGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGc                                                                                                                                        >  2:1346094/1‑150 (MQ=255)
                cTGGACAAGTCGGTCGATAACCTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCa                                                                                                                                       <  1:950946/150‑1 (MQ=255)
                 tGGACAAGTCGGTCGATAACCTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCAc                                                                                                                                      <  1:1786508/150‑1 (MQ=255)
                   gACAAGTCGGTCGATAACCTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTg                                                                                                                                    <  2:2848336/150‑1 (MQ=255)
                     cAAGTCGGTCGATAACCTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGt                                                                                                                                  <  2:1990719/150‑1 (MQ=255)
                      aaGTCGGTCGATAACCTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTg                                                                                                                                 >  1:1019325/1‑150 (MQ=255)
                         tcggtcgATAACCTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTg                                                                                                                              >  2:4525307/1‑150 (MQ=255)
                          cggtcgATAACCTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGCTGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTgg                                                                                                                             <  1:4002289/150‑1 (MQ=255)
                           ggtcgATAACCTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTgga                                                                                                                            <  1:2629186/150‑1 (MQ=255)
                           ggtcgATAACCTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTgga                                                                                                                            <  1:1109179/150‑1 (MQ=255)
                            gtcgATAACCTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTggag                                                                                                                           >  2:538756/1‑150 (MQ=255)
                             tcgATAACCTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTggagg                                                                                                                          <  2:438010/150‑1 (MQ=255)
                               gATAACCTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGAt                                                                                                                        >  2:3758206/1‑150 (MQ=255)
                                 tAACCTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCt                                                                                                                      >  1:4111387/1‑150 (MQ=255)
                                    ccTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGAt                                                                                                                   >  1:4217993/1‑150 (MQ=255)
                                    ccTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGAt                                                                                                                   >  1:2715783/1‑150 (MQ=255)
                                     ttatcgaccgacTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGAtt                                                                                                                  <  1:4513784/146‑1 (MQ=255)
                                     cTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGAtt                                                                                                                  >  2:3397618/1‑150 (MQ=255)
                                         cgaccgacTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCa                                                                                                              >  1:5014625/1‑150 (MQ=255)
                                             cgacTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCt                                                                                                                      <  1:3173010/138‑1 (MQ=255)
                                              gacTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTc                                                                                                         >  2:4721275/1‑150 (MQ=255)
                                              gacTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTc                                                                                                         <  2:1672780/150‑1 (MQ=255)
                                                cTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCt                                                                                                       <  2:3280668/150‑1 (MQ=255)
                                                 tCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTg                                                                                                      <  1:4721275/150‑1 (MQ=255)
                                                 tCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTg                                                                                                      <  1:4282607/150‑1 (MQ=255)
                                                  ccACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGc                                                                                                     >  1:1036104/1‑150 (MQ=255)
                                                  ccACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGc                                                                                                     <  2:1217081/150‑1 (MQ=255)
                                                   cACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCttgc                                                                                                    <  1:4570053/150‑4 (MQ=255)
                                                   cACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCt                                                                                                    >  2:4751957/1‑150 (MQ=255)
                                                   cACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCAACAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCt                                                                                                    <  2:2073401/150‑1 (MQ=255)
                                                    aCTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTc                                                                                                   >  1:5417198/1‑150 (MQ=255)
                                                      tGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCaa                                                                                                 >  2:1255511/1‑150 (MQ=255)
                                                        ccTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAAct                                                                                               <  1:4852286/150‑1 (MQ=255)
                                                         cTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAActc                                                                                              >  2:2932203/1‑150 (MQ=255)
                                                          tCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAActct                                                                                             >  1:1896549/1‑150 (MQ=255)
                                                           cGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTg                                                                                            <  2:3664338/150‑1 (MQ=255)
                                                           cGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTg                                                                                            >  1:215784/1‑150 (MQ=255)
                                                            gATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGt                                                                                           <  2:2663096/150‑1 (MQ=255)
                                                             aTCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGtt                                                                                          <  2:618162/150‑1 (MQ=255)
                                                              tCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTg                                                                                         >  1:2812281/1‑150 (MQ=255)
                                                                  gACCGCCGGTTGGCTGGGTGAGAAGTAGATGGGCATCTGCGTTGGCGGTGTAGACGATGTGCTCAAGGATGGTGCAACCATCCATGACACAGAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCACCTCTGTTGGcta                                                                                     <  2:4197568/150‑3 (MQ=255)
                                                                     cgccgGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGA‑‑TGATTGGCAAGTTCCTg                                                                                                      <  1:3363802/128‑1 (MQ=255)
                                                                      gccgGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGaa                                                                                 <  1:4437432/150‑1 (MQ=255)
                                                                          gTTTGCTTGGTGAGAAGAACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGa                                                                             >  2:1644694/1‑150 (MQ=255)
                                                                           tttGCTTGGTGAGAAGTACATCTGCATCATCGTTGGCGTTGAAGACGAAGTGCTCAAGGATTGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGaa                                                                            <  2:5104410/150‑1 (MQ=255)
                                                                           tttGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGaa                                                                            <  2:1685386/150‑1 (MQ=255)
                                                                           tttGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGaa                                                                            <  1:3758206/150‑1 (MQ=255)
                                                                           tttGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGaa                                                                            <  1:3851268/150‑1 (MQ=255)
                                                                            ttGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTc                                                                                                         <  1:3798333/120‑1 (MQ=255)
                                                                            ttGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTc                                                                                                         >  2:3798333/1‑120 (MQ=255)
                                                                             tGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAaag                                                                          <  1:4015415/150‑1 (MQ=255)
                                                                              gCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAaaga                                                                         >  2:273389/1‑150 (MQ=255)
                                                                               cTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAaagaa                                                                        <  1:925908/150‑1 (MQ=255)
                                                                               cTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAaagaa                                                                        <  1:2644243/150‑1 (MQ=255)
                                                                               cTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAaagaa                                                                        <  2:3486612/150‑1 (MQ=255)
                                                                               cTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAaagaa                                                                        <  2:4003503/150‑1 (MQ=255)
                                                                                ttGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAaagaag                                                                       >  2:361915/1‑150 (MQ=255)
                                                                                 tGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAgc                                                                      <  2:2548351/150‑1 (MQ=255)
                                                                                 tGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAgc                                                                      <  2:1896549/150‑1 (MQ=255)
                                                                                  ggTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAgcg                                                                     <  2:3357545/150‑1 (MQ=255)
                                                                                       gaAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAGCGCAAcc                                                                <  2:3128559/150‑1 (MQ=255)
                                                                                          gTACATCGGCATCAGCGTTGGAGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAGCGCAACCGGc                                                             <  2:1087781/150‑1 (MQ=255)
                                                                                           tACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAGCGCAACCGGCt                                                            <  2:2363258/150‑1 (MQ=255)
                                                                                                cGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAGCGCAACCGGCTAATTa                                                       >  1:1163114/1‑150 (MQ=255)
                                                                                                 ggCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAGCGCAACCGGCTAATTaa                                                      <  1:1185662/150‑1 (MQ=255)
                                                                                                  gCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAGCGCAACCGGCTAATTAAg                                                     <  2:1019325/150‑1 (MQ=255)
                                                                                                   cATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAGCGCAACCGGCTAATTAAgg                                                    >  1:4529709/1‑150 (MQ=255)
                                                                                                    aTCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCTAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAGCGCAACCGGCTAATTAAGGa                                                   <  2:4836062/150‑1 (MQ=255)
                                                                                                    aTCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAGCGCAACCGGCTAATTAAGGa                                                   <  2:4553782/150‑1 (MQ=255)
                                                                                                        gCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAGCGCAACCGGCTAATTAAGGAGTAc                                               <  2:929074/150‑1 (MQ=255)
                                                                                                         cGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAGCGCAACCGGCTAATTAAGGAGTAcc                                              >  1:1721208/1‑150 (MQ=255)
                                                                                                               cGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAGCGCAACCGGCTAATTAAGGAGTACCCATGAt                                        <  2:1527144/150‑1 (MQ=255)
                                                                                                                    aaGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAGCGCAACCGGCTAATTAAGGAGTACCCATGATTTCGa                                   >  1:2330516/1‑150 (MQ=255)
                                                                                                                         gAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAGCGCAACCGGCTAATTAAGGAGTCCCCATGATTTCGATCCTg                              >  2:4993469/1‑150 (MQ=255)
                                                                                                                         gAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAGCGCAACCGGCTAATTAAGGAGTACCCATGATTTCGAt                                  <  2:4224679/146‑1 (MQ=255)
                                                                                                                         gAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAGCGCAACCGGCTAATTAAGGAGTACCCATGATTTCGAt                                  >  1:4224679/1‑146 (MQ=255)
                                                                                                                                tCATGGATGGTTGCACCATCCTTGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAGCGCAACCGGCTAATTAAGGAGTACCCATGAttt                                      <  1:1321033/135‑1 (MQ=255)
                                                                                                                                tCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAGCGCAACCGGCTAATTAAGGAGTACCCATGATTTCGATCCTGCGACGTg                       >  1:2705987/1‑150 (MQ=255)
                                                                                                                                 cAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAGCGCAACCGGCTAATTAAGGAGTACCCATGATTTCGATCCTGCGACGTgg                      >  1:3876166/1‑150 (MQ=255)
                                                                                                                                  aaGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAGCGCAACCGGCTAATTAAGGAGTACCCATGATTTCGATCCTGCGACGTGGc                     >  1:2565611/1‑150 (MQ=255)
                                                                                                                                  aaGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAGCGCAACCGGCTAATTAAGGAGTACCCATGATTTCGATCCTGCGACGTGGc                     >  1:2797724/1‑150 (MQ=255)
                                                                                                                                  aaGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAGCGCAACCGGCTAATTAAGGAGTACCCATGATTTCGATCCTGCGACGTGGc                     >  1:1808736/1‑150 (MQ=255)
                                                                                                                                  aaGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAGCGCAACCGGCTAATTAAGGAGTACCCATGATTTCGATCCTGCGACGTGGc                     >  1:4164212/1‑150 (MQ=255)
                                                                                                                                  aaGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAGCGCAACCGGCTAATTAAGGAGTACCCATGATTTCGATCCTGCGACGTGGc                     >  1:818421/1‑150 (MQ=255)
                                                                                                                                    ggATGGTTGCACCATCCTTGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAGCGCAACCGGCTAATTAAGGAGTAcc                                              >  1:583017/1‑123 (MQ=255)
                                                                                                                                    ggATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGACCCGAAAGAAGCGCAACCGGCTAATTAAGGAGTACCCATGATTTCGATCCTGCGACGTGGCct                   >  2:1317305/1‑150 (MQ=255)
                                                                                                                                          tGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAGCGCAACCGGCTAATTAAGGAGTACCCATGATTTCGATCCTGCGACGTGGCCTGCTGGt             <  2:4436368/150‑1 (MQ=255)
                                                                                                                                              aCCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAGCGCAACCGGCTAATTAAGGAGTACCCATGATTTCGATCCTGCGACGTGGCCTGCTGGtgctg         <  2:3879945/150‑1 (MQ=255)
                                                                                                                                                  tccatGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAGCGCAACCGGCTAATTAAGGAGTACCCATGATTTCGATCCTGCGACGTGGCCTGCTGGTGCTGCTgg     <  1:3587004/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCGGGTGACGTTGCAGCTGGACAAGTCGGTCGATAACCTGCCGACCGACTCCACTGCCTCGATCCTGACCGCCGGTTTGCTTGGTGAGAAGTACATCGGCATCAGCGTTGGCGGTGAAGACGAAGTGCTCAAGGATGGTGCAACCATCCATGACACACAGTCGGCACTGGTGCTGGAGGATCTGATTGGCAAGTTCCTGCTCAACTCTGTTGGCAAGGAACCGAAAGAAGCGCAACCGGCTAATTAAGGAGTACCCATGATTTCGATCCTGCGACGTGGCCTGCTGGTGCTGCTGGCGG  >  NC_002947/1101664‑1101962

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