Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,650,147 G→A 9.7% 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,1470GA 9.7% 268.6 / 13.6 93pseudogene (70/765 nt)ydfEQin prophage; pseudogene;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base G (44/40);  new base A (4/5);  total (48/45)
Fisher's exact test for biased strand distribution p-value = 7.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.30e-01

AAATATCCGGTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGAT  >  NC_000913/1650000‑1650296
                                                                                                                                                         |                                                                                                                                                     
aaataTCCGGTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCATTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTGCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTa                                                                                                                                                     >  2:149796/1‑150 (MQ=255)
 aataTCCGGTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAg                                                                                                                                                    >  1:1499043/1‑150 (MQ=255)
  atatCCGGTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGc                                                                                                                                                   <  1:1453764/150‑1 (MQ=255)
  atatCCGGTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGc                                                                                                                                                   <  1:122024/150‑1 (MQ=255)
   tatCCGGTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGcc                                                                                                                                                  <  1:418888/150‑1 (MQ=255)
   tatCCGGTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGcc                                                                                                                                                  <  2:1511876/150‑1 (MQ=255)
    atCCGGTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCca                                                                                                                                                 >  2:899878/1‑150 (MQ=255)
      ccGGTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTTAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCcaca                                                                                                                                               <  2:1062073/150‑1 (MQ=255)
      ccGGTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCcaca                                                                                                                                               >  2:1558448/1‑150 (MQ=255)
       cGGTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAg                                                                                                                                              >  2:1065596/1‑150 (MQ=255)
        ggTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGc                                                                                                                                             <  2:304226/150‑1 (MQ=255)
        ggTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGCTTAAAACAAATGAGTAAAGGCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGc                                                                                                                                             <  2:829982/150‑1 (MQ=255)
         gTTTCATTCTTATATGCCAGCAATGTCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTAAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTACAGCCACAGcg                                                                                                                                            <  1:672116/150‑2 (MQ=255)
         gTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTACAGCCACAGcg                                                                                                                                            <  1:395252/150‑2 (MQ=255)
         gTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTACAGCCACAGcg                                                                                                                                            <  1:435391/150‑2 (MQ=255)
         gTTTCATTCTTATAAGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGCCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTACAGCCACAGcg                                                                                                                                            <  1:1100593/150‑2 (MQ=255)
          tttCATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGAAGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTACAGCCACAGcgc                                                                                                                                           <  2:413770/150‑3 (MQ=255)
            tCATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCAtt                                                                                                                                         >  2:429718/1‑150 (MQ=255)
             cATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTg                                                                                                                                        <  2:1259680/150‑1 (MQ=255)
                 cTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGc                                                                                                                                    >  1:1126778/1‑150 (MQ=255)
                    atatGCAAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTCCGCCGCTATTCCTGACGAACTGGCTACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTgg                                                                                                                                 <  1:933913/150‑1 (MQ=255)
                     tatGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGc                                                                                                                                >  2:1374145/1‑150 (MQ=255)
                     tatGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGc                                                                                                                                >  1:1005932/1‑150 (MQ=255)
                      atGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGcg                                                                                                                               >  2:604380/1‑150 (MQ=255)
                       tGCCAGCAATGGCAGGTATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGAGCCGCTATTCATGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGACACAGCCATTGAAGCTGGcga                                                                                                                              >  2:1451669/1‑150 (MQ=255)
                                   cAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCTCCGCTATTACTGACGAACTGGCAACAATGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGTCGAACGCCGttc                                                                                                                  >  2:889269/1‑148 (MQ=255)
                                    aGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCt                                                                                                                 >  2:853309/1‑150 (MQ=255)
                                      ggACTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCg                                                                                                               <  2:635501/150‑1 (MQ=255)
                                        aTTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACATCCATTGAAGCTGGCGACGAACGCCGTGCTCGAg                                                                                                             >  2:569638/1‑150 (MQ=255)
                                           tGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCaa                                                                                                          >  2:865655/1‑150 (MQ=255)
                                            gTTCATCCTTAAATCTGACATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTTTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCaaa                                                                                                         >  1:1020531/1‑150 (MQ=255)
                                                 tCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGTCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTc                                                                                                    <  1:709531/150‑1 (MQ=255)
                                                 tCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTc                                                                                                    <  1:483404/150‑1 (MQ=255)
                                                  ccTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAATGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCa                                                                                                   <  2:720913/150‑1 (MQ=255)
                                                     tAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGACGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTg                                                                                                <  1:1525527/150‑1 (MQ=255)
                                                     tAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTg                                                                                                >  2:1421010/1‑150 (MQ=255)
                                                     tAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTg                                                                                                <  1:1239288/150‑1 (MQ=255)
                                                      aaaTCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTgg                                                                                               <  2:790205/150‑1 (MQ=255)
                                                       aaTCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGc                                                                                              >  1:1516830/1‑150 (MQ=255)
                                                          cTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAAt                                                                                           >  1:523131/1‑150 (MQ=255)
                                                              cATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGCGCAACATTTCACTGGCAATTccg                                                                                       >  2:937548/1‑149 (MQ=255)
                                                              cATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCt                                                                                       >  1:150225/1‑150 (MQ=255)
                                                              cATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCt                                                                                       >  1:1026561/1‑150 (MQ=255)
                                                               aTGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTg                                                                                      <  2:1197073/150‑1 (MQ=255)
                                                                tGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTgg                                                                                     <  2:749897/150‑1 (MQ=255)
                                                                   ggTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGAAATTCCTGGaaa                                                                                  >  2:813724/1‑149 (MQ=255)
                                                                     ttAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGCGCTGTGGCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTg                                                                                      <  2:1322594/150‑1 (MQ=255)
                                                                      tAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACAtt                                                                               >  2:1404442/1‑150 (MQ=255)
                                                                       aaaaCAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAAATTCACTGGCAATTCCTGGACCATTa                                                                              >  2:1264076/1‑150 (MQ=255)
                                                                        aaaCAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGCACATTAt                                                                             >  1:293209/1‑150 (MQ=255)
                                                                          aCAAATGAGTAAAGGCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATcc                                                                           <  1:228585/150‑1 (MQ=255)
                                                                            aaaTGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGAAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCgg                                                                         >  2:607744/1‑150 (MQ=255)
                                                                             aaTGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGgc                                                                        <  1:1374145/150‑1 (MQ=255)
                                                                             aaTGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGgc                                                                        <  1:595386/150‑1 (MQ=255)
                                                                                  gTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGCACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGCAACATTATCCGGCTGCTc                                                                   >  2:1520862/1‑150 (MQ=255)
                                                                                  gTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTTGAACATTATCCGGCTGCTc                                                                   <  1:899878/150‑1 (MQ=255)
                                                                                  gTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTc                                                                   <  1:976019/150‑1 (MQ=255)
                                                                                  gTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTc                                                                   >  2:1527553/1‑150 (MQ=255)
                                                                                  gTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTc                                                                   >  2:1345531/1‑150 (MQ=255)
                                                                                     aaGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGTCAATTCCTGGAACATTATCCGGCTGCTCAgg                                                                <  2:72779/150‑1 (MQ=255)
                                                                                       gTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGAc                                                              >  1:1149298/1‑150 (MQ=255)
                                                                                       gTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGAc                                                              >  1:1092841/1‑150 (MQ=255)
                                                                                         cTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACCGCCATTGACGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATACGGCTGCTCAGGACTg                                                            >  2:833387/1‑150 (MQ=255)
                                                                                         cTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTg                                                            <  1:70271/150‑1 (MQ=255)
                                                                                         cTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTg                                                            >  2:360973/1‑150 (MQ=255)
                                                                                         cTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTg                                                            <  1:176233/150‑1 (MQ=255)
                                                                                            tATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGTC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCAATTGAAGCTGGCCACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTgcgc                                                         >  1:926155/1‑150 (MQ=255)
                                                                                               ttGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTta                                                      <  1:63059/150‑1 (MQ=255)
                                                                                                 gcgcCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTtata                                                    <  2:1458195/150‑1 (MQ=255)
                                                                                                   gccgcTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATaaa                                                  <  1:1544671/150‑1 (MQ=255)
                                                                                                       cTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGACGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCAGGCTGCTCAGGAGTGCGCTTATAAATTTa                                              >  2:1277172/1‑150 (MQ=255)
                                                                                                         aTTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTAtt                                            <  1:604380/150‑1 (MQ=255)
                                                                                                          ttCCTGACGAACTGGCAACAACGGAAGAAGGC‑‑‑‑‑‑TCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTg                                           <  2:98147/150‑1 (MQ=255)
                                                                                                            ccTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTc                                         <  1:853309/150‑1 (MQ=255)
                                                                                                               gACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGc                                      <  1:390029/150‑1 (MQ=255)
                                                                                                               gACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGc                                      <  1:1591161/150‑1 (MQ=255)
                                                                                                                    aCTGGCAACAAGTGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGa                                 >  2:731606/1‑150 (MQ=255)
                                                                                                                     cTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCATGGAACATAATCCGGCTGCTCAGGACTGCGCTAATAAATTTATTGTCTGCGAGGAt                                >  2:1710091/1‑150 (MQ=255)
                                                                                                                      tGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑TCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATa                               <  1:248316/150‑1 (MQ=255)
                                                                                                                       ggCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATaa                              <  2:429523/150‑1 (MQ=255)
                                                                                                                       ggCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATaa                              >  1:556808/1‑150 (MQ=255)
                                                                                                                             aaGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTgg                        >  1:931402/1‑150 (MQ=255)
                                                                                                                             aaGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTgg                        <  1:195440/150‑1 (MQ=255)
                                                                                                                               gggAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCCCTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCa                      >  2:987113/1‑150 (MQ=255)
                                                                                                                               gggAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCa                      <  2:1005932/150‑1 (MQ=255)
                                                                                                                                   agaagGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAATTTTCACTGGCAATTCCTGGAACAATATCAGGCTGCTCAGGACTGAGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATAcc                  >  2:869467/1‑150 (MQ=255)
                                                                                                                                      agGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCg               >  2:1107833/1‑150 (MQ=255)
                                                                                                                                            aatggcTGTGGCTACAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCATTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCg               >  1:1539257/5‑150 (MQ=255)
                                                                                                                                            aatggcTGTGGCTACAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCg               >  1:1095358/5‑150 (MQ=255)
                                                                                                                                            aatggcTGTGGCTACAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTAACCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCg               >  1:974722/5‑150 (MQ=255)
                                                                                                                                             atggcTGTGGCTACAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCgc              >  1:1519463/4‑150 (MQ=255)
                                                                                                                                                    tgGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCgccctgccc       <  2:1020531/150‑1 (MQ=255)
                                                                                                                                                     ggCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAATAGTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTTCGCGTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCgccctgcccg      >  1:1388614/1‑149 (MQ=255)
                                                                                                                                                         gTAGACACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGAt  <  1:865655/150‑1 (MQ=255)
                                                                                                                                                         |                                                                                                                                                     
AAATATCCGGTTTCATTCTTATATGCCAGCAATGGCAGGGATTTGTTCATCCTTAAATCTGTCATGAGGTTAAAACAAATGAGTAAAGTCTTTATTTGCGCCGCTATTCCTGACGAACTGGCAACAAGGGAAGAAGGC‑‑‑‑‑‑GCTGTGGCTGTAGCCACAGCCATTGAAGCTGGCGACGAACGCCGTGCTCGAGCAAAATTTCACTGGCAATTCCTGGAACATTATCCGGCTGCTCAGGACTGCGCTTATAAATTTATTGTCTGCGAGGATAAACCTGGCATACCCCGCCCTGCCCTCGAT  >  NC_000913/1650000‑1650296

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