Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,674,003 G→A A229V (GCA→GTA)  yphB ← putative aldose 1‑epimerase YphB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,674,0030GA99.7% 1320.7 / NA 312A229V (GCA→GTA) yphBputative aldose 1‑epimerase YphB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (163/148);  minor base C (0/1);  total (163/149)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 312 total aligned reads displayed.

GATCATCCGGCGCATGGCTCATCGGTTCCAGACAGAAGAAATCAAACGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTGCAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCGGCTGCTCGATGCGGGCCTGACCATTCCATCCGGCAAAACCATTGTTCACCCACTG  >  NC_000913/2673905‑2674100
                                                                                                  |                                                                                                 
 atcatcCGGCGCATGGCTCATCGGTTCCAGACAGAAGAAATCAAACGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAg                                                                                                >  1:4855425/1‑101 (MQ=255)
  tcatcCGGCGCATGGCTCATCGGTTCCAGACAGAAGAAATCAAACGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAgg                                                                                               >  2:4855426/1‑101 (MQ=255)
   catcCGGCGCATGGCTCATCGGTTCCAGACAGAAGAAATCAAACGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGt                                                                                              <  1:4854945/101‑1 (MQ=255)
   catcCGGCGCATGGCTCATCGGTTCCAGACAGAAGAAATCAAACGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGt                                                                                              <  1:4855046/101‑1 (MQ=255)
   catcCGGCGCATGGCTCATCGGTTCCAGACAGAAGAAATCAAACGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGt                                                                                              >  1:4855429/1‑101 (MQ=255)
      ccGGCGCATGGCTCATCGGTTCCAGACAGAAGAAATCAAACGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGc                                                                                           >  1:4855435/1‑101 (MQ=255)
      ccGGCGCATGGCTCATCGGTTCCAGACAGAAGAAATCAAACGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGc                                                                                           >  2:4855434/1‑101 (MQ=255)
        ggCGCATGGCTCATCGGTTCCAGACAGAAGAAATCAAACGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGcgt                                                                                         >  1:4855439/1‑101 (MQ=255)
          cgcATGGCTCATCGGTTCCAGACAGAAGAAATCAAACGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGcgtcg                                                                                       <  1:4854975/101‑1 (MQ=255)
           gcATGGCTCATCGGTTCCAGACAGAAGAAATCAAACGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGcgtcgt                                                                                      <  2:4855001/101‑1 (MQ=255)
            cATGGCTCATCGGTTCCAGACAGAAGAAATCAAACGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGtt                                                                                     >  1:4855443/1‑101 (MQ=255)
             aTGGCTCATCGGTTCCAGACAGAAGAAATCAAACGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGttt                                                                                    >  1:4855446/1‑101 (MQ=255)
               ggCTCATCGGTTCCAGACAGAAGAAATCAAACGCATATCCTTTATCAAACGCAGGGTCGGAAACAAGGATGAAATAACACGGTACAGGTGGCGTCGTTTcc                                                                                  <  1:4854929/101‑1 (MQ=255)
               ggCTCATCGGTTCCAGACAGAAGAAATCAAACGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTcc                                                                                  <  2:4855023/101‑1 (MQ=255)
               ggCTCATCGGTTCCAGACAGAAGAAATCAAACGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTcc                                                                                  <  1:4854941/101‑1 (MQ=255)
               ggCTCATCGGTTCCAGACAGAAGAAATCAAACGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTcc                                                                                  >  2:4855449/1‑101 (MQ=255)
                 cTCATCGGTTCCAGACAGAAGAAATCAAACGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCat                                                                                <  2:4855180/101‑1 (MQ=255)
                   cATCGGTTCCAGACAGAAGAAATCAAACGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCatga                                                                              >  1:4855453/1‑101 (MQ=255)
                   cATCGGTTCCAGACAGAAGAAATCAAACGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCatga                                                                              >  2:4855454/1‑101 (MQ=255)
                   cATCGGTTCCAGACAGAAGAAATCAAACGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCatga                                                                              <  1:4855065/101‑1 (MQ=255)
                        gTTCCAGACAGAAGAAATCAAACGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCatgatgatg                                                                         <  1:4855159/101‑1 (MQ=255)
                        gTTCCAGACAGAAGAAATCAAACGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCatgatgatg                                                                         >  2:4855459/1‑101 (MQ=255)
                         ttCCAGACAGAAGAAATCAAACGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGtt                                                                                     >  1:4855461/1‑88 (MQ=255)
                         ttCCAGACAGAAGAAATCAAACGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATgg                                                                        <  2:4854953/101‑1 (MQ=255)
                           ccAGACAGAAGAAATCAAACGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCa                                                                      <  1:4854816/101‑1 (MQ=255)
                           ccAGACAGAAGAAATCAAACGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCa                                                                      <  1:4854890/101‑1 (MQ=255)
                            cagacagaAGAAATCAAACGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCa                                                                      <  1:4855234/100‑1 (MQ=255)
                             agacagaAGAAATCAAACGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATa                                                                    <  1:4855143/101‑1 (MQ=255)
                             agacagaAGAAATCAAACGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATa                                                                    <  2:4854775/101‑1 (MQ=255)
                             agacagaAGAAATCAAACGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATa                                                                    >  2:4855465/1‑101 (MQ=255)
                               acagaAGAAATCAAACGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAAc                                                                  <  2:4854925/101‑1 (MQ=255)
                                cagaAGAAATCAAACGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAAcc                                                                 <  2:4854989/101‑1 (MQ=255)
                                cagaAGAAATCAAACGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAAcc                                                                 >  1:4855468/1‑101 (MQ=255)
                                cagaAGAAATCAAACGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAAcc                                                                 >  2:4855472/1‑101 (MQ=255)
                                 agaagaAATCAAACGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCt                                                                <  1:4855015/101‑1 (MQ=255)
                                 agaagaAATCAAACGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCt                                                                >  1:4855476/1‑101 (MQ=255)
                                 agaagaAATCAAACGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCt                                                                >  2:4855475/1‑101 (MQ=255)
                                   aagaaATCAAACGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTc                                                              >  1:4855479/1‑101 (MQ=255)
                                    agaaATCAAACGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTcc                                                             >  2:4855481/1‑101 (MQ=255)
                                      aaaTCAAACGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTg                                                           <  2:4855147/101‑1 (MQ=255)
                                      aaaTCAAACGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTg                                                           <  1:4855099/101‑1 (MQ=255)
                                       aaTCAAACGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGc                                                          <  2:4855062/101‑1 (MQ=255)
                                       aaTCAAACGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGc                                                          >  1:4855485/1‑101 (MQ=255)
                                         tCAAACGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCgg                                                        >  2:4855488/1‑101 (MQ=255)
                                         tCAAACGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCgg                                                        >  1:4855494/1‑101 (MQ=255)
                                         tCAAACGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCgg                                                        >  1:4855491/1‑101 (MQ=255)
                                           aaaCGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCGgct                                                      <  2:4855132/101‑1 (MQ=255)
                                           aaaCGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCGgct                                                      >  1:4855498/1‑101 (MQ=255)
                                           aaaCGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCGgct                                                      <  1:4855084/101‑1 (MQ=255)
                                           aaaCGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCGgct                                                      >  2:4855499/1‑101 (MQ=255)
                                           aaaCGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCGgct                                                      <  2:4854966/101‑1 (MQ=255)
                                             aCGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCt                                                            >  2:4855505/1‑93 (MQ=255)
                                             aCGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCGgctgc                                                    <  1:4855080/101‑1 (MQ=255)
                                              cGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCGgctgct                                                   >  1:4855506/1‑101 (MQ=255)
                                                  tatCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCGGCTGCTCGAt                                               >  2:4855510/1‑101 (MQ=255)
                                                   atCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGcgtcgt                                                                                      >  2:4855514/1‑61 (MQ=255)
                                                    tCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCGGCTGCTCGATGc                                             >  1:4855516/1‑101 (MQ=255)
                                                      cTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCGGCTGCTCGATGc                                             <  2:4855091/99‑1 (MQ=255)
                                                         tATCAAACGCAGGGTCGGAAACAAAGATGATATAACACGGTACAGGTGGCGTCGTTTCCATAATGATGGCATAACCTTCCTGCGGCTGCTCGATGCGGGcc                                        >  1:4855520/1‑101 (MQ=255)
                                                          aTCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCGGCTGCTCGATGCGGGCCt                                       >  2:4855522/1‑101 (MQ=255)
                                                          aTCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCGGCTGCTCGATGCGGGCCt                                       <  2:4855215/101‑1 (MQ=255)
                                                           tCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCGGCTGCTCg                                                 >  2:4855525/1‑90 (MQ=255)
                                                           tCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCGGCTGCTCGATGCGGGCCTg                                      <  2:4855195/101‑1 (MQ=255)
                                                             aaaCGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCGGCTGCTCGATGCGGGCCTGAc                                    >  2:4855529/1‑101 (MQ=255)
                                                               aCGCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCGGCTGCTCGATGCGGGCCTGACCa                                  >  2:4855532/1‑101 (MQ=255)
                                                                 gCAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCGGCTGCTCGATGCGGGCCTGACCAtt                                >  1:4855534/1‑101 (MQ=255)
                                                                  cAGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCGGCTGCTCGATGCGGGCCTGACCATTc                               >  2:4855537/1‑101 (MQ=255)
                                                                   aGGGTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCGGCTGCTCGATGCGGGCCTGACCATTcc                              <  1:4855216/101‑1 (MQ=255)
                                                                     ggTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCGGCTGCTCGATGCGGGCCTGa                                     <  1:4854679/92‑1 (MQ=255)
                                                                     ggTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCGGCTGCTCGATGCGGGCCTGACCATTCCAt                            >  1:4855543/1‑101 (MQ=255)
                                                                      gTCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCGGCTGCTCGATGCGGGCCTGACCATTCCATc                           >  2:4855544/1‑101 (MQ=255)
                                                                       tCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCGGCTGCTCGATGCGGGCCTGACCATTCCATcc                          <  1:4854988/101‑1 (MQ=255)
                                                                       tCGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCGGCTGCTCGATGCGGGCCTGACCATTCCATcc                          <  1:4855185/101‑1 (MQ=255)
                                                                        cGGAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCatg                                                                               <  1:4854967/47‑1 (MQ=255)
                                                                         ggAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCGGCTGCTCGATGCGGGCCTGACCATTCCATCCgg                        >  2:4855551/1‑101 (MQ=255)
                                                                         ggAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCGGCTGCTCGATGCGGGCCTGACCATTCCATCCgg                        >  1:4855549/1‑101 (MQ=255)
                                                                          gAAACAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCGGCTGCTCGATGCGGGCCTGACCATTCCATCCGGc                       <  1:4855034/101‑1 (MQ=255)
                                                                           aaaCAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCGGCTGCTCGATGCGGGCCTGACCATTCCATCCGGCa                      >  1:4855555/1‑101 (MQ=255)
                                                                            aaCAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCGGCTGCTCGATGCGGGCCTGACCATTCCATCCGGCaa                     >  2:4855557/1‑101 (MQ=255)
                                                                              cAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCGGCTGCTCGATGCGGGCCTGACCATTCCATCCGGCaaaa                   <  2:4855116/101‑1 (MQ=255)
                                                                              cAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCGGCTGCTCGATGCGGGCCTGACCATTCCATCCGGCaaaa                   <  1:4855171/101‑1 (MQ=255)
                                                                              cAAAGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCGGCTGCTCGATGCGGGCCTGACCATTCCATCCGGCaaaa                   >  1:4855560/1‑101 (MQ=255)
                                                                               aaaGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCGGCTGCTCGATGCGGGCCTGACCATTCCATCCg                         >  2:4855565/1‑94 (MQ=255)
                                                                               aaaGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCGGCTGCTCGATGCGGGCCTGACCATTCCATCCGGCAAAAc                  <  1:4855241/101‑1 (MQ=255)
                                                                               aaaGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCGGCTGCTCGATGCGGGCCTGACCATTCCATCCGGCAAAAc                  >  2:4855562/1‑101 (MQ=255)
                                                                                aaGATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCGGCTGCTCGATGCGGGCCTGACCATTCCATCCGGCAAAAcc                 <  2:4855038/101‑1 (MQ=255)
                                                                                  gATGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCGGCTGCTCGATGCGGGCCTGACCATTCCATCCGGCAAAACCAt               >  2:4855569/1‑101 (MQ=255)
                                                                                   aTGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCGGCTGCTCGATGCGGGCCTGACCATTCCATCCGGCAAAACCAtt              <  2:4855104/101‑1 (MQ=255)
                                                                                    tGAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCGGCTGCTCGATGCGGGCCTGACCATTCCATCCGGCAAAACCATTg             <  1:4855121/101‑1 (MQ=255)
                                                                                     gAAATAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCGGCTGCTCGATGCGGGCCTGACCATTCCATCCGGCAAAACCATTGt            >  1:4855573/1‑101 (MQ=255)
                                                                                      aaaTAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCGGCTGCTCGATGCGGGCCTGACCATTCCATCCGGCAAAACCATTGtt           >  2:4855574/1‑101 (MQ=255)
                                                                                        aTAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGc                                                                       >  2:4855578/1‑39 (MQ=255)
                                                                                         tAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCGGCTGCTCGATGCGGGCCTGACCATTCCATCCGGCAAAACCATTGTTCAc        <  2:4855288/101‑1 (MQ=255)
                                                                                         tAACACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCGGCTGCTCGATGCGGGCCTGACCATTCCATCCGGCAAAACCATTGTTCAc        <  2:4855164/101‑1 (MQ=255)
                                                                                          aaCACGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCGGCTGCTCGATGCGGGCCTGACCATTCCATCCGGCAAAACCATTGTTCAcc       >  1:4855581/1‑101 (MQ=255)
                                                                                            cacGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCGGCTGCTCGATGCGGGCCTGACCATTCCATCCGGCAAAACCATTGTTCACCCa     <  2:4855081/101‑1 (MQ=255)
                                                                                            cacGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCGGCTGCTCGATGCGGGCCTGACCATTCCATCCGGCAAAACCATTGTTCACCCa     <  1:4855250/101‑1 (MQ=255)
                                                                                              cGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCGGCTGCTCGATGCGGGCCTGACCATTCCATCCGGCAAAACCATTGTTCACCCACt   <  2:4855225/101‑1 (MQ=255)
                                                                                              cGGTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCGGCTGCTCGATGCGGGCCTGACCATTCCATCCGGCAAAACCATTGTTCACCCACt   >  1:4855585/1‑101 (MQ=255)
                                                                                               ggTACAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCGGCTGCTCGATGCGGGCCTGACCATTCCATCCGGCAAAACCATTGTTCACCCAc    <  1:4855319/99‑1 (MQ=255)
                                                                                                  |                                                                                                 
GATCATCCGGCGCATGGCTCATCGGTTCCAGACAGAAGAAATCAAACGCATATCCTTTATCAAACGCAGGGTCGGAAACAAAGATGAAATAACACGGTGCAGGTGGCGTCGTTTCCATGATGATGGCATAACCTTCCTGCGGCTGCTCGATGCGGGCCTGACCATTCCATCCGGCAAAACCATTGTTCACCCACTG  >  NC_000913/2673905‑2674100

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