Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,627,689 T→C R125R (CGT→CGC 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,6890TC97.9% 148.4 / ‑4.3 48R125R (CGT→CGCnadRDNA‑binding transcriptional repressor/NMN adenylyltransferase NadR
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (24/23);  total (25/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.49e-01

AGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTGCCGGATCGTCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATC  >  NC_000913/4627444‑4627932
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   
aGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGtt                                                                                                                                                                                                                                                <  1:364009/251‑1 (MQ=255)
    cgcAAAAGCTGTAGGCTTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCTGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTAtta                                                                                                                                                                                                                                            <  2:412184/251‑2 (MQ=255)
      cAAAAGCTTGAGGCGTTGCACCGTTTGTTGGGGTTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCAGACCGTACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATAGCATTATGGGTTTTGACGATACCCGCGAGCGCGCGGTGTTTGAAGACAGCGCTATGGCGCAGCAGCCCACCGTGCCGGAGCGTCTGCGCTGGTTAttaca                                                                                                                                                                                                                                          <  2:182654/251‑4 (MQ=255)
        aaaGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACaaa                                                                                                                                                                                                                                        <  1:133570/251‑1 (MQ=255)
                            gTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAat                                                                                                                                                                                                                    >  1:274716/1‑251 (MQ=255)
                                   ggggCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGc                                                                                                                                                                                                              <  2:274716/250‑1 (MQ=255)
                                           aGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCt                                                                                                                                                                                                     <  2:455328/251‑1 (MQ=255)
                                                               accattggTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATgg                                                                                                                                                                                  >  1:388537/7‑250 (MQ=255)
                                                                                cGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAata                                                                                                                                                                                                                   <  1:133615/200‑1 (MQ=255)
                                                                                cGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAata                                                                                                                                                                                                                   >  2:133615/1‑200 (MQ=255)
                                                                                       ttCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGCACTGCATATCATTATGGGATTTGACGCTACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATGGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATACGCCCGGCTGGGAtgg                                                                                                                                                         >  1:466039/1‑250 (MQ=255)
                                                                                           aCCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGTCGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGAtgtgtg                                                                                                                                                      <  2:441660/250‑1 (MQ=255)
                                                                                                                         tCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGa                                                                                                                        >  1:124254/1‑250 (MQ=255)
                                                                                                                                              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)
                                                                                                                                                                                                 tgttTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGc                                                                                                                                                                                                      >  1:445933/1‑100 (MQ=255)
                                                                                                                                                                                                 tgttTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGc                                                                                                                                                                                                      <  2:445933/100‑1 (MQ=255)
                                                                                                                                                                                                  gtttgAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGa                                              <  2:124254/247‑1 (MQ=255)
                                                                                                                                                                                                      gAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGcc                                                                                                         <  1:337811/188‑1 (MQ=255)
                                                                                                                                                                                                      gAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGcc                                                                                                         >  2:337811/1‑188 (MQ=255)
                                                                                                                                                                                                       aaGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGCTGCGCCACAGTATATGGAACATCTGGGGATCgat                                         >  1:165132/1‑250 (MQ=255)
                                                                                                                                                                                                               gctaTGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGgaaga                                                                                 <  1:314674/200‑1 (MQ=255)
                                                                                                                                                                                                               gctaTGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGgaaga                                                                                 >  2:314674/4‑203 (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)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   
AGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTGCCGGATCGTCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATC  >  NC_000913/4627444‑4627932

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