Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,650,147 G→A 15.2% pseudogene (70/765 nt) ydfE → Qin prophage; pseudogene;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,650,1470GA15.2% 282.5 / 35.0 112pseudogene (70/765 nt)ydfEQin prophage; pseudogene;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base G (44/51);  new base A (8/9);  total (52/60)
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.87e-01

Only 100 of 112 total aligned reads displayed.

TCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTAT  >  NC_000913/1650049‑1650245
                                                                                                  |                                                                                                  
tCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGt                                                                                                   >  1:641614/1‑100 (MQ=255)
 ccTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTa                                                                                                  >  1:1121693/1‑100 (MQ=255)
 ccTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTa                                                                                                  <  2:641617/100‑1 (MQ=255)
 ccTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTa                                                                                                  <  2:373872/100‑1 (MQ=255)
 ccTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTa                                                                                                  >  2:1173562/1‑100 (MQ=255)
  cTTAACTCTGTCATGAGGTTAATACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCCGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGt                                                                                                   >  1:1689941/1‑98 (MQ=255)
  cTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGtcg                                                                                                 >  1:1771284/1‑98 (MQ=255)
  cTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAg                                                                                                 <  1:410890/100‑1 (MQ=255)
  cTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAg                                                                                                 >  1:519689/1‑100 (MQ=255)
  cTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAg                                                                                                 <  2:1142931/100‑1 (MQ=255)
  cTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAg                                                                                                 >  2:2692016/1‑100 (MQ=255)
   ttAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGc                                                                                                <  1:1082955/100‑1 (MQ=255)
   ttAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGc                                                                                                <  1:667753/100‑1 (MQ=255)
   ttAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGc                                                                                                <  2:879622/100‑1 (MQ=255)
    tAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGcc                                                                                               <  1:673678/100‑1 (MQ=255)
    tAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGcc                                                                                               <  1:2137207/100‑1 (MQ=255)
    tAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCCGTAGcc                                                                                               >  1:1853758/1‑100 (MQ=255)
     aaaTCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCca                                                                                              >  1:1992914/1‑100 (MQ=255)
       aTCTGTCATGAGGTTAAAACAATTGAGTAAAGTCTTTTTTTGCGCCGCTATGCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCcaca                                                                                            <  1:1964193/100‑1 (MQ=255)
       aTCTGTCATGAGGTTAAAACAACTGAGTAAAGTCTTTGTTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCcaca                                                                                            <  1:1964238/100‑1 (MQ=255)
       aTCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCcaca                                                                                            <  1:1645661/100‑1 (MQ=255)
        tCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAg                                                                                           >  2:2117897/1‑100 (MQ=255)
          tGTCATGAGGTTAAAACAAATTAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTACAGCCACAGcg                                                                                         <  2:2666360/100‑2 (MQ=255)
          tGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGcc                                                                                         <  2:2808275/100‑1 (MQ=255)
          tGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGcc                                                                                         >  2:1300752/1‑100 (MQ=255)
          tGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGcc                                                                                         <  2:956437/100‑1 (MQ=255)
          tGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTACAGCCACAGcg                                                                                         <  1:1637994/100‑2 (MQ=255)
           gTCATGGTGTTAAAACAAATGAGTAAAGTCTTTATTTGGGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAgn                                                                                          <  1:1073480/98‑2 (MQ=255)
           gTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTACAGCCACAGcgc                                                                                        <  2:2529265/100‑3 (MQ=255)
           gTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTACAGCCACAGcgc                                                                                        <  2:2895970/100‑3 (MQ=255)
           gTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTACAGCCACAGcgc                                                                                        <  2:375644/100‑3 (MQ=255)
             cATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCAtt                                                                                      >  2:1925768/1‑100 (MQ=255)
             cATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTACAGCCACAGcgcct                                                                                      >  2:2494130/1‑96 (MQ=255)
             cATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTACAGCCACAGcgcct                                                                                      <  1:1595948/100‑5 (MQ=255)
              aTGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTACAGCCACAGcgcctt                                                                                     <  1:1324915/100‑6 (MQ=255)
               tGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGa                                                                                    >  2:1332415/1‑100 (MQ=255)
               tGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTACAGCCACAGcgccttc                                                                                    <  1:2494115/100‑7 (MQ=39)
                 atgTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAg                                                                                  <  2:240878/98‑1 (MQ=255)
                  ggTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGc                                                                                 >  2:836863/1‑100 (MQ=255)
                   gTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTg                                                                                     >  1:2684945/1‑95 (MQ=255)
                     tAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTgg                                                                              <  1:2903525/100‑1 (MQ=255)
                      aaaaCAAATGAGTAAAGTCTTTATCTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGc                                                                             <  1:703002/100‑1 (MQ=255)
                         aCAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGcgac                                                                          <  1:1332410/100‑1 (MQ=255)
                          cAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGcgacg                                                                         >  1:2100205/1‑100 (MQ=255)
                          cAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGcgac                                                                          >  1:239106/1‑99 (MQ=255)
                                 gTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCg                                                                  <  1:2749571/100‑1 (MQ=255)
                                 gTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCg                                                                  <  2:2100218/100‑1 (MQ=255)
                                  tAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGt                                                                 <  2:1399522/100‑1 (MQ=255)
                                    aaGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGc                                                               <  2:2720760/100‑1 (MQ=255)
                                     aGTTTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTAGCGACGAACGCCg                                                                  <  2:239107/96‑1 (MQ=255)
                                      gTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTc                                                             >  2:75961/1‑100 (MQ=255)
                                        cTTTATGTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGa                                                           >  1:2433869/1‑100 (MQ=255)
                                         tttATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAg                                                          <  2:1087843/100‑1 (MQ=255)
                                         tttATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAg                                                          <  2:264899/100‑1 (MQ=255)
                                         tttATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAg                                                          >  1:8131/1‑100 (MQ=255)
                                         ttaaTTTGCGCCGCTATTCCTGACGAACTGGCCACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAg                                                          <  2:636059/97‑1 (MQ=255)
                                           tATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCa                                                        >  2:2393052/1‑100 (MQ=255)
                                           tATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCa                                                        >  2:1980425/1‑100 (MQ=255)
                                                gcgcCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAAtt                                                   <  2:1131569/100‑1 (MQ=255)
                                                  gccgcTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTc                                                 <  2:1271802/100‑1 (MQ=255)
                                                  gccgcTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTc                                                 <  2:2745723/100‑1 (MQ=255)
                                                    cgcTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCAc                                               <  2:2542599/100‑1 (MQ=255)
                                                        aTTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTgg                                            >  2:2731019/1‑99 (MQ=255)
                                                        aTTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGc                                           <  2:1992925/100‑1 (MQ=255)
                                                        aTTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGc                                           <  1:1341012/100‑1 (MQ=255)
                                                        aTTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGc                                           <  1:1065895/100‑1 (MQ=255)
                                                           ccTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAAt                                        >  1:1757924/1‑100 (MQ=255)
                                                              gACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAAtt                                       >  1:2290169/1‑98 (MQ=255)
                                                                cGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTg                                   <  1:2393037/100‑1 (MQ=255)
                                                                cGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCGCGAGCAAAATTTCACTGGCAATTCCTg                                   >  2:1573269/1‑100 (MQ=255)
                                                                 gAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTgg                                  <  2:1304802/100‑1 (MQ=255)
                                                                 gAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTgg                                  >  2:381362/1‑100 (MQ=255)
                                                                     tGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACa                              <  2:53985/100‑1 (MQ=255)
                                                                      ggCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACAt                             <  1:1949427/100‑1 (MQ=255)
                                                                        caacaaGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTa                           >  2:79186/1‑100 (MQ=255)
                                                                        caacaaGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTa                           >  2:2363848/1‑100 (MQ=255)
                                                                        caacaaGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTa                           >  2:2860079/1‑100 (MQ=255)
                                                                           caaGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATcc                        >  1:999741/1‑100 (MQ=255)
                                                                             aGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCgg                      >  2:167966/1‑100 (MQ=255)
                                                                              gggAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGgc                     >  2:2754974/1‑100 (MQ=255)
                                                                                gaagaaGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGgctg                   >  2:2609736/1‑100 (MQ=255)
                                                                                   gaagGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTc                <  1:2217692/100‑1 (MQ=255)
                                                                                    aagGCGCTGTGGCCGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCa               >  2:1958272/1‑100 (MQ=255)
                                                                                     agGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGCCTGCTCAg              >  1:2633915/1‑100 (MQ=255)
                                                                                     aatggcTGTGGCTACAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAg              >  1:401183/5‑100 (MQ=255)
                                                                                      atggcTGTGGCTACAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAgg             >  2:2426715/4‑100 (MQ=255)
                                                                                      atggcTGTGGCTACAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATGATCCGGCTGCTCAgg             >  2:2869944/4‑100 (MQ=255)
                                                                                       gcgcTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGa            >  1:2519412/1‑100 (MQ=255)
                                                                                        ggcTGTGGCTACAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATACGGCTGCTCAGGAc           >  1:258167/2‑100 (MQ=255)
                                                                                        ggcTGTGGCTACAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCCCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGAc           >  2:837295/2‑100 (MQ=255)
                                                                                         gcTGTGGCTACAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACt          >  1:2101836/1‑100 (MQ=255)
                                                                                         gcTGTGGCTACAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACt          >  2:2268130/1‑100 (MQ=255)
                                                                                             tgGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTgcgn      <  2:1278247/100‑2 (MQ=255)
                                                                                             tgGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTgcgc      >  2:130603/1‑100 (MQ=255)
                                                                                              ggCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCt     <  2:1757932/100‑1 (MQ=255)
                                                                                               gCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCtt    >  1:1395232/1‑100 (MQ=255)
                                                                                                cTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTta   <  1:999553/100‑1 (MQ=255)
                                                                                                cTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTta   <  2:1985858/100‑1 (MQ=255)
                                                                                                cTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATATCACTTGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCt     >  1:1131314/1‑98 (MQ=255)
                                                                                                 tGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTtat  >  1:817207/1‑100 (MQ=255)
                                                                                                  |                                                                                                  
TCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTAT  >  NC_000913/1650049‑1650245

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: