Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 658,153 A→G 51.5% pseudogene (123/195 nt) ybeM → pseudogene, C‑N hydrolase superfamily homology;putative enzyme; Not classified; putative amidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913658,1530AG51.5% 18.5 / 205.7 151pseudogene (123/195 nt)ybeMpseudogene, C‑N hydrolase superfamily homology;putative enzyme; Not classified; putative amidase
Reads supporting (aligned to +/- strand):  ref base A (38/35);  new base G (48/29);  total (86/65)
Fisher's exact test for biased strand distribution p-value = 2.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 151 total aligned reads displayed.

CTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGT  >  NC_000913/658056‑658251
                                                                                                 |                                                                                                  
cTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGCCGCATCGCTGTTTGCCCTGCCGGAAGCAtt                                                                                                  >  2:3283686/1‑100 (MQ=255)
cTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCAtt                                                                                                  >  2:3516694/1‑100 (MQ=255)
  gTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGc                                                                                                >  1:962488/1‑100 (MQ=255)
  gTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGc                                                                                                >  1:1112822/1‑100 (MQ=255)
   ttACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCt                                                                                               <  1:305164/100‑1 (MQ=255)
   ttACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCt                                                                                               <  2:2289496/100‑1 (MQ=255)
    tACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCt                                                                                               >  1:2838325/1‑99 (MQ=255)
     aCATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTgg                                                                                             <  1:269528/100‑1 (MQ=255)
        tCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCGCTGCCGGAAGCATTGCTgg                                                                                             >  2:2767535/1‑97 (MQ=255)
        tCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGgcgc                                                                                          >  1:1257698/1‑100 (MQ=255)
        tCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGgcgc                                                                                          >  2:221016/1‑100 (MQ=255)
         cTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGgcgcg                                                                                         >  1:3481460/1‑100 (MQ=255)
         cTGTGTGGGAAAAGAACGCTGAGATTGGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGCCGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGgcgcg                                                                                         <  2:4095721/100‑1 (MQ=255)
           gtgtgGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGcgcgcg                                                                                       <  2:3478891/100‑1 (MQ=255)
             gtgGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCgat                                                                                     >  1:247412/1‑100 (MQ=255)
              tgGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGACCTGCCGGAAGCATTGCTGGCGCGCgatg                                                                                    >  2:2478231/1‑100 (MQ=255)
               gggAAAAGAACGCTGAGATTTGTGCCTCGTTGGTGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCgatga                                                                                   >  2:3863270/1‑100 (MQ=255)
               gggAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCgatga                                                                                   >  1:745759/1‑100 (MQ=255)
                ggAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCgatgat                                                                                  >  1:442991/1‑100 (MQ=255)
                 gAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATc                                                                                 >  1:3841799/1‑100 (MQ=255)
                 gAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATc                                                                                 >  1:1977301/1‑100 (MQ=255)
                 gAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATc                                                                                 >  1:593439/1‑100 (MQ=255)
                 gAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATc                                                                                 >  1:595725/1‑100 (MQ=255)
                  aaaaGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCa                                                                                <  1:3974523/100‑1 (MQ=255)
                   aaaGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCat                                                                               <  2:2931663/100‑1 (MQ=255)
                    aaGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCatg                                                                              <  1:995953/100‑1 (MQ=255)
                    aaGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCatg                                                                              <  2:419390/100‑1 (MQ=255)
                     aGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCatga                                                                             <  1:1975589/100‑1 (MQ=255)
                      gAACGCTGAGATTTGTGCCTCGTTGATGGCGCCGGCTGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCat                                                                               <  2:4101195/97‑1 (MQ=255)
                      gAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCatgat                                                                            >  2:2099407/1‑100 (MQ=255)
                        aCGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGc                                                                          <  1:3373621/100‑1 (MQ=255)
                        aCGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCACCGCGGTTTGCCCTGCCGGAAGCGGTGCTGGCGCGCGATGATCATGATGc                                                                          >  2:2030397/1‑100 (MQ=255)
                         cggtGAGGTTTGTGCCTCGTTTATGGCGAAGGTGGCGGAAAATGACGCATCGCTGTTTGGCCTTCCGGACGCATTGCTAGCGCGCGATGATCATGATGc                                                                          <  2:2732707/96‑1 (MQ=255)
                         cGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCa                                                                         >  2:1244596/1‑100 (MQ=255)
                         cGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCa                                                                         <  2:104105/100‑1 (MQ=255)
                         cGCTGAGATTTGTGCCGCGTTGATGGCGCCGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGc                                                                          <  2:3045727/99‑1 (MQ=255)
                             gagaTTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATc                                                                     <  1:1288587/100‑1 (MQ=255)
                               gaTTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTa                                                                   >  2:1386264/1‑100 (MQ=255)
                                 tttGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATc                                                                 <  1:1241765/100‑1 (MQ=255)
                                  ttGTGCCTCGTTGATGGCGCAGTCGGCGGAAAACGACGCATCGTTTTTTGCCCTGTCTGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATc                                                                 <  2:1870952/99‑1 (MQ=255)
                                   tgtgCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCgg                                                               <  2:3768356/100‑1 (MQ=255)
                                    gtgCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGGATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCgg                                                               <  2:3767026/99‑1 (MQ=255)
                                     tgCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGtt                                                             <  2:1113477/100‑1 (MQ=255)
                                      gCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTa                                                            <  1:2577847/100‑1 (MQ=255)
                                      gCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTa                                                            >  1:2166172/1‑100 (MQ=255)
                                      gCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTa                                                            >  2:591470/1‑100 (MQ=255)
                                      gCCTCGTTGATGGAGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGCTGATCATGATGCAGATCTATCGGCTa                                                            >  2:2562079/1‑100 (MQ=255)
                                       ccTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTaa                                                           >  2:2184267/1‑100 (MQ=255)
                                         tCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAAt                                                         >  1:1200449/1‑100 (MQ=255)
                                          cGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATc                                                        <  2:3607473/100‑1 (MQ=255)
                                            ttGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCTTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAg                                                      <  2:1567430/100‑1 (MQ=255)
                                            ttGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCan                                                      <  1:1297212/100‑2 (MQ=255)
                                             tGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGc                                                     <  1:202297/100‑1 (MQ=255)
                                             tGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGc                                                     <  2:1200441/100‑1 (MQ=255)
                                               aTGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGcac                                                   <  1:3164873/100‑1 (MQ=255)
                                               aTGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGcac                                                   <  1:29243/100‑1 (MQ=255)
                                                tGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGcaca                                                  >  2:918007/1‑100 (MQ=255)
                                                  gcgcAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAgc                                                <  1:498569/100‑1 (MQ=255)
                                                       ggcggcGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGctgctg                                           <  1:1493563/100‑1 (MQ=255)
                                                       ggcggcGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGctgctg                                           <  1:2283161/100‑1 (MQ=255)
                                                         cggcggAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGa                                         >  1:326366/1‑100 (MQ=255)
                                                          ggcggAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGaa                                        >  1:2621877/1‑100 (MQ=255)
                                                           gcggAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAg                                       >  1:3633492/1‑100 (MQ=255)
                                                            cggAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAgg                                      >  1:3017393/1‑100 (MQ=255)
                                                            cggAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAgg                                      >  2:4061627/1‑100 (MQ=255)
                                                              gAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCg                                    <  1:4107333/100‑1 (MQ=255)
                                                                aaaCGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGaa                                  >  2:2295696/1‑100 (MQ=255)
                                                                  acgacgCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAAtt                                <  2:3668546/100‑1 (MQ=255)
                                                                   cgacgCATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTc                               <  1:862346/100‑1 (MQ=255)
                                                                      cgCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTc                            >  2:2564164/1‑100 (MQ=255)
                                                                        cATCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCgg                          <  2:1574128/100‑1 (MQ=255)
                                                                         aTCGCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCgg                          >  1:1405304/1‑99 (MQ=255)
                                                                          tCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGCTGCGGATCTATCGGTTAAATCAGCACAGCTGCTGGAGGGCGAATTCCTCg                           >  2:2416288/1‑97 (MQ=255)
                                                                          tCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGAc                        <  2:428831/100‑1 (MQ=255)
                                                                            gCTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACtt                      >  1:553187/1‑100 (MQ=255)
                                                                            gCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACtt                      >  1:42961/1‑100 (MQ=255)
                                                                             cTGTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACttt                     >  2:3109215/1‑100 (MQ=255)
                                                                               gTTTGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTAc                   <  2:728977/100‑1 (MQ=255)
                                                                               gTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTAc                   >  2:579581/1‑100 (MQ=255)
                                                                                 ttGCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACgg                 >  2:1507696/1‑100 (MQ=255)
                                                                                   gCCCTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCg               <  2:323607/100‑1 (MQ=255)
                                                                                   gCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCCTGATCCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTGATCGCACTTTACGGCg               >  2:1770608/1‑100 (MQ=255)
                                                                                   gCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCg               >  1:705292/1‑100 (MQ=255)
                                                                                   gCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCg               >  2:4003614/1‑100 (MQ=255)
                                                                                    cccTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCga              >  1:3116608/1‑100 (MQ=255)
                                                                                    cccTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCga              >  1:1326822/1‑100 (MQ=255)
                                                                                     ccTGCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCgag             >  1:1315705/1‑100 (MQ=255)
                                                                                       tGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCGAGaa           >  1:3679518/1‑100 (MQ=255)
                                                                                        gCCGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCGAGaaa          >  1:1871243/1‑100 (MQ=255)
                                                                                         ccGGAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCGAGAAAg         <  1:1605844/100‑1 (MQ=255)
                                                                                         ccGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCGAGAAAg         <  2:4032225/100‑1 (MQ=255)
                                                                                         ccGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCGAGAAAg         <  1:3439468/100‑1 (MQ=255)
                                                                                           ggAAGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCGAGAAAGTa       <  1:1726341/100‑1 (MQ=255)
                                                                                           ggAAGCGTTGCTGGCGCGCGATGATCAGGATGCAGATCTATGGGTTAAATCAGCAAAGCTGGTGGAAGGCGAATTCCTCGGACTTTACGGCGAGAAAGTa       >  2:728007/1‑100 (MQ=255)
                                                                                           ggAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCGAGAAAGTa       >  1:4054857/1‑100 (MQ=255)
                                                                                           ggAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCGAGAAAGTa       <  1:2284983/100‑1 (MQ=255)
                                                                                             aaGCGTTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCGAGAAAGTaaa     >  1:3615065/1‑100 (MQ=255)
                                                                                               gCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCGAGAAAGTAAACg   <  2:2959120/100‑1 (MQ=255)
                                                                                                cATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGt  >  1:1843790/1‑100 (MQ=255)
                                                                                                cATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGt  >  1:412528/1‑100 (MQ=255)
                                                                                                 |                                                                                                  
CTGTTACATCTGTGTGGGAAAAGAACGCTGAGATTTGTGCCTCGTTGATGGCGCAGGCGGCGGAAAACGACGCATCGCTGTTTGCCCTGCCGGAAGCATTGCTGGCGCGCGATGATCATGATGCAGATCTATCGGTTAAATCAGCACAGCTGCTGGAAGGCGAATTCCTCGGACTTTACGGCGAGAAAGTAAACGT  >  NC_000913/658056‑658251

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