Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,660,200 C→G 6.8% A711G (GCC→GGC)  ynfE → putative selenate reductase, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,660,2000CG 6.8% 164.9 / 4.9 59A711G (GCC→GGC) ynfEputative selenate reductase, periplasmic
Reads supporting (aligned to +/- strand):  ref base C (31/24);  new base G (2/2);  total (33/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

GAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATG  >  NC_000913/1660052‑1660348
                                                                                                                                                    |                                                                                                                                                    
gAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGctgc                                                                                                                                                      <  2:360192/149‑1 (MQ=255)
  cgTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGATGGAACTCCCCTGAACGTAGAACCTTCCCACTGAAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGctgcct                                                                                                                                                    <  1:867928/148‑1 (MQ=255)
  aGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGctgcct                                                                                                                                                    <  1:334037/149‑1 (MQ=255)
      aTAAGCCAAATGCCAGTTTATGCCTCTTCCTTTGAAGGCTGGATATCACCTGATCGTAGAACCTTCCCACTGCAACTGTGTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGACg                                                                                                                                                <  2:807873/149‑1 (MQ=255)
          gCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGATCCTTCCCACTGCAACTTGTTGGTTTCCATTACAAATCCCGTACTCACTCGACCAACGGCAATATTGATCTCCAGAAGGCTGCCTGCCGTCag                                                                                                                                            >  1:515326/1‑149 (MQ=255)
          ccccATAGCCAGTTTAAGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTACCATTACAACACCCGACCTCACTCGACCTACGGCAATATTGATCTCCTGACGGCAGGCAGCCTTCag                                                                                                                                            <  1:819378/148‑1 (MQ=255)
           cccATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCagg                                                                                                                                           >  1:39515/1‑149 (MQ=255)
            ccATTGCCAGTTTATTCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACTTAGAACCTTACCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTTCGGCAATATTGATCTCCTGACGGCAGGCAGCCTTCagga                                                                                                                                          <  2:102958/149‑1 (MQ=255)
             ccttGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGACCCTTCCCACTGCAACTGTTTGGTTTCCATTACAACTCCCGTACTCGCTCTCCCTCCGGCAAGATTGCTCTCCTGACGGCTGCCTGCCGTCcggag                                                                                                                                         >  2:213836/3‑149 (MQ=255)
                 gCCAGTTTATGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGgtgt                                                                                                                                     <  2:150306/149‑1 (MQ=255)
                  ccAGTTTATGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGACGGCAGGCAGCCTTCaggagatcaa                                                                                                                                    >  1:288404/1‑144 (MQ=255)
                   cAGTTTATGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACATACCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTgg                                                                                                                                   >  1:270605/1‑149 (MQ=255)
                          tGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAAGCTTCCCTCTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCTCTCGACCTACTGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCATcc                                                                                                                            <  2:1143740/149‑1 (MQ=255)
                          tGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAAcc                                                                                                                            <  2:1159987/149‑1 (MQ=255)
                          cgCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAAcc                                                                                                                            <  2:1143739/148‑1 (MQ=255)
                           gCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAAccc                                                                                                                           <  1:501092/149‑1 (MQ=255)
                              tCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAGCCCtat                                                                                                                        >  1:32075/1‑149 (MQ=255)
                              tCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCtat                                                                                                                        <  2:418485/149‑1 (MQ=255)
                              tCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCtat                                                                                                                        >  1:274864/1‑149 (MQ=255)
                              tCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCtat                                                                                                                        >  1:1109831/1‑149 (MQ=255)
                                      tGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATgcgc                                                                                                                <  1:37166/149‑1 (MQ=255)
                                      tGAAGGCTGGAACTCCCCTGAACGTAGAACCATCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATgcgc                                                                                                                <  2:873556/149‑1 (MQ=255)
                                              ggAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGt                                                                                                        >  1:518393/1‑149 (MQ=255)
                                                   tCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATATCCTGAAGGCTGCCTGCCGTCAGTAGGTGTGGATCAACACTATAGATGCGCAGAAACGTGGGAt                                                                                                   >  1:510064/1‑149 (MQ=255)
                                                    gcccTGAACTTAGAACCATCGATCTTCTACTGTTTGCTGTCAAAAACACTTCCCGTACTCACTCGACCTACGGCACTATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCATATAGATGCGCAGAAACGTGGGAtt                                                                                                  <  2:1096942/148‑1 (MQ=255)
                                                    ccccTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGAtt                                                                                                  <  2:88854/149‑1 (MQ=255)
                                                     ccttGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTAAAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCTGAAACGTGGGTTTg                                                                                                 <  2:52714/146‑1 (MQ=255)
                                                     cccTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATATCCTGAAGGCTGCCTGCCGGCAGGAGGTGTGGAGCAACCCTATAGATGCGCAGAAACGTGGGATTg                                                                                                 >  1:378415/1‑149 (MQ=255)
                                                       cTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGcc                                                                                               >  1:734799/1‑149 (MQ=255)
                                                       cTGAACCTAGAACCTTCCCACTGCAACTGTATGGTTTCCATTACAAATAACGTACTCACTCGCCCTACGGCAATATTGATCTCCTGAAGGCTGCCTCCCGTCAGGAGGTGTCGATCACCCCTATAGATGCGCAGAAACGTGGGATTTcc                                                                                               >  2:1117998/1‑149 (MQ=255)
                                                         gAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTCCTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCaa                                                                                             <  2:553610/149‑1 (MQ=255)
                                                            cGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCAGTACTCACTCGACCTACGGCAATATTGATCTGCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAAAGTGGGATTGCCAACgg                                                                                          <  2:1023438/149‑1 (MQ=255)
                                                              gaGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTCCTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGCTGCGCAGAAACGTGGGATTGCCAACGGCg                                                                                        <  2:337365/148‑1 (MQ=255)
                                                                   ccTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAAGCCTATAGCTGCGCAGAAACGTGGTAATGCCAACGGCGATATg                                                                                   >  2:1045758/1‑149 (MQ=255)
                                                                   ccTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACCTGGGATTGCCAACGGCGATATg                                                                                   >  1:100593/1‑149 (MQ=255)
                                                                       cccACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGc                                                                               >  2:820212/1‑149 (MQ=255)
                                                                          aCTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCggg                                                                            >  1:335560/1‑149 (MQ=255)
                                                                              cAACTGTTTGGTTTCCATTACAAATCCCGTCCGCCATCGACCTCCGGCAATATTGATCTCCTGCAGGCTGCCTGCCGTCAGGGGGGGTGGATCCCCCCTATAGATGCGCCGCAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGtt                                                                        >  1:893555/1‑149 (MQ=255)
                                                                              cAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGtt                                                                        <  2:270605/149‑1 (MQ=255)
                                                                                  tGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAAc                                                                    >  1:304651/1‑149 (MQ=255)
                                                                                               ttACAAATCCCGTACTCACTCCACCTACGGCAATATTGATCTCCTGAAGGCTGCCTCCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGCCGATATTGTGCGGGTCTTTAACCATCGTGGCGAAg                                                       >  1:121434/1‑149 (MQ=255)
                                                                                                     aTCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGACGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGc                                                 <  1:131717/149‑1 (MQ=255)
                                                                                                     aTCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGc                                                 <  1:217101/149‑1 (MQ=255)
                                                                                                      tCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACACTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGGGCGGGTGTTTACCCATCGTGGCGAAGTTCGGCt                                                >  1:796084/1‑149 (MQ=255)
                                                                                                        ccGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCCACCCTATAGATGCGCAGAACCGTGGGAATGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTAc                                              >  2:687637/1‑149 (MQ=255)
                                                                                                          gaagtcactcGTCCTACGGCAATTTTGATCTCCTGATGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGCTTTTGTCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCa                                            <  2:438861/145‑1 (MQ=255)
                                                                                                          gTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCa                                            >  1:152291/1‑149 (MQ=255)
                                                                                                           tACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAg                                           <  1:870504/149‑1 (MQ=255)
                                                                                                                actcGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCaaa                                      >  1:1024060/1‑149 (MQ=255)
                                                                                                                 ctcGACCTACGGCAATAATGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAg                                     <  2:39515/149‑1 (MQ=255)
                                                                                                                    gACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTaa                                  <  1:820212/149‑1 (MQ=255)
                                                                                                                     aCCTACCGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTCCCAACGGCGATATGGTGCGGGTGTTTAAACATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAac                                 <  2:510064/149‑1 (MQ=255)
                                                                                                                        tACGGCAATATTGATCTCCTGAAGGCTGCATGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGACAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAAcacc                              >  2:71024/1‑149 (MQ=255)
                                                                                                                              aaTATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTAt                        >  2:137870/1‑149 (MQ=255)
                                                                                                                               atatTGATCTCCTGAAGTCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGGGTTTAACCATCGTTGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTAtt                       <  1:298173/149‑1 (MQ=255)
                                                                                                                                   tGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTcc                   >  2:1133794/1‑149 (MQ=255)
                                                                                                                                       ctcCTGACGGCAGGCAGCCTTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGAATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTgg               >  1:282754/1‑149 (MQ=255)
                                                                                                                                       ctcCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAAATGAGCATAAACGTGGGATTGCCAACGTCCATATGGTGCGGGTGTTTAACAATCGAGGCCAAGTTCGGCTACCAGACAAATGAACACCACGTATTCTCCCtag               >  2:495903/1‑147 (MQ=255)
                                                                                                                                        tcCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGa              >  2:15319/1‑149 (MQ=255)
                                                                                                                                          cTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGt            >  1:416876/1‑149 (MQ=255)
                                                                                                                                          cTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGt            >  1:1105143/1‑149 (MQ=255)
                                                                                                                                            gAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGCAGttt          >  2:231090/1‑148 (MQ=255)
                                                                                                                                                   gcctgccGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTAt   >  1:837527/1‑149 (MQ=255)
                                                                                                                                                    cctgccGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCCCGTATTCTCCCTGGAGTTAGCGCTATg  >  2:855857/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GAAGTGATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATG  >  NC_000913/1660052‑1660348

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