Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,627,833 T→C D173D (GAT→GAC nadR → DNA‑binding transcriptional repressor/NMN adenylyltransferase NadR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,627,8330TC94.1% 139.0 / ‑0.4 52D173D (GAT→GACnadRDNA‑binding transcriptional repressor/NMN adenylyltransferase NadR
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (28/20);  total (31/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.43e-01

TTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTGCCGGATCGTCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCA  >  NC_000913/4627586‑4628071
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              
ttGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGgacctg                                                                                                                                                                                                                                             >  1:19079/1‑251 (MQ=255)
  gACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAATGGGATTCAGCCGgacctgat                                                                                                                                                                                                                                           <  1:407397/251‑1 (MQ=255)
  gACGAACTGCATATCATTATGGGTTTGGTCGATGCCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGACGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGATGAGGGCATGGAGCCGTAGCCGCACGGCTGGGATGTGTGGAGCAACGGCATCTAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGgacctgat                                                                                                                                                                                                                                           <  1:394879/251‑1 (MQ=255)
   aCGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGgacctgatc                                                                                                                                                                                                                                          <  2:360037/251‑1 (MQ=255)
       acagCATATCATTATGGGATTTGACGAAACACGCAATAACGCGTTGTATGAAGACAGCGCTACGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGGTATTAAAAACTTTTAAATAACAGAAAAATATGCGCATTCATGCTTTCAACGAAGAGGGAGTGGCGCCGTATCCGCAAGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAACAGGGGATTCAGCCGAACCTGATCTacc                                                                                                                                                                                                                                      >  1:124312/4‑250 (MQ=255)
          gCATATCATTATGGGTTTTGACGATACCCGCGCTCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCt                                                                                                                                                                                                                                   <  2:40445/251‑1 (MQ=255)
          gCATATCATTATGGGTTATGTCGATTCCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTGTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGGGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACAcc                                                                                                                                                                                                                                    <  2:117279/250‑1 (MQ=255)
              atccttATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAACTATCAGAAAACTATTCGCATTCATGCTTTCAACGAAGAGGGCCTGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGCCCTGGTCTCCACCTCggt                                                                                                                                                                                                                               >  2:121861/5‑250 (MQ=255)
               acATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACTTGAGCTACACCTCGgaa                                                                                                                                                                                                                              <  1:345835/250‑1 (MQ=255)
                cATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGga                                                                                                                                                                                                                               >  1:398325/1‑249 (MQ=255)
                 aTTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGCCAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGga                                                                                                                                                                                                                               >  1:279857/1‑248 (MQ=255)
                 aTTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAAACAGCCCTATGTAGCCGAAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGCAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGAAGCGGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCCTACATAAGTTTATGGCTGAAACGGGGATTCAGCCGTAACTGCTCTACCACTCGga                                                                                                                                                                                                                               >  1:427900/1‑248 (MQ=255)
                         ttttGACGATACCCGCGATCGCGCGTTGATTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGATTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGACCAGGGGATTCAGCCGGACCTGGACTACACCTCGGAAGAAACCGATg                                                                                                                                                                                                                    >  2:480659/1‑251 (MQ=255)
                            tGACGATACCCGTGACCGCGCGTAGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTGCCGGATCGTCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACCCCTCGGAAGAAGCCGATgcgc                                                                                                                                                                                                                 >  2:153531/1‑251 (MQ=255)
                                             gcgcgTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAACTATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAACAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAACCGAt                                                                                                                                                                                                                     >  2:415227/1‑230 (MQ=255)
                                             gcgcgTTGGTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGAt                                                                                                                                                                                                                     <  1:415227/230‑1 (MQ=255)
                                              cgcgTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGCTGCGCCACAGTATATGGAACATc                                                                                                                                                                                               >  2:378398/1‑251 (MQ=255)
                                                    gtttgAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGa                                                                                                                                                                                         <  2:124254/247‑1 (MQ=255)
                                                         aaGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGCTGCGCCACAGTATATGGAACATCTGGGGATCgat                                                                                                                                                                                    >  1:165132/1‑250 (MQ=255)
                                                                 gctaTGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGgaaga                                                                                                                                                                                                                            >  2:314674/4‑203 (MQ=255)
                                                                 gctaTGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGgaaga                                                                                                                                                                                                                            <  1:314674/200‑1 (MQ=255)
                                                                         gcagcaGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAACCTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGCGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCgctc                                                                                                                                                                    >  2:141205/1‑248 (MQ=255)
                                                                         acgacgattttccgatctggATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATc                                                                                                                                                                    <  2:222744/233‑1 (MQ=255)
                                                                          cagcagCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGCAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATcc                                                                                                                                                                   >  1:397446/1‑251 (MQ=255)
                                                                             cagcCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAACTATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTac                                                                                                                                                                                                                                       >  2:274692/1‑180 (MQ=255)
                                                                             cagcCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTac                                                                                                                                                                                                                                       <  1:274692/180‑1 (MQ=255)
                                                                              agcCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCAGGGATCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGaaa                                                                                                                                                               <  1:479447/251‑1 (MQ=255)
                                                                                 ccaCCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTACATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGCAGCCGATGCGCCACCGTATATTGACCATATGGGGATCGCTACGGTGATGGTCGATCCGacaagt                                                                                                                                                            >  2:427413/3‑246 (MQ=255)
                                                                                  caccttGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTa                                                                                                                                                           <  1:135070/246‑1 (MQ=255)
                                                                                           ggATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATc                                                                                                                                                                    >  1:222744/1‑233 (MQ=255)
                                                                                                  cTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATAGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGTTCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATc                                                                                                                                             >  2:83939/1‑249 (MQ=255)
                                                                                                    gcgctGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAg                                                                                                                                           <  2:388537/245‑1 (MQ=255)
                                                                                                          gTTATTACAAACTTTTAAATATCAGAAAAATCTTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTgcgc                                                                                                                                   >  2:116136/1‑251 (MQ=255)
                                                                                                            tATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGgaagaa                                                                                                                                                                                                                           <  2:47327/161‑1 (MQ=255)
                                                                                                            tATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGgaagaa                                                                                                                                                                                                                           >  1:47327/1‑161 (MQ=255)
                                                                                                                acAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATGCAGCCCGCCCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATc                                                                                                                              <  2:477519/249‑1 (MQ=255)
                                                                                                                                cAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCAGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAACCGGACCTTTATGTGTATCAGCGGTGCGCAGCTCCGCGAAAACCCGATCCg                                                                                                             >  2:122377/1‑251 (MQ=255)
                                                                                                                                     aaaTATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACt                                                                                                        >  1:80947/1‑251 (MQ=255)
                                                                                                                                       atatTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAGAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTgg                                                                                                      >  1:343130/1‑251 (MQ=255)
                                                                                                                                                 acATGATTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGCAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCTATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGATAACCCGTTCCGCTACTGGGATGATATTc                                                                                            >  2:95368/2‑251 (MQ=255)
                                                                                                                                                        tgcAACGAAGAGGGTATGGAGCCGGATCCGCACTGCGGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATAGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGAGCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCTGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTAcgg                                                                                      <  1:222220/248‑3 (MQ=255)
                                                                                                                                                                           aGCCGTTTCCGCCCGGCTGGGATGGGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTtgta                                                                  <  2:397446/251‑2 (MQ=255)
                                                                                                                                                                                   ccGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGt                                                          >  1:400409/1‑251 (MQ=255)
                                                                                                                                                                                          gCTGGGATGTGTGGAGCAACGGCATAAAAAAGTTTATGGCTGACAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCCGGGGATCGCTACGGAGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAAACCGTTCCGCTAATGGGAATATATTCCTACGGAAGTGAAGACGTTCTTTGTACGTACCGTgg                                                        >  2:73849/1‑246 (MQ=255)
                                                                                                                                                                                                              ggCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCttt                                                                                                                                                      >  1:115674/1‑132 (MQ=255)
                                                                                                                                                                                                              ggCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCttt                                                                                                                                                      <  2:115674/132‑1 (MQ=255)
                                                                                                                                                                                                                cATCAAAAAGGTTGTGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATCTGGAACATCTGGGGATCGATACGGTGCTGGTCGCTCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGAGCAGCCGTTCTTTGGACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCggg                              >  2:208931/1‑249 (MQ=255)
                                                                                                                                                                                                                      aaaGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATcc                                                                                                                                                                   >  1:185700/1‑111 (MQ=255)
                                                                                                                                                                                                                      aaaGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATcc                                                                                                                                                                   <  2:185700/111‑1 (MQ=255)
                                                                                                                                                                                                                       aaGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATacac                      >  2:58279/1‑251 (MQ=255)
                                                                                                                                                                                                                       aaGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCAc                      <  1:153531/251‑1 (MQ=255)
                                                                                                                                                                                                                           ttATGGCTGAAAAGGGGATTCAGCCGGACATGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTCCCTTTATGAGTATCCGCGGTGCGCAGATCCGCGAAAACCCGATCCGCTACTGGGAATATATTCCTACGGAAGGGAAGCCGTTCTTTGTACGTAACGTGGCGGTCCTTGGTGGCGAGTGGAGAGGAAGATCCACCCTg                  >  2:472683/1‑251 (MQ=255)
                                                                                                                                                                                                                             aTGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGt                <  1:83939/251‑1 (MQ=255)
                                                                                                                                                                                                                                    gaaaaGGGAATTCAGCCGGACCTGCTCTACACCTCGGAAGAAGCCGATGCGCAACAGTATATGGAACATCTGGGGATCGCTACGGTGCTGGTCGATCCGCAAAGTCCCTTTATGAGTATCAGCGGTGCGCAGATCCGCGACAACCCGTTACGCTACAGGGCATATATTCCTACGGACGTGAAGCCGTTCTTTGTACGTACCGTGGCGCTCATTGGTGGCGAGTCGATCGGTAAATCAACCCTGGAaaaca           >  2:412514/2‑250 (MQ=255)
                                                                                                                                                                                                                                       aggggATTCAGCCGGACCGGATCTACTCCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACtt      <  2:398325/249‑1 (MQ=255)
                                                                                                                                                                                                                                            gagtcAGCCTGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCa  <  2:379211/247‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              
TTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTGCCGGATCGTCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCA  >  NC_000913/4627586‑4628071

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