Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,627,774 T→A 30.1% W154R (TGG→AGG)  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,7740TA30.1% 250.7 / 150.2 178W154R (TGG→AGG) nadRtrifunctional protein: nicotinamide mononucleotide adenylyltransferase, ribosylnicotinamide kinase, transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (47/76);  new base A (30/23);  total (78/100)
Fisher's exact test for biased strand distribution p-value = 3.12e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.62e-01

Only 100 of 182 total aligned reads displayed.

CCGGATCGTCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTA  >  NC_000913/4627675‑4627871
                                                                                                   |                                                                                                 
ccGGATCGTCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCt                                                                                                   >  2:13620/1‑100 (MQ=255)
     tCGTCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGAt                                                                                              >  1:3492766/1‑100 (MQ=255)
      cGTCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGTCTGGGAtg                                                                                             >  2:2765107/1‑100 (MQ=255)
      cGTCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGAtg                                                                                             <  2:4151683/100‑1 (MQ=255)
      cGTCTGCGTTGGTTATTGCAAACTTGTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGAtg                                                                                             <  2:1707118/100‑1 (MQ=255)
       gTCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCAGGGAtgt                                                                                            >  1:3216973/1‑100 (MQ=255)
         cTGCGTTGGTTATTGCAAACTTGTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGAtgtgt                                                                                          >  2:2329817/1‑100 (MQ=255)
          tGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCAGGGAtgtgt                                                                                          >  1:4063522/1‑99 (MQ=255)
           gCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTgg                                                                                        <  1:1762220/100‑1 (MQ=255)
           gCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCAGGGATGTGTgg                                                                                        <  2:4068027/100‑1 (MQ=255)
           gCGTTGGTTATTGCAAACTTGTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTgg                                                                                        <  1:3105734/100‑1 (MQ=255)
           gCGTTGGTTATTGCAAACTTGTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCACACGGCTGGGATGTGTgg                                                                                        >  2:875054/1‑100 (MQ=255)
            cGTTGGTTATTGCAAACTTGTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGa                                                                                       <  2:1606409/100‑1 (MQ=255)
             gTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAg                                                                                      <  1:3646321/100‑1 (MQ=255)
             gTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCAGGGATGTGTGGAg                                                                                      <  2:2187908/100‑1 (MQ=255)
              ttGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGc                                                                                     <  2:1827917/100‑1 (MQ=255)
              ttGGTTATTGCAAACTTGTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGc                                                                                     <  2:1346558/100‑1 (MQ=255)
                  ttATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCAGGGATGTGTGGAGCAACg                                                                                 <  2:1489733/100‑1 (MQ=255)
                   tATTGCAAACTTGTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACgg                                                                                <  2:1190967/100‑1 (MQ=255)
                    aTTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGc                                                                               <  2:2039625/100‑1 (MQ=255)
                    aTTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCAGGGATGTGTGGAGCAACGGc                                                                               <  1:2092688/100‑1 (MQ=255)
                    aTTGCAAACTTGTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGc                                                                               <  2:896524/100‑1 (MQ=255)
                    aTTGCAAACTTGTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGc                                                                               >  2:1658167/1‑100 (MQ=255)
                     ttGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCAGGGATGTGTGGAGCAACGGCa                                                                              >  1:2413739/1‑100 (MQ=255)
                       gCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATc                                                                            <  1:3881147/100‑1 (MQ=255)
                        cAAACTTGTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCa                                                                           <  2:3337209/100‑1 (MQ=255)
                            cTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCaaaa                                                                        <  1:2410754/99‑1 (MQ=255)
                            cTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCAGGGATGTGTGGAGCAACGGCATCaaaaa                                                                       >  2:3887639/1‑100 (MQ=255)
                            cTTGTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCaaaaa                                                                       >  2:61830/1‑100 (MQ=255)
                              tttAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCAGGGATGTGTGGAGCAACGGCATCAAAAAGt                                                                     >  2:2489681/1‑100 (MQ=255)
                                  aaTATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTAt                                                                 <  1:4108587/100‑1 (MQ=255)
                                    tatCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATgg                                                               <  2:3990472/100‑1 (MQ=255)
                                    tatCAGAAAAATATTCGCATTCATGCTCTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCAGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATgg                                                               <  1:3342322/100‑1 (MQ=255)
                                      tCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCAGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCt                                                             >  2:962535/1‑100 (MQ=255)
                                      tCAGAAAAATATTCGCATTCATACTTTCAACGAAGAGGTCATGGAGCCGTATCCGCAAGGAATGGATGTGTGGAGCAGCGGCATCAAAACGTTTATGGc                                                              >  1:340554/1‑99 (MQ=255)
                                       cAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCAGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTg                                                            >  1:4223941/1‑100 (MQ=255)
                                        aGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGa                                                           >  2:3776019/1‑100 (MQ=255)
                                         gAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGaa                                                          <  2:1738507/100‑1 (MQ=255)
                                            aaaTATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGaaaaa                                                       <  2:4228623/100‑1 (MQ=255)
                                              atatTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAgg                                                     <  2:3648710/100‑1 (MQ=255)
                                              atatTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAgg                                                     <  1:1381182/100‑1 (MQ=255)
                                                 ttCTCATTCATGCTTTGAACGAAGAGGGAAGGGAGCCGTATCCGCACGGCAGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAATGGAt                                                  <  1:1272149/100‑1 (MQ=255)
                                                 ttCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCAGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGAt                                                  <  2:4265930/100‑1 (MQ=255)
                                                 ttCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCAGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGAt                                                  >  1:2393321/1‑100 (MQ=255)
                                                   cGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATcc                                                <  1:1627016/100‑1 (MQ=255)
                                                   cGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATcc                                                <  2:1230555/100‑1 (MQ=255)
                                                   cGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGCCTGGGATGTGTGGAGCAACGGCATCAAAAAGCTTATGGCTGAAAAAGGGATcc                                                >  2:567573/1‑100 (MQ=255)
                                                     cATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCAGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAg                                              >  1:3569451/1‑100 (MQ=255)
                                                      aTTCATGCTTTCAACGAAGCGGGCATGGAGCCGTATCAGCACGGCAGGGATGTGTGGAGCAACGGCATCAAAAAGGTACTGGCTGAAAAAGGGATCCAGc                                             >  1:606960/1‑100 (MQ=255)
                                                      aTTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGc                                             >  2:3523642/1‑100 (MQ=255)
                                                       ttCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCAGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGcc                                            >  2:2062018/1‑100 (MQ=255)
                                                        tCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCAGGGATGTGTGGAGCAACGGCATCAAAAAGAATAAAAAAAAAAAAGGGATCCAGCCg                                           <  1:1558706/100‑1 (MQ=255)
                                                         cATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCgg                                          >  1:3194597/1‑100 (MQ=255)
                                                         cATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCgg                                          >  1:3188053/1‑100 (MQ=255)
                                                            gCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGgatc                                       <  2:2640453/100‑1 (MQ=255)
                                                            gCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGgatc                                       <  1:271531/100‑1 (MQ=255)
                                                               ttCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGGGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGgatctga                                    >  1:1331387/1‑100 (MQ=255)
                                                               ttCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCAGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGgatctga                                    >  2:613045/1‑100 (MQ=255)
                                                                tCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGgatctgat                                   <  1:2978205/100‑1 (MQ=255)
                                                                tCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCAGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGgatctgat                                   >  2:3771448/1‑100 (MQ=255)
                                                                 cAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGgatctgttc                                  <  2:3856323/100‑1 (MQ=255)
                                                                 cAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGGGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCAAGCCGgatctgatc                                  >  1:428565/1‑100 (MQ=255)
                                                                  aaCGAAGAGGGCATGGAGTCGTATCCGCGCGGCTGGAATGTGTGGAGCAACGTCATCAAAAAGTTTATAGATAAAAAAGGGATCCGGCTGgatccctt                                   <  1:3511659/98‑5 (MQ=255)
                                                                  aaCGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGgatctgatct                                 <  2:2223534/100‑1 (MQ=255)
                                                                  aaCGAAGAGGGCATGGAGCCGTATCCGCACCGCTGGGATTTGTGGAGCGACGGCATCTAAAAAAAAAAAAAAAAAAAAGGGATCCAGCCGgatctgatc                                  <  1:2763254/99‑1 (MQ=255)
                                                                   aCGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTa                                >  2:3097841/1‑100 (MQ=255)
                                                                   aCGAAGAGGGCATGGAGCCGTATCCGCACGGCAGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTa                                >  1:2766348/1‑100 (MQ=255)
                                                                    cGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCGGCTGgatccctttttc                               <  2:3121223/100‑9 (MQ=255)
                                                                    cGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTac                               >  2:3165312/1‑100 (MQ=255)
                                                                    cGAAGAGGGCATGGAGCCGTATCCGCACGGCCGGGATGTGTGGCGCAACGGCATCAAAAAGTTTAAGGCTGAAACAGGGATCCAGCCGGACATGATCTac                               >  1:299657/1‑100 (MQ=255)
                                                                    cGAAGAGGGCATAGAGCCGTATCCGCACGACAGGGATGTGTGGAGCAACGGCATCAAAAAGAAAAAAAAAAAAAAAGGGATCCAGCCGGATCTGATCTac                               <  1:1469413/100‑1 (MQ=255)
                                                                     gAAGAGGGCATGGAGCCGTATCCGCACGGCAGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTaca                              <  1:903712/100‑1 (MQ=255)
                                                                       agagGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACAcc                            <  2:3358532/100‑1 (MQ=255)
                                                                       agagGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGTTCTACAcc                            <  2:3195072/100‑1 (MQ=255)
                                                                          gggCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTATATAGATAAAAAAGGGATCCAGCCGGATCTGATCTACACCTCg                         <  1:2580730/100‑1 (MQ=255)
                                                                           ggCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCgg                        >  2:3455325/1‑100 (MQ=255)
                                                                           ggCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCgg                        <  2:3819211/100‑1 (MQ=255)
                                                                           ggCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGAGACCCAGCCCGCTCTGATCTACACCTCgg                        <  1:3912741/100‑1 (MQ=255)
                                                                              aTGGAGCCGTATCCGCACGGCAGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGgaag                     <  2:2872642/100‑1 (MQ=255)
                                                                               tGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGgaaga                    <  2:3119519/100‑1 (MQ=255)
                                                                               tGGAGCCGTATCCGCACGGCAGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGgaaga                    >  1:1428478/1‑100 (MQ=255)
                                                                                 gAGCCGTATCCGCACGGCAGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGaagaag                  >  1:2827412/1‑100 (MQ=255)
                                                                                  agtcGTATCCGCACAGCTGAGATGTGTGGAGCAACGGCATCAAAAAGAAAAAAAATAAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGaagaag                  <  1:1843720/96‑1 (MQ=255)
                                                                                  aGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGc                 >  2:3607947/1‑100 (MQ=255)
                                                                                   gCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGcc                >  1:3153419/1‑100 (MQ=255)
                                                                                      gTATCCGCACGGCAGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGAt             >  1:107470/1‑100 (MQ=255)
                                                                                         tCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATgcg          <  1:3607933/100‑1 (MQ=255)
                                                                                         tCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATgcg          <  2:3268632/100‑1 (MQ=255)
                                                                                           cGCACGGCAGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGcc        >  1:3724375/1‑100 (MQ=255)
                                                                                             cACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCcac      <  2:942513/100‑1 (MQ=255)
                                                                                              aCGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCcaca     <  2:3354626/100‑1 (MQ=255)
                                                                                               cGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAg    >  2:2653445/1‑100 (MQ=255)
                                                                                               cGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAg    >  2:1058911/1‑100 (MQ=255)
                                                                                               cGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAg    <  1:1152067/100‑1 (MQ=255)
                                                                                               cGGCAGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATTTTTTTTAATGATACGGCGACCACCGAGATCTACACCTCGGAAGAAGCCGATGCGCCACAg    <  1:1485059/100‑1 (MQ=255)
                                                                                                ggCTGGGATGTGTGGAGCAACGGCATCAAAAAGAATATAACTAAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAg    <  1:1701190/99‑1 (MQ=255)
                                                                                                 gCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGta  <  2:19857/100‑1 (MQ=255)
                                                                                                 gCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGta  <  2:1946919/100‑1 (MQ=255)
                                                                                                 gCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGAGCTACACCTCGGAAGAAGCCGATGCGCCACAGta  >  1:557418/1‑100 (MQ=255)
                                                                                                 gCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGACGGCGCCGGATGAGATCTACACCTCGGAAGAAGCCGATGCGCCACAGta  <  1:2267153/100‑1 (MQ=255)
                                                                                                   |                                                                                                 
CCGGATCGTCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTA  >  NC_000913/4627675‑4627871

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: