Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,276,422 C→T P318S (CCG→TCG)  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,4220CT98.0% 149.6 / NA 51P318S (CCG→TCG) pdeCc‑di‑GMP phosphodiesterase PdeC
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (27/22);  minor base G (0/1);  total (28/23)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTT  >  NC_000913/4276179‑4276667
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     
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)
                                               ggggATTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACt                                                                                                                                                                                                  <  2:170185/250‑1 (MQ=255)
                                               ggggAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTg                                                                                                                                                                                                 <  1:79725/251‑1 (MQ=255)
                                               ggggAGTCTTGTTCTGCTATTACACTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGACCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGCAGGGGTTGATAGAACAGATAACTg                                                                                                                                                                                                 >  2:113626/1‑251 (MQ=255)
                                                   aGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGAtt                                                                                                                                                                                              >  1:404623/1‑250 (MQ=255)
                                                          ttCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTac                                                                                                                                                                                       <  2:137782/250‑2 (MQ=255)
                                                           tCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTAtt                                                                                                                                                                                      <  1:70155/250‑1 (MQ=255)
                                                              gCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTGTGGCTTTATTACCAGCCAAGAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTGATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGa                                                                                                                                                                                    <  1:500886/249‑1 (MQ=255)
                                                                   tACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTACTGATAATGTc                                                                                                                                                                             >  1:381587/1‑251 (MQ=255)
                                                                    aCTCTAGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATACTCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGGGTTTATTCCGCTGACAGCAAAGGAGGGGATGATAAAACAGATAACTGATTGTGTTATTGATAATGTCt                                                                                                                                                                            >  1:350654/1‑251 (MQ=255)
                                                                                       cAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTACGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGAtt                                                                                                                                                          >  1:120290/1‑250 (MQ=255)
                                                                                         aaaCTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACTTCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTAcc                                                                                                                                                       >  1:188592/1‑251 (MQ=255)
                                                                                              aTTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCa                                                                                                                                                  <  1:93354/251‑1 (MQ=255)
                                                                                              aTTTATCTCCCAAACGTAAATTGCAACGCCCCCTCGCAAAACATCAACTTTGTCTTTATTATCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTTCCTgg                                                                                                                                                    <  2:377059/249‑1 (MQ=255)
                                                                                                       ccAAACGGAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCATGGTGAGCAGGGGCAAATAATGACTTCGGCAGAGTTTATTCCGCTGGCAGACAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGAGATTACCTGGAAACACATGc                                                                                                                                          >  2:408753/1‑250 (MQ=255)
                                                                                                           aCGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGAAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAggac                                                                                                                                     >  2:301492/1‑248 (MQ=255)
                                                                                                                aatgtaaaccatcacccaaaaaaCATCTATTGTGTCTTTATTACCAGCCAATAATCGAGATCAAAACAGAAAAATGTACCGGCGCTGATGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGGCTTCCGCGATCTGGGCGATTACCGGGCAACACATGCAGATCGCttt                                                                                                                                <  2:350654/233‑3 (MQ=255)
                                                                                                                               ccGAAAAACATCAACTTTGTCTTTAATACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTg                                                                                                                 <  2:420208/250‑1 (MQ=255)
                                                                                                                                cGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCt                                                                                                                  <  2:133187/249‑1 (MQ=255)
                                                                                                                                 aaaaaaCATCTACTTTGTCTTTATTACCAGCCAATAATCGAAATCAAAACAGAAAAATGGATCGGCGCTGAAGCGTTGTGACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCATAAAAGGAGGGGATGATAGAACAGATAACTGATGATGTTATTGATAATGTCTTCCGCTATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTTCCCTGTc                                                                                                               <  2:265401/250‑1 (MQ=255)
                                                                                                                                             cTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGCAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACCCATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGAttt                                                                                                   >  1:272939/1‑251 (MQ=255)
                                                                                                                                             cTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGAttt                                                                                                   >  2:464324/1‑251 (MQ=255)
                                                                                                                                                       attaACAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGCGTTTTTTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTCTGTTTCTATTAACCTGTCGGCCTCCGATTTTAATACGTc                                                                                          >  2:132514/1‑250 (MQ=255)
                                                                                                                                                        ttaacAGCCAATAATCGCTATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGGGCGAAGGGGCAACTAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGTTATGTTTCTTTTAACCTGTCGGCCTCCGAtttttgtacgtgac                                                                                        >  2:458922/5‑241 (MQ=255)
                                                                                                                                                             aGCCAATAATCGATATCAAAACAGAAAAATGTATTGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGAAGGAGCAAATAATGAATTCGGCAGAGTTTATTACGCTGGCAGAAAAGGAGGAGCTGATAGAACAGATAACTGATTATGTTATTGATAATGACTTCCACGATATGGGCGATTACCTGGCAACACATGAGGATCGCTATGTATCTATTAACCTGTCGGCTTACGATTTTCATaagatccggttg                                                                                   >  2:276573/1‑240 (MQ=255)
                                                                                                                                                             aGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACa                                                                                                                                                                                                         >  2:450636/1‑133 (MQ=255)
                                                                                                                                                             aGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACa                                                                                                                                                                                                         <  1:450636/133‑1 (MQ=255)
                                                                                                                                                               ccAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGAt                                                                                 >  2:341698/1‑251 (MQ=255)
                                                                                                                                                                     aacGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGcgcg                                                                           <  1:306894/249‑1 (MQ=255)
                                                                                                                                                                          tatCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGCATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTACTAGCGCGaatc                                                                       >  2:485709/1‑250 (MQ=255)
                                                                                                                                                                           atCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGaatca                                                                      >  1:41326/1‑250 (MQ=255)
                                                                                                                                                                             cAAAATAGAAAAATGTAGCGGCTCGAAAGCGTTGGTACGTTGGCTTGGTGAGCAGGGGCAAATAATGAATGCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTAGTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCCGCCTCCTATTTTCATACGGCACGGTTGATAGCGCGaatcaat                                                                    <  2:183133/250‑1 (MQ=255)
                                                                                                                                                                                     aaaaaTGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACa                                                           <  1:357624/251‑1 (MQ=255)
                                                                                                                                                                                       aaaTGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATGTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACaga                                                         <  2:69595/251‑1 (MQ=255)
                                                                                                                                                                                         aTGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGCTTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGc                                                       >  1:239147/1‑251 (MQ=255)
                                                                                                                                                                                          tGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGc                                                       >  1:172562/1‑250 (MQ=255)
                                                                                                                                                                                           gTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCa                                                      >  1:139082/1‑250 (MQ=255)
                                                                                                                                                                                                  cgcTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCgg                                              >  2:357831/1‑251 (MQ=255)
                                                                                                                                                                                                              gttACGTTGGCCGGGTGAGCAGGGGCAAATAAGGAATTCGGCAGATTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTTATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCgcag                                   <  2:276690/250‑1 (MQ=255)
                                                                                                                                                                                                                      ggCCTGGTGAGCAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCaaattaa                           <  1:405063/250‑1 (MQ=255)
                                                                                                                                                                                                                                 cAGGGGCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGAc               <  2:47092/251‑1 (MQ=255)
                                                                                                                                                                                                                                   ggggCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTg             >  1:58127/1‑251 (MQ=255)
                                                                                                                                                                                                                                     ggCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGaatca                                                                      <  2:240006/192‑1 (MQ=255)
                                                                                                                                                                                                                                     ggCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGaatca                                                                      >  1:240006/1‑192 (MQ=255)
                                                                                                                                                                                                                                     ggCAAATAATGAATTCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAg           <  2:139082/251‑1 (MQ=255)
                                                                                                                                                                                                                                               gAATTCGGCAGAGGTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAAGGAAATTTGAAGTGACTGAGCATGCAttt  <  2:33514/250‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     
ACCCATTTTTATAATCATCTTATTTTTGCGTTGCCCGCCGGTATTTTGGGGAGTCTTGTTCTGCTATTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGATCGCTATGTTTCTATTAACCTGTCGGCCTCCGATTTTCATACGTCACGGTTGATAGCGCGAATCAATCAGAAAACAGAGCAATACGCGGTGCGTCCGCAGCAAATTAAATTTGAAGTGACTGAGCATGCATTT  >  NC_000913/4276179‑4276667

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: