Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,276,212 C→T P248S (CCC→TCC)  pdeC → c‑di‑GMP phosphodiesterase PdeC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,276,2120CT92.3% 109.2 / ‑1.1 39P248S (CCC→TCC) pdeCc‑di‑GMP phosphodiesterase PdeC
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (23/13);  total (24/15)
Fisher's exact test for biased strand distribution p-value = 5.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

ACCCGACATTGCAATGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGG  >  NC_000913/4275976‑4276455
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                   
aCCCGACATTGCAATGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGAAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGCCAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATATCCCATTTTTATAATCATCTTATTTTTGCGTTGTCCGCAGGTATTTTg                                                                                                                                                                                                                                       >  2:136640/1‑251 (MQ=255)
                                        aaCCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTTCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGTCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGAc                                                                                                                                                                                               <  2:71149/251‑1 (MQ=255)
                                         aCCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGTCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGAc                                                                                                                                                                                               >  1:276690/1‑250 (MQ=255)
                                             ttttCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGTCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACaaa                                                                                                                                                                                            <  1:233026/249‑1 (MQ=255)
                                             ttttCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGTCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAAc                                                                                                                                                                                          >  1:420208/1‑251 (MQ=255)
                                                     ttAAGCAAAGAGGCCTCGGCAGCAACGCTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGCCAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGTCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTAGACAAAACTATTTAtc                                                                                                                                                                                  >  1:183133/1‑251 (MQ=255)
                                                        aGCAAAGAGGCCTCGGAGACAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGTCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTcc                                                                                                                                                                               <  2:29795/251‑1 (MQ=255)
                                                         gCAAAGAGGCCTAGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGAGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGTCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTATGGCTACGTATTCGACAAAACTATTTAtctcac                                                                                                                                                                              >  1:171221/1‑249 (MQ=255)
                                                          cAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGTCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCa                                                                                                                                                                             <  1:10669/251‑1 (MQ=255)
                                                                      ggCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGCTTTAACCGTACAAAGCAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGTTCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGAAAATTGc                                                                                                                                                                 >  2:82555/1‑251 (MQ=255)
                                                                         agcaACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAACCGCCCAATGGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGTCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAAc                                                                                                                                                              <  1:385885/251‑1 (MQ=255)
                                                                          agaaaCGTTTTCTCCGCTGATTAATTTTAAGGATTTAACCGATCAAAAAAATGGTTTTTAATAGGCGACAGTTTATGCGACAAAACGCCAAATGGCAGCCATTGTTGCGACTTCATATCAACGTCTTAGACCCTATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAAc                                                                                                                                                              <  1:223330/248‑1 (MQ=255)
                                                                                gTTTTCGCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGTCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGccc                                                                                                                                                        <  1:289299/251‑1 (MQ=255)
                                                                                  tttCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGTCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCt                                                                                                                                                       <  1:181575/250‑1 (MQ=255)
                                                                                    tctcCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGTCCGCCGGTATTTTGGGGAGTCTTGCTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCg                                                                                                                                                     <  2:89061/250‑1 (MQ=255)
                                                                                         gCTGATTCATTTGAAGGCTTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGTCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTAt                                                                                                                                                                                   >  2:455319/1‑215 (MQ=255)
                                                                                         gCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGTCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTAt                                                                                                                                                                                   <  1:455319/215‑1 (MQ=255)
                                                                                                     gAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGTCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTc                                                                                                                                   >  1:47092/1‑251 (MQ=255)
                                                                                                     gAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGTCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGACCTAGAAAAACATCAACTTTGTc                                                                                                                                   >  1:298958/1‑251 (MQ=255)
                                                                                                         attttAACCTTACAAAGAAAGGGTTATTTATATGCGACATTTTATCCCAAAAAACGCCCATTTGCAGCCATTGTTGCGCCTTCATATCAACGTCTTATAACCCATTTTTATAATCTTTTTTTTTTGGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTtta                                                                                                                               <  1:45007/249‑1 (MQ=255)
                                                                                                                                      atatGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGTCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAAc                                                                                                    >  1:102559/1‑249 (MQ=255)
                                                                                                                                        atGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGTCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGa                                                                                                 >  1:52408/1‑250 (MQ=255)
                                                                                                                                                 tttATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTATTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGTCCGCCGGTCTTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCACTAATCGATATCAAAACAGAAAAATGTATc                                                                                       >  1:33514/1‑251 (MQ=255)
                                                                                                                                                        gACAAAACGCGCAATTGCAGCCATTGTTACGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGTCCGCCGGTATTTTGGGGATTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGAAACGCGCCCTAGAAAAACATAAACTTTGTCTTTATTACCTGCCAATAATCGATATCAAAACAGAAAAATGAATCGGCGCTg                                                                                >  2:506511/1‑251 (MQ=255)
                                                                                                                                                          cAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGTCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAACCATCAACTTTGTCTTTATTACAAGCCAATAATCGTTATCAAAACAGAAAAATGAATCGGCGCtaa                                                                               >  2:405063/1‑248 (MQ=255)
                                                                                                                                                             aactccaaaTTGCAGCCATTGTTGCGAGTTCATATCAACGTCTTATAACCCTTTTTTATAATCATCTTATTTTTGCGTTGTCCGCCGGTATTTTGGGGAGTCTTGTTTAGCTATTACTCTGGCTACGTATTCGACAAATCTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCg                                                                           <  2:345512/244‑1 (MQ=255)
                                                                                                                                                               cGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGTCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCgtt                                                                         >  1:226640/1‑251 (MQ=255)
                                                                                                                                                                          aGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGTCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGc                                                              >  1:168878/1‑251 (MQ=255)
                                                                                                                                                                             cATTGTTGCGACTTCATATCAACGTCTTATACCCCATTTTTATATTCATCTTATTTTTGCGTTGTCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTg                                                           <  2:102559/251‑1 (MQ=255)
                                                                                                                                                                             cATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGTCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTg                                                           >  1:130630/1‑251 (MQ=255)
                                                                                                                                                                                                 aaCGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGTCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGTCCTGGTGAGCGGGGGCAaataat                                        >  2:227441/1‑250 (MQ=255)
                                                                                                                                                                                                  aCGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGTCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCaataat                                                                                                                <  2:150767/177‑1 (MQ=255)
                                                                                                                                                                                                  aCGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGTCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCaataat                                                                                                                >  1:150767/1‑177 (MQ=255)
                                                                                                                                                                                                            aCCCATTTTTATAATCATCTTATTTTTGCGTTGTCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGc                               >  1:28194/1‑248 (MQ=255)
                                                                                                                                                                                                                       taATCATCTTATTTTTGCGTTGTCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGc                 >  1:396135/1‑251 (MQ=255)
                                                                                                                                                                                                                        aaTCATCGTATTTTTGCGTTGTCCGCCGGTATTTTGGGGAGGCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACGATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGGCTTTATTCCCAGCCAATAAGCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCTTTGGTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCt                <  1:34449/251‑1 (MQ=255)
                                                                                                                                                                                                                               ttATTTTTGCGTTGTCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCACTAATCGATATCAAAACAGAAATATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGa          >  2:2555/1‑250 (MQ=255)
                                                                                                                                                                                                                                    tttGCGTTGTCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAagg     >  1:129257/1‑250 (MQ=255)
                                                                                                                                                                                                                                      tGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAaggagg  >  1:73360/1‑251 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                   
ACCCGACATTGCAATGGGGTGTGTATGATACGGTGACGAAAACCTTTTTCTCGTTAAGCAAAGAGGCCTCGGCAGCAACGTTTTCTCCGCTGATTCATTTGAAGGATTTAACCGTACAAAGAAATGGCTATTTATATGCGACAGTTTATTCGACAAAACGCCCAATTGCAGCCATTGTTGCGACTTCATATCAACGTCTTATAACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGG  >  NC_000913/4275976‑4276455

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