Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,628,114 A→C Y267S (TAT→TCT)  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,628,1140AC96.7% 629.6 / 5.0 185Y267S (TAT→TCT) nadRtrifunctional protein: nicotinamide mononucleotide adenylyltransferase, ribosylnicotinamide kinase, transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (3/3);  new base C (78/100);  total (82/103)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 185 total aligned reads displayed.

GCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGG  >  NC_000913/4628017‑4628209
                                                                                                 |                                                                                               
gCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTg                                                                                               <  2:2031321/100‑1 (MQ=255)
  gATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTc                                                                                             >  1:5286646/1‑100 (MQ=255)
  gATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTc                                                                                             >  2:4373652/1‑100 (MQ=255)
    tCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCtt                                                                                           <  1:102106/100‑1 (MQ=255)
    tCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCtt                                                                                           >  1:3126439/1‑100 (MQ=255)
      cTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCtttt                                                                                         >  1:2276535/1‑100 (MQ=255)
      cTTGGCGGCGAGTCGAGCGGTAAATCCACCCTCGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCtttt                                                                                         >  2:4944302/1‑100 (MQ=255)
      cTTGGCGCCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATGTTCAACACCACCAGTGCGTGGGAATATGGCCGCGGTTCTGTCttgt                                                                                         <  2:5103973/100‑3 (MQ=255)
       ttGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTc                                                                                        <  1:441544/100‑1 (MQ=255)
       ttGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTc                                                                                        <  2:1715776/100‑1 (MQ=255)
         ggcggcGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTcac                                                                                      <  2:5214249/100‑1 (MQ=255)
          gcggcgAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTcaca                                                                                     <  2:486565/100‑1 (MQ=255)
                 gtagAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGCATATGGCCGCGCTTCTGTCTTTTCACACCTcggc                                                                              >  2:2530406/4‑100 (MQ=255)
                 gTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTcggc                                                                              >  1:5054653/1‑100 (MQ=255)
                   cGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTcggcgg                                                                            <  2:2404621/100‑1 (MQ=255)
                   cGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTcggcgg                                                                            >  2:1365913/1‑100 (MQ=255)
                    gAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGt                                                                           <  1:555141/100‑1 (MQ=255)
                    gAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGt                                                                           <  1:250491/100‑1 (MQ=255)
                     aGCGGTAAATCCACCCTGGTAAACAAACTTGCCCATGTCTTCAACGCCGCCAGTGCGTGGGAACAGGGCGGCGACTCTGTCTTTTCACACCTCGGCGGtg                                                                          <  2:4158266/100‑1 (MQ=255)
                     aGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGtg                                                                          <  1:5207173/100‑1 (MQ=255)
                      gCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGtga                                                                         <  1:351617/100‑1 (MQ=255)
                       cGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGtgat                                                                        <  2:1658551/100‑1 (MQ=255)
                        ggTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGtgatg                                                                       >  1:5487632/1‑100 (MQ=255)
                        ggTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGtgatg                                                                       >  1:762205/1‑100 (MQ=255)
                        ggTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGtgatg                                                                       <  1:4020418/100‑1 (MQ=255)
                         gTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGtgatga                                                                      <  2:2725410/100‑1 (MQ=255)
                           aaaTCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATgaga                                                                    <  1:425855/100‑1 (MQ=255)
                            aaTCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGAt                                                                   >  2:772044/1‑100 (MQ=255)
                            aaTCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGAt                                                                   >  1:3862763/1‑100 (MQ=255)
                            aaTCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGAt                                                                   >  1:769677/1‑100 (MQ=255)
                              tCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCg                                                                 <  2:417755/100‑1 (MQ=255)
                              tCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCg                                                                 <  2:3463136/100‑1 (MQ=255)
                               ccACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGc                                                                >  2:3992097/1‑100 (MQ=255)
                               ccACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGc                                                                >  2:3183380/1‑100 (MQ=255)
                                cACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCa                                                               >  1:441535/1‑100 (MQ=255)
                                 aCCCTGGTAAACAAACTTGCCAATATCTTCAACACCGCCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACTCCTCGGCGGTGATGAGATCGCAt                                                              <  2:5466129/100‑1 (MQ=255)
                                 aCCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCAt                                                              <  1:4644278/100‑1 (MQ=255)
                                  cccTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCAtt                                                             >  1:387247/1‑100 (MQ=255)
                                   ccTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTg                                                            <  1:105338/100‑1 (MQ=255)
                                    cTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGc                                                           <  1:3269399/100‑1 (MQ=255)
                                    cTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGc                                                           <  1:1294934/100‑1 (MQ=255)
                                      ggTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAg                                                         >  2:4217152/1‑100 (MQ=255)
                                           acaaacTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTc                                                    >  2:2206965/1‑100 (MQ=255)
                                           acaaacTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTc                                                    >  1:1365138/1‑100 (MQ=255)
                                            caaacTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCt                                                   <  1:2555554/100‑1 (MQ=255)
                                             aaacTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTg                                                  <  2:1036680/100‑1 (MQ=255)
                                                cTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACt                                               <  2:3420729/100‑1 (MQ=255)
                                                  tGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTAc                                             >  1:1015534/1‑100 (MQ=255)
                                                   gCCAATCTCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACg                                            <  2:5329414/100‑1 (MQ=255)
                                                   gCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACg                                            >  2:5051081/1‑100 (MQ=255)
                                                   gCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACg                                            <  1:3720141/100‑1 (MQ=255)
                                                      atcttCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATa                                         <  2:5502145/96‑1 (MQ=255)
                                                      aaTATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATa                                         >  1:345842/1‑100 (MQ=255)
                                                      aaTATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATa                                         <  2:116358/100‑1 (MQ=255)
                                                      aaTATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATa                                         >  2:2278984/1‑100 (MQ=255)
                                                      aaTATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATa                                         >  1:2376175/1‑100 (MQ=255)
                                                         atCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATaaaa                                      >  2:3258578/1‑100 (MQ=255)
                                                          tCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAAt                                     <  1:3779995/100‑1 (MQ=255)
                                                           cTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATc                                    >  2:2981842/1‑100 (MQ=255)
                                                           cTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATc                                    >  2:3528354/1‑100 (MQ=255)
                                                           cTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATc                                    >  2:2010680/1‑100 (MQ=255)
                                                            ttCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACGCCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATcg                                   <  1:3956078/100‑1 (MQ=255)
                                                            ttCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATcg                                   >  1:902064/1‑100 (MQ=255)
                                                             tCAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATcgc                                  <  1:2578263/100‑1 (MQ=255)
                                                              cAACACCACCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATcgcg                                 <  1:749889/100‑1 (MQ=255)
                                                                 caccacCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTg                              >  2:2927826/1‑100 (MQ=255)
                                                                 caccacCAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTg                              >  1:2981835/1‑100 (MQ=255)
                                                                  accaccAGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTgg                             >  1:196722/1‑100 (MQ=255)
                                                                   ccaccaGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTggg                            >  2:2054479/1‑100 (MQ=255)
                                                                   ccaccaGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTggg                            <  1:4207602/100‑1 (MQ=255)
                                                                   ccaccaGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTggg                            >  2:4383648/1‑100 (MQ=255)
                                                                     accaGTGCTTGGGAATCTGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGcc                          <  1:4577083/100‑1 (MQ=255)
                                                                     accaGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGcc                          <  2:4845790/100‑1 (MQ=255)
                                                                     accaGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGcc                          <  1:3556160/100‑1 (MQ=255)
                                                                        aGTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACg                       >  1:2218779/1‑100 (MQ=255)
                                                                         gTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGc                      >  2:613954/1‑100 (MQ=255)
                                                                         gTGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAATATTCTGACTACGATAAAATCTCGCTGGGCCACGc                      >  2:1320382/1‑100 (MQ=255)
                                                                         gTGCATGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGc                      >  2:4650785/1‑100 (MQ=255)
                                                                          tGCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCt                     <  1:3824385/100‑1 (MQ=255)
                                                                           gCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTc                    >  1:4667544/1‑100 (MQ=255)
                                                                           gCGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTc                    <  2:3082955/100‑1 (MQ=255)
                                                                            cGTGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCa                   >  2:834661/1‑100 (MQ=255)
                                                                              tGGGAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAAt                 >  2:2544637/1‑100 (MQ=255)
                                                                               gggAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATa                <  1:5355478/100‑1 (MQ=255)
                                                                                ggAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATAc               >  2:223858/1‑100 (MQ=255)
                                                                                ggAATATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATAc               <  1:2076561/100‑1 (MQ=255)
                                                                                  aaTATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACAt             >  1:4978960/1‑100 (MQ=255)
                                                                                  aaTATGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACAt             <  1:43717/100‑1 (MQ=255)
                                                                                    tctGGTCGCGTTTATGTCTTTTCACTACCCTGGGGTGATGATATCGCACTGCATTATTCTTACTACGATAAAATCCCGCTGGGCCACGCTCAATACATTg           <  2:1303898/98‑1 (MQ=255)
                                                                                     atGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGa          <  2:800900/100‑1 (MQ=255)
                                                                                      tGGCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGAt         <  1:2743639/100‑1 (MQ=255)
                                                                                       ggCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGAtt        >  2:1821984/1‑100 (MQ=255)
                                                                                       ggCCGCGATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGAtt        >  1:1881519/1‑100 (MQ=255)
                                                                                       ggCCGAGATTCCGTCTGTTCACCCCGCGGAGGGGATGAGATCGCATTTCAGTATTCTGACTACGATTAAATCGCGCTGGGCCACGCTCAATACATTGAtt        <  2:4547597/100‑1 (MQ=255)
                                                                                          cgcgATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTg     <  1:3009981/100‑1 (MQ=255)
                                                                                           gcgATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGc    <  1:2444523/100‑1 (MQ=255)
                                                                                            cgATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCg   >  1:683065/1‑100 (MQ=255)
                                                                                            cgATTCTGTCTTTTCACACCTCGGCGGAGATGAGATCGCATTGCAGTATTCTGACTACGATAAACTCGCGCTGGGCCACGCTCAATACATTGATTTTGCg   >  2:5414215/1‑100 (MQ=255)
                                                                                             gATTCTGTCTTTTCACCCCTCGGCGGTGATCAGATCGCACTGCAGTATTCTGACTACGATACAATCGCGCTGAGCCACGCTCCCTACATTGCTTTTGccg  <  2:1122218/100‑3 (MQ=255)
                                                                                             gATTCTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCgg  <  1:796084/100‑1 (MQ=255)
                                                                                                 |                                                                                               
GCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGG  >  NC_000913/4628017‑4628209

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