Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,650,147 G→A 16.9% 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,1470GA16.9% 290.5 / 48.0 130pseudogene (70/765 nt)ydfEQin prophage; pseudogene;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base G (57/51);  new base A (11/11);  total (68/62)
Fisher's exact test for biased strand distribution p-value = 8.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.44e-01

Only 100 of 135 total aligned reads displayed.

TCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATA  >  NC_000913/1650049‑1650246
                                                                                                  |                                                                                                   
tCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGt                                                                                                    <  1:294581/100‑1 (MQ=255)
  cTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTa                                                                                                   >  1:1097305/1‑99 (MQ=255)
  cTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAg                                                                                                  >  2:3939466/1‑100 (MQ=255)
  cTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAg                                                                                                  <  2:3044379/100‑1 (MQ=255)
    tAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGcc                                                                                                <  2:665449/100‑1 (MQ=255)
    tAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGcc                                                                                                <  1:2877298/100‑1 (MQ=255)
    tAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGcc                                                                                                <  2:1174129/100‑1 (MQ=255)
    tAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGCCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGcc                                                                                                <  1:818642/100‑1 (MQ=255)
       aTCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTgagatcgga                                                                                             >  1:2448291/1‑92 (MQ=255)
        tCTGTCATGAGCTTAAAACGAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAg                                                                                            <  1:1898495/100‑1 (MQ=255)
          tGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGcc                                                                                          <  2:2430647/100‑1 (MQ=255)
          tGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGcc                                                                                          <  1:4204544/100‑1 (MQ=255)
          tGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTACAGCCACAGcg                                                                                          <  2:2872460/100‑2 (MQ=255)
          tGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTACAGCCACAGcg                                                                                          <  2:1095438/100‑2 (MQ=255)
          tGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTACAGCCACAGcg                                                                                          <  2:1327320/100‑2 (MQ=255)
          tGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTAACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCGACAGcc                                                                                          >  1:2093593/1‑100 (MQ=255)
           gTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCa                                                                                         >  1:2498778/1‑100 (MQ=255)
           gTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTACAGCCACAGcgc                                                                                         <  1:918360/100‑3 (MQ=255)
            tCATTAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTACAGCCACAGc                                                                                           <  1:2662092/97‑1 (MQ=255)
            tCATGAGGTTAAAACAAATGAGTAAAGTCTTGATTTGCGCCGCTATTGCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCAACAgaaa                                                                                         >  1:3458626/1‑96 (MQ=255)
            tCATGAGGTGAAAACAAATGAGTAAAGTCTTTATGTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCAt                                                                                        <  1:215437/100‑1 (MQ=255)
              aTGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGACCTGGCAACAAGGGAAGAAGGCGCTGTGGCTACAGCCACAGcgcct                                                                                       <  1:4067200/99‑5 (MQ=255)
              aTGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCttt                                                                                       >  2:1435415/1‑99 (MQ=255)
               tGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGa                                                                                     >  1:2839391/1‑100 (MQ=255)
                    ttAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTg                                                                                <  2:1685284/100‑1 (MQ=255)
                     tAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTgg                                                                               >  2:3716941/1‑100 (MQ=255)
                     tAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTgg                                                                               >  2:2361668/1‑100 (MQ=255)
                      aaaaCAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGc                                                                              >  2:1550069/1‑100 (MQ=255)
                         aCAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGcgac                                                                           <  1:26661/100‑1 (MQ=255)
                          cAAATTAGTATAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTTGCTGTATCCACAGCCATTGAAGCTGGcgacg                                                                          <  1:3907771/100‑1 (MQ=255)
                           aaaTGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGcgacg                                                                          <  1:2507811/99‑1 (MQ=255)
                            aaTGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGaa                                                                        >  2:2114217/1‑100 (MQ=255)
                             aTGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAAc                                                                       >  2:3615961/1‑100 (MQ=255)
                              tGAGGAAAGTCTTTATTTGCGCCGCTATTCCTGGCGAATCGGCGACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACg                                                                      >  1:487895/1‑100 (MQ=255)
                               gAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGc                                                                     >  1:2776247/1‑100 (MQ=255)
                                aGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGc                                                                     >  1:1834657/1‑99 (MQ=255)
                                aGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAACTGGCGACGAACGcc                                                                    >  2:4341438/1‑100 (MQ=255)
                                 gTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCg                                                                   <  2:2946436/100‑1 (MQ=255)
                                   aaaGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTg                                                                 >  1:1449937/1‑100 (MQ=255)
                                    aaGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGc                                                                <  1:2209485/100‑1 (MQ=255)
                                    aaGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGc                                                                >  2:1078946/1‑100 (MQ=255)
                                    aaGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGc                                                                >  1:3433515/1‑100 (MQ=255)
                                    aaGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGc                                                                >  2:4176296/1‑100 (MQ=255)
                                    aaGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGc                                                                <  2:551476/100‑1 (MQ=255)
                                    aaGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGc                                                                >  2:3421093/1‑100 (MQ=255)
                                     aGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCt                                                               >  2:1709354/1‑100 (MQ=255)
                                         tttATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAg                                                           <  2:2003726/100‑1 (MQ=255)
                                         tttATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAg                                                           <  1:3053307/100‑1 (MQ=255)
                                         tttATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAg                                                           <  2:2112396/100‑1 (MQ=255)
                                           tATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCAAGTGAAGCTGGCGACGAACGCCGTGCTCGa                                                            <  2:648761/97‑1 (MQ=255)
                                            aTTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCaa                                                        <  2:2085386/100‑1 (MQ=255)
                                              ttGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCaaaa                                                      <  1:3696829/100‑1 (MQ=255)
                                              tgacgccgcGATTCCTGTCGTCCTGTCAACAGTGGAAGAAGGCGCTGTGGCGGTGGCCACAGCCATTGAACCTTGCCACGAACGCCGTTCTCGAg                                                           <  1:1374993/92‑1 (MQ=255)
                                               tGCGCCGCTATTCCTGACGCACTTGCAACAAGGGAGGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTCGCGACGAACGCCGTGCTCGAg                                                           <  1:3538859/94‑1 (MQ=255)
                                                    cgcTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCAc                                                <  2:401003/100‑1 (MQ=255)
                                                    cgcTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCAc                                                <  2:3120869/100‑1 (MQ=255)
                                                    cgcTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCAc                                                >  1:912674/1‑100 (MQ=255)
                                                     gcTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACt                                               <  2:1877395/100‑1 (MQ=255)
                                                        aTTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGc                                            <  1:3615941/100‑1 (MQ=255)
                                                        aTTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGc                                            <  1:3649917/100‑1 (MQ=255)
                                                        aTTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGc                                            <  2:216411/100‑1 (MQ=255)
                                                         ttCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCa                                           <  2:1681242/100‑1 (MQ=255)
                                                            cTGACGAACTGGCACCAAGGGAAGAAGGCACTGTGGCTGTAGCCCCAGCCAATGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAAATTCACTGGCa                                           >  1:471677/1‑97 (MQ=255)
                                                             tGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTc                                       <  2:1418738/100‑1 (MQ=255)
                                                             tGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTc                                       >  2:236668/1‑100 (MQ=255)
                                                              gACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTcc                                      >  1:4052704/1‑100 (MQ=255)
                                                                cGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTg                                    >  2:2770811/1‑100 (MQ=255)
                                                                cGAACTGGCAACAAGGGAAGAAAGCGCTATGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGGGCGCGAGCAAAATTTCACTGGCAATTCCTg                                    >  1:3386329/1‑100 (MQ=255)
                                                                    cTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAAc                                >  2:3352689/1‑100 (MQ=255)
                                                                    cTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAAc                                >  1:4273615/1‑100 (MQ=255)
                                                                     tGGCAACAAGGGAAGGAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTTGCGACGAACGCCGTGCTTGAGCAAAATTTCACTGGCAATTCCTGGAACa                               <  2:4326996/100‑1 (MQ=255)
                                                                     tGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACa                               <  1:4176266/100‑1 (MQ=255)
                                                                      ggCAACAAGGGAAGAAGGCGCTGTGGCTGTAGTCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAAGTCCTGGAACAt                              >  1:589078/1‑100 (MQ=255)
                                                                      ggCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACAt                              <  2:3078044/100‑1 (MQ=255)
                                                                      ggCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACAt                              >  2:1550557/1‑100 (MQ=255)
                                                                      ggCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACAt                              >  2:3264377/1‑100 (MQ=255)
                                                                         aacaaGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTAt                           >  1:1719244/1‑100 (MQ=255)
                                                                           caaGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATcc                         <  1:274294/100‑1 (MQ=255)
                                                                              gggAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGgc                      <  1:3312000/100‑1 (MQ=255)
                                                                              gggAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGgc                      <  2:2770031/100‑1 (MQ=255)
                                                                              gggAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGgc                      >  2:2322040/1‑100 (MQ=255)
                                                                               ggAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTGTCCGgct                     >  1:3297240/1‑100 (MQ=255)
                                                                                gcagaagGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGgctg                    >  2:2976267/3‑100 (MQ=255)
                                                                                 aagaagGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGgctgc                   >  1:1908532/1‑100 (MQ=255)
                                                                                  ttcaatggcTGTGGCTACAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGgct                     >  1:3602106/8‑97 (MQ=38)
                                                                                    aagGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCa                <  1:245912/100‑1 (MQ=255)
                                                                                     agGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAg               <  1:3264361/100‑1 (MQ=255)
                                                                                     aatggcTGTGGCTACAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACAATATCCGGCTGCTCAg               >  2:4336440/5‑100 (MQ=255)
                                                                                      atggcTGTGGCTACAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAACCTTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAg               >  1:4085183/4‑99 (MQ=255)
                                                                                        ggcTGTGGCTACAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGAc            >  2:2229262/2‑100 (MQ=255)
                                                                                        ggcTGTGGCTACAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGAc            >  1:100425/2‑100 (MQ=255)
                                                                                        cgcTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGAc            >  1:1032746/1‑100 (MQ=255)
                                                                                          cTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTg          >  1:1732322/1‑100 (MQ=255)
                                                                                          cTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTg          >  2:622302/1‑100 (MQ=255)
                                                                                           tgtgGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTgc         >  1:3603528/1‑100 (MQ=255)
                                                                                            gtgGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCCGGACTgcg        <  1:3885249/100‑1 (MQ=255)
                                                                                            gtgGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTgcg        <  2:585872/100‑1 (MQ=255)
                                                                                               gCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTgcgc       <  1:1456090/98‑1 (MQ=255)
                                                                                                cTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTta    >  2:3794755/1‑100 (MQ=255)
                                                                                                  gTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTtata  <  2:1869103/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
TCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATA  >  NC_000913/1650049‑1650246

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