Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,570,550 C→T S342N (AGC→AAC)  yjiR ← fused putative DNA‑binding transcriptional regulator/putative aminotransferase YjiR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,570,5500CT94.7% 117.5 / ‑5.8 39S342N (AGC→AAC) yjiRfused putative DNA‑binding transcriptional regulator/putative aminotransferase YjiR
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (20/16);  total (22/17)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

GCTGCCGCGCGACGCAGACCATATCGACCTGTTCAGGCAATTCTATCCACAGTAAAAATCCGCCTTTCGGGCGCGTAATACAGATTTCGCAGGGAAAATATTCCCGTATCCAGCAGGTATAAAGCGCCAAATTGCGCTGATAGATCTGCCGCATCCGCCGGATATGGCGATGATAGTGACCTTCCAGCACAAACGTTGCCGCCGCCATTTGCGTGGACGGCACATTAAAGCTGCTGATGGCGTATTTCATATGCATCAGTTTATCGTGATAACGCCCCGGTGCGACCCAACCCACGCGCAGGCCTGGTGCAATACTTTTACTGAACGAGCTGCACAACAGCACTCGCCCGTCGATATCCCAGGAATGAATGGTCCGCGGGCGCGGATACTCCGTCGCCAGTTCGCCATAGACATCATCTTCAAAAATCACAATATCATGACGCTGAGCGAGAGAGAGAACGGCCCGTTTGCGCGCGT  >  NC_000913/4570320‑4570796
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
gCTGCCGCGCGACGCAGACCAGTCCGACCTGTTCAGGCAATGCTATCCACAGTAAAAATCCGCCTTTCGGTCGCGTAATTCAGATTTCGCAGGGAAAAGATTCCCGTATCCAGCAGGTATAAAGCCCTAAATTGCGCGGATAGATCTGCCGCATCCGCCGGAAATGGCGATGATAGTGACCTTCCAGCACAAACGTTACCGCCGCCATTTGCGTGGACGTCGCATTAAAGGTGCTGATTGCGTATTTCata                                                                                                                                                                                                                                    <  2:390680/251‑1 (MQ=255)
                      atCGACCTGTTCAGGCAATTCTATCCACAGTAAAAATCCGCCTTTCGGGCGCGTAATACAGATTTCGCAGGGAAAATATACCCGTATCCAGCAGGTATAAAGCGCCAAATTGCGCTGATAGATCTGCCGCATCCGCCGGATATGGCGATGATAGTGACCTTCCAGCACAAACGTTGCCGCCGCCATTTGCGTGGACGGCACATTAAAGTTGCTGATGGCGTATTTCATATGCATCAGTTTATCGTGATAAc                                                                                                                                                                                                              <  1:64071/251‑1 (MQ=255)
                              gTTCAGGCAATTCTATCCACAGTAAAAATCCGCCTTTCGGGCGCGTAATACAGATTTCGCAGGGAAAATATTCCCGTATCCAGCAGGTATAAAGCGCCAAATTGCGCTTATAGATCTGCCGCATCCGCCGGATATGGCGATGATAGTGACCTTCCATCACAAACGTTGCCGCCGCCATTTGCGTGGACGGCACATTAAAGTTGCTGATGGCGTATTTCATATGCATCAGTTTATCGTGATAACGCCCCg                                                                                                                                                                                                        <  1:400020/249‑1 (MQ=255)
                                       aTTCTATCCACAGTAAAAATCCGCCTTTCGGGCGCGTAATACAGATTTCGCAGGGAAAATATTCCCGTATCCAGCAGGTATAAAGCGCCAAATTGCGCTGATAGATCTGCCGCATCCGCCGGATATGGCGATGATAGTGACCTTCCAGCACAAACGTTGCCGCCGCCATTTGCGTGGACGGCACATTAAAGTTGCTGATGGCGTATTTCATATGCATCAGTTTATCGTGATAACGCCCCGGTGCGacccaa                                                                                                                                                                                             <  1:83761/251‑1 (MQ=255)
                                                                     ggCGCGTAATACAGATTTCGCAGGGAAAATATTCCCGTATCCAGCAGGTATAAAGCGCCAAATTGCGCTGATAGATCTGCCGCATCCGCCGGATATGGCGATGATAGTGACCTTCCAGCACAAACGTTGCCGCCGCCATTTGCGTGGACGGCACATTAAAGTTGCTGATGGCGTATTTCATATGCATCAGTTTATCGTGATAACGCCCCGGTGCGACCCAACCCACGCGCAGGCCTGGTGCAATACtttt                                                                                                                                                                >  1:361192/1‑250 (MQ=255)
                                                                            aaTACAGATTTCGCAGGGAAAATATTCCCGTATCCAGCAGGTATAAAGCGCCAAATTGCGCTGATAGATCTGCCGCATCCGCCGGATATGGCGATGATAGTGACCTTCCAGCACAAACGTTGCCGCCGCCATTTGCGTGGACGGCACATTAAAGTTGCTGATGGCGTATTTCATATGCATCAGTTTATCGTGATAACGCCCCGGTGCGACCCAACCCACGCGCAGGCCTGGTGCAATACTTTTACTGAACg                                                                                                                                                        >  1:155765/1‑251 (MQ=255)
                                                                             aTACAGATTTCGCAGGGAAAATATTCCCGTATCCAGCAGGTATAAAGCGCCAAATTGCGCTGATAGATCTGCCGCATCCGCCGGATATGGCGATGATAGTGACCTTCCAGCACAAACGTTGCCGCCGCCATTTGCGTGGACGGCACATTAAAGTTGCTGATGGCGTATTTCATATGCATCAGTTTATCGTGATAACGCCCCGGTGCGACCCAACCCACGCGCAGGCCTGGTGCAATACTTTTACTGAACg                                                                                                                                                        >  1:309247/1‑250 (MQ=255)
                                                                                cAGATTTCGCCGGGAAAATCTTCCCGTATCCAGCAGGTATAAAGCGCCAAATTGCGCTGATAGATCTGCCGCATCCGCCGGATATGGCGATGATAGTGACCTTCCAGCACAAACGTTGCCGCCGCCATTTGAGTGGACGGCACATTAAAGTTGCTGATGGCGTATTTCATATGCATCAGTTTATCGTGATAACGCCCCGGTGCGACCCAACCCACGCGCAGGCCTGGTGCAATACTTTTACTGAACGGGCt                                                                                                                                                    >  2:58962/1‑251 (MQ=255)
                                                                                     ttCGCAGGGAAAATATTCCCGTATCCAGCAGGTATAAAGCGCCAAATTGCGCTGATAGATCTGCCGCATCCGCCGGATATGGCGATGATAGTGACCTTCCAGCACAAACGTTGCCGCCGCCATTTGCGTGGACGGCACATTAAAGTTGCTGATGGCGTATTTCATATGCATCAGTTTATCGTGATAACGCCCCGGTGCGACCCAACCCACGCGCAGGCCTGGTGCAATACTTTTACTGAACGAGCTGcac                                                                                                                                                <  1:75986/250‑1 (MQ=255)
                                                                                       cGCAGGGAAAATATTCCCGTATCCAGCAGGTATAAAGCGCCAAATTGCGCTGATAGATCTGCCGCATCCGCCGGATATGGCGATGATAGTGACCTTCCAGCACAAACGTTGCCGCCGCCATTTGCGTGGACGGCACATTAAAGTTGCTGATGGCGTATTTCATATGCATCAGTTTATCGTGATAACGCCCCGGTGCGACCCAACCCACGCGCAGGCCTGGTGCAATACTTTTACTGAACGAGCTGCacaac                                                                                                                                             >  1:115172/1‑251 (MQ=255)
                                                                                       cGCAGGGAAAATATTCCCGTATCCAGCAGGTATAAAGCGCCAAATTGCGCTGATAGATCTGCCGCATCCGCCGGATATGGCGATGATAGTGACCTTCCAGCACAAACGTTGCCGCCGCCATTTGCGTGGACGGCACATTAAAGTTGCTGATGGCGTATTTCATATGCATCAGTTTATCGTGATAACGCCCCGGTGCGACCCAACCCACGCGCAGGCCTGGTGCAATACTTTTACTGAACGAGCTGCacaac                                                                                                                                             <  1:150904/251‑1 (MQ=255)
                                                                                                   atTCCCGTATCCAGCAGGTATAAAGCGCCAAATTGCGCTGATAGATCTGCCGCATCCGCCGGATATGGCGATGATAGTGACCTTCCAGCACAAACGTTGCCGCCGCCATTTGCGTGGACGGCACATTAAAGTTGCTGATGGCGTATTTCATATGCATCAGTTTATCGTGATAACGCCCCGGTGCGACCCAACCCACGCGCAGGCCTGGTGCAATACTTTTACTGAACGAGCTGCACAACAGCACTCGCCCg                                                                                                                                 >  1:111617/1‑251 (MQ=255)
                                                                                                     tCCCGTATCCAGCAGGTATAAAGCGCCAAATTGCGCTGATAGATCTGCCGCATCCGCCGGATATGGCGATGATAGTGACCTTCCAGCACAAACGTTGCCGCCGCCATTTGCGTGGACGGCACATTAAAGTTGCTGAt                                                                                                                                                                                                                                                 >  1:8515/1‑137 (MQ=255)
                                                                                                     tCCCGTATCCAGCAGGTATAAAGCGCCAAATTGCGCTGATAGATCTGCCGCAGCCGCCGGATATGGCGATGATAGTGACCTTCCAGCACAAACGTTGCCGCCGCCATTTGCGTGGACGGCACATTAAAGTTGCTGAt                                                                                                                                                                                                                                                 <  2:8515/137‑1 (MQ=255)
                                                                                                           aTCCAGCAGGTATAAAGCGCCAAATTGCGCTGATAGATCTGCCGCATCCGCCGGATATGGCGATGATAGTGACCTTCCAGCACAAACGTTGCCGCCGCCATTTGCGTGGACGGCACATTAAAGTTGCTGATGGCGTATTTCATATGCATCAGTTTATCGTGATAACGCCCCGGTGCGACCCAACCCACGCGCAGGCCTGGGGCAATACTTTTACTGAACGAGATGCACAACAGCACTCGCCCGTCGATATc                                                                                                                         >  1:256204/1‑251 (MQ=255)
                                                                                                                             gcCAAATTGCGCTGATAGATCTGCCGCATCCGCCGGATATGGCGATGATAGTGACCTTCCAGCACAAACGTTGCCGCCGCCATTTGCGTGGACGGCACATTAAAGTTGCTGATGGCGTATTTCATATGCATCAGTTTATCGTGATAACGCCCCGGTGCGACCCAACCCACGCGCAGGCCTGGTGCAATACTTTTACTGAACGAGCTGCACAACAGCACTCGCCCGTCGATATCCCAGGAATGAATGGTCcg                                                                                                       >  1:255525/1‑251 (MQ=255)
                                                                                                                                        cTGATAGATCTGCCGCATTCGCCGGATATGGCGATGATAGTGACCTTCCAGCACAAACGTTGCCGCCGCAATTTGCGGGGACGGCACATTAAAGTTGCTGATGGCGTATTTCATATGCATCAGTTTATCGTGATAACGCACCGGTGCGACCCAACCCACGCGCAGGCCTGGGGCAATACTTTTCCTGAACGAGCTGCACAACATCACTCGCACGTCGATATCCATGGAATGAATGGTCCGCGGGCGAGGAt                                                                                            >  2:64404/1‑251 (MQ=255)
                                                                                                                                             aGATCTGCCGCATCCGCCGGATATGGCGATGATAGTGACCTTCCAGCACAAACGTTGCCGCCGCCATTTGCGTGGACGGCACATTAAAGTTGCTGAGGGCGTATTTCATATGCATCAGTTTATCGTGATAACGCCCCGGTGCGACCCAACCCACGCGCAGGCCTGGTGCAATACTTTTACTGAACGAGCTGCACAACAGCACTCGCCCGTCGATATCCCAGGAATGAATGGTCCGCGGGCGCGGATACTcc                                                                                       <  2:305564/251‑1 (MQ=255)
                                                                                                                                               aTCTGCCGCATCCGCCGGATATGGCGATGATAGTGACCTTCCAGCACAAACGTTGCCGCCGCCATTTGAGTGGACGGCACATTAAAGTTGCTGATAGCGTATTTCATCTGCCTCAGTTTATCGTGATAACGCACCGGTGCACCCCAACCCACGCGCAGGCCTGGTGCAATACTTTTCCTGAACGAGCTGAACAACAGCACTCGCCCGTCGATATCCCAGGAATGAATGGTCCGCGGGCGCGGATACTCCGt                                                                                     >  1:231348/1‑251 (MQ=255)
                                                                                                                                                   gCCGCATCCGCCGGATATGGCGATGATAGTGACCTTCCAGCACAAACGTTGCCGCCGCCATTTGCGTGGACGGCACATTAAAGTTGCTGATGGCGTATTTCATATGCATCAGTTTATCGTGATAACGCCCCGGTGCGACCCAACCCACGCGCAGGCCTGGTGCAATACTTTTACTGAACGAGCTGCACAACAGCACTCGCCCGTCGATATCCCAGGAATGAATGGTCCGCGGGCGCGGATACTCCGTCGcc                                                                                 >  1:246738/1‑251 (MQ=255)
                                                                                                                                                   gCCGCATCCGCCGGATATGGCGATGATAGTGACCTTCCAGCACAAACGTTGCCGCCGCCATTTGCGTGGACGGCACATTAAAGTTGCTGATGGCGTATTTCATATGCATCAGTTTATCGTGATAACGCCCCGGTGCGACCCAACCCACGCGCAGGCCTGGTGCAATACTTTTACTGAACGAGCTGCACAACAGCACTCGCCCGTCGATATCCCAGGAATGAATGGTCCGCGGGCGCGGATACTCCGTCGcc                                                                                 <  1:394159/251‑1 (MQ=255)
                                                                                                                                                             ccgGACATGGCGATGATAGTGACCTTCCAGCACAAACGCTGCCGCCGCCATTTGCGTGGACGGCACATTAAAGTTGCTGATCGCCTATTTCATCTGCATCAGTTTATCGTGATAACGCCCCGGTGCGACCCCACCCACCCGCAGGCCTGGTGCAATCCTTTTACTGAACGAGCTGCACAACAGCACTCGCCCGTCGATATCCCAGGAATGAATGGTTCGCGGGCGCGGGTCCTCTGTCGCCatgtctacc                                                                        >  2:371864/1‑242 (MQ=255)
                                                                                                                                                              cgGATTTGGCGAGGATAGTGACCTTCCAGCACAAACGTTGCCGCCGCCATTTGCGTGGACGGCACATTAAAGTTGCTGATGGCGTATTTCATATGCAGCAGTTTATCGTGATAACGCCCCGGTGCGACCCAACCCACGCGCAGGCCTGGTGCAATACTTTTACTGAACGAGCTGCACAACAGCACTCGCCCGTCGATATCCCAGGAATGAATGGTCCGCGGGCGCGGATACTCCGTCGCCAGTTCGCCATa                                                                      <  2:361192/251‑1 (MQ=255)
                                                                                                                                                                       cGATGATAGTGACCTTCCAGCACAAACGTTGCCGCCGCCATTTGCGTGGACGGCACATTAAAGTTGCTGATGGCGTACTTCATATGCATCAGTTTATCGTGATAACGCCCCGGTGCGACCCAACCCACGCGCAGGCCTGGTGCAATACTTTTACTGAACGAGCTGCACAACAGCACTCGCCCGTCGATATCCCAGGAATGAATGGTCCGCGGGCGCGGCTACTCCGTCGCCAGTTCGCCACAGCCatcatc                                                             >  2:110531/1‑251 (MQ=255)
                                                                                                                                                                              aGTGACCTTCCAGCATAAACGTTGCCGCCTCCAATTGCGTGGACGGCACATTAAAGTTGCTGATGGCGTATTTCATATGCATCAGTTTATCGTGATAACGCCCCGGTGCGACCCAACCCACGCGCAGGCCTTGTGCAATACGTTTACTGAACGATCTGCACAACAGCACTCGCCCGTCGATATCCCAGGAATGAATGGTCCGCGGGCGCGTATACTCCGTCGTTAGTTCGCCATAGACATCATCTTCaaaa                                                      <  2:401608/251‑1 (MQ=255)
                                                                                                                                                                               gggACCTTCCAGCACAAACGTTGCCTCCGCCATTTGCGTGGACGGCACATTAAAGTTGCTGATTGCGTATTTCATATGCATCAGTTTATCGTGATAACGCCCCGGTGCGACCCAACCCACGCGCAGGCCTGGTGCAATACTTTTACTGAACGAGCTGCACAACAGCACTCGCCCGTCGATATCCCAGGAATGAATTGTCCGCGGGCGCTGATACTCCGTCGCCAGTTCTCCATAGACATCATCGTCaaaaa                                                     <  2:167413/249‑1 (MQ=255)
                                                                                                                                                                                 ggccTTCCAGCTCAAACGGTGCCGCCGCCATTTGCGTGGACGGCACATTAAAGTTGCTGATGGCGTATTTCATATGCATCAGTTTATCGTGATAACGCCCCGGTGCGACCCAACCCACGCGCAGGCCTGGTGCAATACTTTGACTGAACGAGCTGCACAACAGCACTCGCCCGTCGATATCCCAGGAATGAATTGTCCGCGGGCGTGGATTCTCCGTCGCCAGTTCGCCAGAGTCATCATCTTCAAAAATc                                                   <  2:20881/249‑1 (MQ=255)
                                                                                                                                                                                  aCCTTCCAGCACAAACGTTGCCGCCGCCATTTGCGTGGACGGCACATTAAAGTTGCTGATGGCGTATTTCATATGCATCAGTTTATCGTGATAACGCCCCGGTGCGACCCAACCCACGCGCAGGCCTGGTGCAATACTTTTACTGAACGAGCTGCACAACAGCACTCGCCCGTCGATATCCCAGGAATGAATGGTCCGCGGGCGCGGATACTCCGTCGCCAGTTCGCCATAGACATCATCTTCAAAAATc                                                   >  2:96760/1‑250 (MQ=255)
                                                                                                                                                                                   ccTTCCAGCACAAACGTTGCCGCCGCCATTTGCGTGGCCGGCACATTAAAGCTGCTGATGGCGTATTTCATATGCATCAGTTTATCGTGATAACGCCCCGGTGCGACCCAACCCACGCGCAGGCCTGGTGCAATACTTTTACTGAACGAGCTGCACAACAGCACTCGCCCGTCGATATCCCAGGAATGACTGGTCCGCGGGCGCGGATACTCCGTCGCCAGTTCGCCATAGACATCATCTTCAAAAATcac                                                 >  2:62930/1‑251 (MQ=255)
                                                                                                                                                                                       ccAGCACAAACGTTGCCGCCGCCATTTGCGTGGACGGCACATTAAAGCTGCTGATGGCGTATTTCATCTGCATCAGATTATCGTGATACCGACCCGGTGCGACCCCACCCACGCGCAGGCCTGGTGCAATACTTTTACTGAACGAGCTGCACCACAGCACTCGCCCGTCGATATCCCAGGAATGAATGGTCAGCGGGCGCGGATACTCCGTCGCCAGTTCGACATAGACATCATCTTCAAAAATCACAat                                              >  1:24880/1‑250 (MQ=255)
                                                                                                                                                                                         aGCAAAAACGTTGCCGCCGCCATTTGCGTGGACGTCACATTAAAGTTGCTGATGGCGTATTTCATATGCATCAGTTTATCGTGATAACGCCCCGGTGCGACCCAACCCACGCGCAGGCCTGGTGCAATACTTTTACTGAACGAGCTGCACAACAGCACTCGCCCGTCGATATCCCAGGAATAAATGGTCCGCGGGCGCGGATACTCCGTCGCCAGTTCGCCATAGACATCATCTTCAAAAATCACAATATc                                           <  2:118855/251‑1 (MQ=255)
                                                                                                                                                                                           cacaAACGTTGCCGCCGCCATTTGCGTGGACGGCACATTAAAGTTGCTGATGGCGTATTTCATATGCATCAGTTTATCGTGATAACGCCCCGGTGCGACCCAACCCACGCGCAGGCCTGGTGCAATACTTTTACTGAACGAGCTGCACAACAGCACGCGCCCGTCGATATCCCAGGAATGAATGGTCCGCGGGCGCGGATACTCCGTCGCCAGTTCGCCATAGACATCATCTTCAAAAATCACAATATCAt                                         <  1:264751/251‑1 (MQ=255)
                                                                                                                                                                                              aaaCGTTGCCGCCGCCATTTGCGTGGACGGCACATTAAAGTTGCTGATGGCGTATTTCATATGCATCAGTTTATCGTGATAACGCCCCGGTGCGACCCAACCCACGCGCAGGCCTGGTGCAATACTTTTACTGAACGAGCTGCACAACAGCACTCGCCCGTCGATATCCCAGGAATGAATGGTCCGCGGGCGCGGATACTCCGTCGCCAGTTCGCCATAGACATCATCTTCAAAAATCACAATATCATGAc                                      >  1:162195/1‑251 (MQ=255)
                                                                                                                                                                                                    tGCCGCCGCCATTTGCGTGGACGGCACATTAAAGTTGCTGATGGCGTATTTCATATGCATCAGTTTATCGTGATAACGCCCCGGTGCGACCCAACCCACGCGCAGGCCTGGTGCAATACTTTTACTGAACGAGCTGCACAACAGCACTCGCCCGTCGATATCCCAGGAATGAATGGTCCGCGGGCGCGGATACTCCGTCGCCAGTTCGCCATAGACATCATCTTCAAACATCACAATATCATGACGCTg                                  >  1:279994/1‑249 (MQ=255)
                                                                                                                                                                                                             cATTTGCGTGGACGGCACATTAAAGTTGCTGATGGCGTATTTCATATGCATCAGTTTATCGTGATAACGCCCCGGTGCGACCCAACCCACGCGCAGGCCTGGTGCAATACTTTTACTGAACGAGCTGCACAACAGCACTCGCCCGTCGATATCCCAGGAATGAATGGTCCGCGGGCGCGGATACTCCGTCGCCAGTTCGCCATAGACATCATCTTCAAAAATCACAATATCATGACGCTGAGCgagag                          >  1:342893/1‑248 (MQ=255)
                                                                                                                                                                                                                 tgcgtcgtcggcggaTTAAAGTTGCTGATGGCGTAGTTCATATGCATCAGTTTATCGTGATAACGCCCCGGTGCGACCCAACCCACGCGCAGGCCTGGTGCAATACTTTTACTGAACGAGCTGCACAACAGCACTCGCCCGTCGATATCCCAGGAATGAATGGTCCGCGGGCGCGGATGCGCCGTGGCCCGATCGCCATAGACATCGTCTTCAAAAATCAGAATGTCATGACGGGGAGCGAGAGAGAGaa                    <  2:312162/236‑1 (MQ=255)
                                                                                                                                                                                                                               aTTAAAGTTGCTGATGGCGTATTTCATATGCATCAGTTTATCGTGATAACGCCCCGGTGCGACCCAACCCACGCGCAGGCCTGGTGCAATACTTTTACTGAACGAGCTGCACAACAGCACTCGCCCGTCGATATCCCAGGAATGAATGGTCCGCGGGCGCGGATACTCCGTCGCCAGTTCGCCATAGACATCATCTTCAAAAATCACAATATCATGACGCTGAGCGAGAGAGAGAACGGCCCGTTTgcgcg     >  1:248750/1‑251 (MQ=255)
                                                                                                                                                                                                                                  aaaGTTGCTGATGGCGTATTTCATATGCATCAGTTTATCGTGATAACGCCCCGGTGCGACCCAACCCACGCGCAGGCCTGGTGCAATACTTTTACTGAACGAGCTGCACAACAGCACTCGCCCGTCGATATCCCAGGAATGAATGGTCCGCGGGCGCGGATACTCCGTCGCCAGTTCGCCATAGACATCATCTTCAAAAATCACAATATCATGACGCTGAGCGAGAGAGAGAACGGCCCGTTTGCGCGCGt  <  1:62930/251‑1 (MQ=255)
                                                                                                                                                                                                                                  aaaGTTGCTGATGGCGTATTTCATATGCATCAGTTTATCGTGATAACGCCCCGGTGCGACCCAACCCACGCGCAGGCCTGGTGCAATACTTTTACTGAACGAGCTGCACAACAGCACTCGCCCGTCGATATCCCAGGAATGAATGGTCCGCGGGCGCGGATACTCCGTCGCCAGTTCGCCATAGACATCATCTTCAAAAATCACAATATCATGACGCTGAGCGAGAGAGAGAAAGGACCGTTTGCGCGCGt  >  1:294946/1‑251 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
GCTGCCGCGCGACGCAGACCATATCGACCTGTTCAGGCAATTCTATCCACAGTAAAAATCCGCCTTTCGGGCGCGTAATACAGATTTCGCAGGGAAAATATTCCCGTATCCAGCAGGTATAAAGCGCCAAATTGCGCTGATAGATCTGCCGCATCCGCCGGATATGGCGATGATAGTGACCTTCCAGCACAAACGTTGCCGCCGCCATTTGCGTGGACGGCACATTAAAGCTGCTGATGGCGTATTTCATATGCATCAGTTTATCGTGATAACGCCCCGGTGCGACCCAACCCACGCGCAGGCCTGGTGCAATACTTTTACTGAACGAGCTGCACAACAGCACTCGCCCGTCGATATCCCAGGAATGAATGGTCCGCGGGCGCGGATACTCCGTCGCCAGTTCGCCATAGACATCATCTTCAAAAATCACAATATCATGACGCTGAGCGAGAGAGAGAACGGCCCGTTTGCGCGCGT  >  NC_000913/4570320‑4570796

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