Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,627,837 A→T 38.1% I175F (ATC→TTC)  nadR → trifunctional protein: nicotinamide mononucleotide adenylyltransferase, ribosylnicotinamide kinase, transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,627,8370AT38.1% 88.3 / 112.0 103I175F (ATC→TTC) nadRtrifunctional protein: nicotinamide mononucleotide adenylyltransferase, ribosylnicotinamide kinase, transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (23/40);  new base T (21/18);  total (44/59)
Fisher's exact test for biased strand distribution p-value = 1.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

Only 100 of 120 total aligned reads displayed.

TTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCG  >  NC_000913/4627738‑4627936
                                                                                                   |                                                                                                   
ttCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGgatctga                                                                                                     >  2:942182/1‑100 (MQ=255)
ttCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGgatctga                                                                                                     >  2:2544662/1‑100 (MQ=255)
  cAATGAAGAGGGCATGGAGCCATATCCGCACGGCTGGGATGTATGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATCCAGCCGgacctgatc                                                                                                   <  1:1168265/100‑1 (MQ=255)
  cAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGgatctgttc                                                                                                   >  1:2915537/1‑100 (MQ=255)
  cAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGgatctgatc                                                                                                   >  1:1182819/1‑100 (MQ=255)
  cAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGgatctgatc                                                                                                   >  2:1662002/1‑100 (MQ=255)
   aaCGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGgatctgatct                                                                                                  >  2:1262847/1‑100 (MQ=255)
    aCGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGTCGgatctgttcaa                                                                                                 <  1:2062957/100‑3 (MQ=255)
     cGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTac                                                                                                <  1:2024099/100‑1 (MQ=255)
      gAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTaca                                                                                               <  2:901964/100‑1 (MQ=255)
        agagGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGTAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACAcc                                                                                             <  1:449480/100‑1 (MQ=255)
        agagGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACAcc                                                                                             <  1:1611206/100‑1 (MQ=255)
            ggCATGGAGCCGTATCCGCACGTTTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAATGGGATCGCGCCGGATGTGATCTACACCTCgg                                                                                         <  2:713011/100‑1 (MQ=255)
            ggCATGGAGCCGTATCCGCACGGCTGGGATTTGTTGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCgg                                                                                         <  1:814303/100‑1 (MQ=255)
            ggCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCgg                                                                                         <  1:1749219/100‑1 (MQ=255)
             gCATGGAGCCGTGACCGCACGGCTGAGATGTGTGGAGCAACGGCTTCTAAACGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTACTGCACCTCGga                                                                                        >  1:3032461/1‑100 (MQ=255)
             gCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGga                                                                                        <  1:188296/100‑1 (MQ=255)
              cATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGgaa                                                                                       >  1:80890/1‑100 (MQ=255)
              cATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGgaa                                                                                       <  2:1260189/100‑1 (MQ=255)
               aTGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGTATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGgaag                                                                                      <  1:1010732/100‑1 (MQ=255)
                tGGAGCAGTATCCGCACGGCTGGGATGGGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGgaag                                                                                      <  2:787460/99‑1 (MQ=255)
                 ggAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGgaagaa                                                                                    >  2:1006729/1‑100 (MQ=255)
                  gAGCCGTATCCGCACGGCGGGGATGTGTGGAGCAACGGCATCAAAAAGTATATGAATGAAAAAGGGATCCAGCCGGATCTGCTCTACACCTCGga                                                                                        <  2:1879305/95‑1 (MQ=255)
                    gCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGTATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGcc                                                                                 <  1:2130664/100‑1 (MQ=255)
                    gCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGcc                                                                                 <  1:761419/100‑1 (MQ=255)
                      cGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGa                                                                               >  1:2595748/1‑100 (MQ=255)
                           ccGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATgcgc                                                                          >  2:424832/1‑100 (MQ=255)
                            cGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGTTGAAAAAGGGACCCAGCCCGATCTGATCTACACCTCGGAAGAAGCCGATGCGcc                                                                         <  2:2289220/100‑1 (MQ=255)
                            cGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATgcgc                                                                          >  1:411695/1‑99 (MQ=255)
                            cGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGCATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGcc                                                                         >  1:1361041/1‑100 (MQ=255)
                             gCACGGCTGGTATGTGTAGAAAAAAGAAATAAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGcc                                                                         <  2:1380016/99‑1 (MQ=255)
                                cGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAg                                                                     <  2:3032428/100‑1 (MQ=255)
                                cGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAg                                                                     <  2:970072/100‑1 (MQ=255)
                                    tGGGATGTGTGGAGCAACGGTATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGtat                                                                  <  1:753666/99‑1 (MQ=255)
                                     gggATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTatat                                                                <  1:993194/100‑1 (MQ=255)
                                      ggATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATg                                                               >  1:1971534/1‑100 (MQ=255)
                                        aTGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGa                                                             <  1:2731358/100‑1 (MQ=255)
                                         tgtgtgGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATggga                                                            <  1:2053439/100‑3 (MQ=255)
                                         tgtgtgGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGaa                                                            <  1:527140/100‑1 (MQ=255)
                                          gtgtgGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAAc                                                           >  1:2308029/1‑100 (MQ=255)
                                          gtgtgGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTAAACTACACCTCGGAAGAAGCCAATTCGCCACAGTATATGGAAc                                                           >  2:248268/1‑100 (MQ=255)
                                           tgtgGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACa                                                          <  1:2478036/100‑1 (MQ=255)
                                           tgtgGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACa                                                          >  2:1465904/1‑100 (MQ=255)
                                             tgGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATc                                                        >  2:1488257/1‑100 (MQ=255)
                                              ggAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCt                                                       >  2:2335766/1‑100 (MQ=255)
                                              ggAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCt                                                       >  1:1640611/1‑100 (MQ=255)
                                               gAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTg                                                      <  1:1672954/100‑1 (MQ=255)
                                                aGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTgg                                                     >  1:2384246/1‑100 (MQ=255)
                                                 gCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTggg                                                    >  2:1897946/1‑100 (MQ=255)
                                                 gCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTggg                                                    <  1:1262859/100‑1 (MQ=255)
                                                  cAACGGCATAAAAAAGTTTATGGCTGAAAAAGGTATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTgggg                                                   <  1:2284713/100‑1 (MQ=255)
                                                   aaCGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGa                                                  <  2:2330613/100‑1 (MQ=255)
                                                   aaCGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGa                                                  >  1:2922425/1‑100 (MQ=255)
                                                   aaCGGCATCAAAAAGTTTATGGCTGAAAAAAGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGACCagatcggaa                                                  >  2:1994418/1‑92 (MQ=255)
                                                    aCGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGAt                                                 >  1:1943142/1‑100 (MQ=255)
                                                     cGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACATAGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATc                                                >  1:1041190/1‑100 (MQ=255)
                                                      ggCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACATAGGAAGAAGACGATGCGCCACAGTATATGGAACATCTGGGGATCg                                               >  1:143363/1‑100 (MQ=255)
                                                        cATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCgag                                             <  2:1557022/100‑1 (MQ=255)
                                                         aTCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTAGGGATCgaga                                            >  1:1414021/1‑100 (MQ=255)
                                                          tgaaaaaGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGAc                                           <  1:2968634/98‑1 (MQ=255)
                                                             aaaaGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGt                                        <  2:2167542/100‑1 (MQ=255)
                                                              aaaGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTg                                       >  1:2461477/1‑100 (MQ=255)
                                                               aaGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGc                                      <  1:2544693/100‑1 (MQ=255)
                                                                aGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCt                                     <  2:105106/100‑1 (MQ=255)
                                                                 gTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGc                                      >  2:830653/1‑98 (MQ=255)
                                                                 gTTTATGGCTGAAAAAGGGATCCAGCCGGATATGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTg                                    >  2:403657/1‑100 (MQ=255)
                                                                  tttATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTgg                                   <  2:2915505/100‑1 (MQ=255)
                                                                    tATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTc                                 <  1:2837763/100‑1 (MQ=255)
                                                                    tATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTc                                 <  1:2606922/100‑1 (MQ=255)
                                                                       ggCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGAt                              >  2:177892/1‑100 (MQ=255)
                                                                        gCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATc                             <  2:2714714/100‑1 (MQ=255)
                                                                          tGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCg                           <  1:2491022/100‑1 (MQ=255)
                                                                          tGAAAAAGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCg                           <  1:690732/100‑1 (MQ=255)
                                                                           gAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCAGGGGATCGAGACGGTGCTGATCGATCCGa                          >  1:869720/1‑100 (MQ=255)
                                                                             agaacagATCCGGCTGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATAGAAAATCTAGGAATCGAGACGGTGCTGGTCGATCCGaaa                        >  1:1875929/7‑100 (MQ=255)
                                                                              aaaGGGATCCAGCCGGATCTTTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAAc                       <  1:2066010/100‑1 (MQ=255)
                                                                              aaaGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAAc                       <  1:2520769/100‑1 (MQ=255)
                                                                               aacagATCCGGCTGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACg                      >  2:79903/5‑100 (MQ=255)
                                                                               aaGGGATCCAGCCAGATCTGTTCTACACCTCGGAATAAGCCGATGCGCCACAGTATATTGAACATCTGGGGATCGATACGGTGCTCGTCGATCCAAAACg                      <  2:2528857/100‑1 (MQ=255)
                                                                                acagATCCGGCTGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATATGGGGATAGAGACGGTGCTGGTAGATACGAAACGt                     >  1:1970413/4‑100 (MQ=255)
                                                                                aGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGt                     >  1:46438/1‑100 (MQ=255)
                                                                                aGGGATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGt                     <  1:2064973/100‑1 (MQ=255)
                                                                                 gggATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTa                    <  2:2300681/100‑1 (MQ=255)
                                                                                  ggATCCAGCCGGATCTGATCTACTCCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTAc                   <  1:942188/100‑1 (MQ=255)
                                                                                   gATCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTAcc                  >  1:2742954/1‑100 (MQ=255)
                                                                                    aTCCAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTAcc                  >  1:1290976/1‑99 (MQ=255)
                                                                                     tCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCTCAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCtt                <  2:1943119/100‑1 (MQ=255)
                                                                                       cAGCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTa              >  2:116380/1‑100 (MQ=255)
                                                                                       cAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTa              >  2:2907056/1‑100 (MQ=255)
                                                                                         gCCGGATCTGTTCTACACCTCGTAAGACGCCGGTGCGCGACATTATATGGAAAATCTGTGTATCGATTTTGTTCTGGTCGATCCGAAACTTACCTTTATg            <  2:665189/100‑1 (MQ=255)
                                                                                         gCCGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATg            <  1:1568477/100‑1 (MQ=255)
                                                                                         gCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATg            <  1:1488273/100‑1 (MQ=255)
                                                                                          ccGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCttt               <  2:2618248/96‑1 (MQ=255)
                                                                                           cGGATCTGTTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAg          >  2:2314661/1‑100 (MQ=255)
                                                                                             gatctgttctACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTa        <  1:2133166/100‑1 (MQ=255)
                                                                                             gatctgttctACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATATGGGGATCGAGACGGTGCTAGTCGATCCGAAACGTACCTTTATGAGTa        >  1:2259239/1‑100 (MQ=255)
                                                                                             gatctgatctACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTa        >  2:912775/1‑100 (MQ=255)
                                                                                                 tgatctACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAg    >  1:2966464/1‑100 (MQ=255)
                                                                                                  gatctACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGc   <  2:791138/100‑1 (MQ=255)
                                                                                                  gatctACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGc   <  1:987650/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
TTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCG  >  NC_000913/4627738‑4627936

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: