Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,423,383 A→G 5.7% intergenic (+189/+40) yhdZ → / ← rrfF putative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,423,3830AG 5.7% 658.8 / 17.3 207intergenic (+189/+40)yhdZ/rrfFputative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon
Reads supporting (aligned to +/- strand):  ref base A (107/85);  new base G (9/3);  total (119/88)
Fisher's exact test for biased strand distribution p-value = 2.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 218 total aligned reads displayed.

GAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAG  >  NC_000913/3423114‑3423654
                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                               
gAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCCACCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACCCGGGAGGGGGTTCACCGACAAACAACAGATAAAACAAAAg                                                                                                                                                                                                                                                                             >  2:432317/1‑274 (MQ=255)
                atgaATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTTCTTCTCGGGGGAGTGTTCACCGACAAACAACAGCTAAAACAAAAGGCCCAGGCTTCCGACTGAGcc                                                                                                                                                                                                                                                        >  1:87848/1‑279 (MQ=255)
                           tgcaagctgaaaGCAGGGAGTAGGAGGGCAGTTGTATCGGAGGTGATCCAGTAATTGAAGGTGAATTCGAAAAGCAAAAAGGCCATCCTTGCGGATGGCCTTTGGCTTGATTTGATGTCTGTCCGTTTGTGGGGGGGGACATGCCCGCCGCCCTCCGGGCCGTTGGTTCGCAACGTTCAAATCCGCTCCCGGCTGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGtt                                                                                                                                                                                                                                                 <  4:290094/266‑1 (MQ=17)
                                       tgggggggtagaaGGGCATTTTAATAGAAGGTNANTCATTAAATGAATGTTAGTTGGAAAAGCAAAAAGGCCATCCGTTCNGGTTGCCTTGCGCTTGATTTGATGTATGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGNTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCTGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGc                                                                                                                                                                                                                                 <  3:87848/267‑1 (MQ=17)
                                           agCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAAACTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACCACAGATAAAACAAAAGGCCCAGTCTTCCCACTGAGCCTTTTGTTTTTTTTGATGTCttg                                                                                                                                                                                                                                  >  2:161485/1‑272 (MQ=255)
                                                                                                                    ttCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCAt                                                                                                                                                        <  2:286508/275‑1 (MQ=255)
                                                                                                                            gCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGCCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACATTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCAACGACAAACAACAGATAAAACaaaa                                                                                                                                                                                                                                                                              <  6:1412455/149‑1 (MQ=255)
                                                                                                                              cTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAgg                                                                                                                                                                                                                                                                            >  6:1702498/1‑149 (MQ=255)
                                                                                                                                ttCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGcc                                                                                                                                                                                                                                                                          <  8:1585775/149‑1 (MQ=255)
                                                                                                                                     ttgatttgatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTc                                                                                                                                                                                                                                                                     >  7:1161953/1‑149 (MQ=255)
                                                                                                                                     agaatttgcctggcggccgtagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCTCAGTc                                                                                                                                                                                                                                                                     >  8:193457‑M1/130‑149 (MQ=255)
                                                                                                                                       gatttgatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCtt                                                                                                                                                                                                                                                                   >  5:1020140/1‑149 (MQ=255)
                                                                                                                                          ttgatGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCg                                                                                                                                                                                                                                                                <  8:29396/149‑1 (MQ=40)
                                                                                                                                              tGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTg                                                                                                                                                                                                                                                            >  6:1258345/1‑149 (MQ=39)
                                                                                                                                               gTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGa                                                                                                                                                                                                                                                           <  7:569653/149‑1 (MQ=39)
                                                                                                                                                 cTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGc                                                                                                                                                                                                                                                         >  6:1606582/1‑149 (MQ=38)
                                                                                                                                                  tGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCTCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACCAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTTATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGCCGCTACGGCGTTTCCCTTCCGAGt                                                                                                                                                               >  1:109090/1‑238 (MQ=255)
                                                                                                                                                  tGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGcc                                                                                                                                                                                                                                                        >  5:1588654/1‑149 (MQ=38)
                                                                                                                                                   ggCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCt                                                                                                                                                                                                                                                       >  8:1041705/1‑149 (MQ=37)
                                                                                                                                                   ggCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCCCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCCGTTCCCTACTACCGCATGGGGGGACCCCACACTACCACCCGCCCTCCCGGCTTTCACTTCCTGGTTCGGCGTGGGGTCAGGGGGGGCCCCcgctacacttccc                                                                                                                         >  4:537609/1‑265 (MQ=11)
                                                                                                                                                    gcatgtaaTGGCGGCCGCCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTgccgc                                                                                                                        <  2:404412/268‑1 (MQ=255)
                                                                                                                                                    gcagtagcgcgctggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCtt                                                                                                                                                                                                                                                      >  5:1320199‑M1/115‑149 (MQ=255)
                                                                                                                                                    gCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCtt                                                                                                                                                                                                                                                      >  5:970847/1‑149 (MQ=37)
                                                                                                                                                      aGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCtttt                                                                                                                                                                                                                                                    <  6:1588654/149‑1 (MQ=40)
                                                                                                                                                       gtagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTTCGACTGAGcctttcg                                                                                                                                                                                                                                                   >  7:1222728‑M1/112‑147 (MQ=255)
                                                                                                                                                         ttATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGtt                                                                                                                                                                                                                                                 <  8:368350/149‑1 (MQ=40)
                                                                                                                                                         agcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctcctcatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTCCGACTGAGCCTTTCGtt                                                                                                                                                                                                                                                 >  5:1670913‑M1/110‑149 (MQ=255)
                                                                                                                                                            gcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTTCGACTGAgcctttcgt                                                                                                                                                                                                                                                  <  8:487328‑M1/39‑4 (MQ=255)
                                                                                                                                                             ggcgggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTAt                                                                                                                                                                                                                                             >  6:1294440/1‑149 (MQ=40)
                                                                                                                                                              gcgggcgTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTAtt                                                                                                                                                                                                                                            <  5:1606582/149‑1 (MQ=40)
                                                                                                                                                                     tCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTc                                                                                                                                                                                                                                     >  8:1663343/1‑149 (MQ=255)
                                                                                                                                                                        tGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTgg                                                                                                                                                                                                                                  >  7:281541/1‑149 (MQ=255)
                                                                                                                                                                        tGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAAt                                                                                                    <  2:537613/275‑1 (MQ=255)
                                                                                                                                                                              ccccccTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGGGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCACCAGACAAATTCTTTTCTAATCTGcc                                                                                               <  2:163418/271‑1 (MQ=255)
                                                                                                                                                                              cccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATGAAACAAAAGGCCCAGTCTTCCGACTGAGCCt                                                                                                                                                                                                                                                       <  7:1229102‑M1/34‑1 (MQ=255)
                                                                                                                                                                              cccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCt                                                                                                                                                                                                                                                       >  8:1229102‑M1/89‑122 (MQ=255)
                                                                                                                                                                                aCCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTccc                                                                                                                                                                                                                          >  7:1090905/1‑149 (MQ=255)
                                                                                                                                                                                 cccTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCTCAGTCGGAAGACTGGGCCttt                                                                                                                                                                                                                                                     <  7:86686/121‑1 (MQ=14)
                                                                                                                                                                                    tCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGc                                                                                                                                                                                                                                                                           <  7:637443/96‑1 (MQ=255)
                                                                                                                                                                                          ccGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCa                                                                                                                                                                                                                >  8:163979/1‑149 (MQ=255)
                                                                                                                                                                                            agaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTGAGc                                                                                                                                                                                                                                                         >  5:12710‑M1/75‑106 (MQ=255)
                                                                                                                                                                                                cTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGga                                                                                                                                                                                                          >  7:1007179/1‑149 (MQ=255)
                                                                                                                                                                                                     cAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGAccc                                                                                                                                                                                                     >  5:1075857/1‑149 (MQ=255)
                                                                                                                                                                                                     cAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAACCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGAccc                                                                                                                                                                                                     >  5:1216442/1‑149 (MQ=255)
                                                                                                                                                                                                      aaCGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGAcccc                                                                                                                                                                                                    >  8:1249435/1‑149 (MQ=255)
                                                                                                                                                                                                      aaCGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGAcccc                                                                                                                                                                                                    >  7:387569/1‑149 (MQ=255)
                                                                                                                                                                                                           tCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACt                                                                                                                                                                                               >  6:511505/1‑149 (MQ=255)
                                                                                                                                                                                                             aaaTCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTAc                                                                                                                                                                                             >  6:80150/1‑149 (MQ=255)
                                                                                                                                                                                                              aaTCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACt                                                                                                                                                                                                                                                             >  8:1430977/1‑84 (MQ=255)
                                                                                                                                                                                                              aaTCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACt                                                                                                                                                                                                                                                             <  7:1430977/84‑1 (MQ=255)
                                                                                                                                                                                                               agccacTCCCGGCGGATTTGTCCCAACACGGAGGGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGGCGTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCa                                                              <  4:272602/269‑1 (MQ=17)
                                                                                                                                                                                                                tCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCa                                                                                                                                                                                                                                >  8:1627837/1‑111 (MQ=255)
                                                                                                                                                                                                                tCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCa                                                                                                                                                                                                                                <  7:1627837/111‑1 (MQ=255)
                                                                                                                                                                                                                tCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCacac                                                                                                                                                                                                <  5:418309/143‑1 (MQ=255)
                                                                                                                                                                                                                tCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCacac                                                                                                                                                                                                >  6:418309/1‑143 (MQ=255)
                                                                                                                                                                                                                tCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCAt                                                                                                                                                                                          <  6:1034680/149‑1 (MQ=255)
                                                                                                                                                                                                                tCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCg                                                        <  1:280655/279‑1 (MQ=255)
                                                                                                                                                                                                                   gCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCgg                                                                                                                                                                                       >  7:972339/1‑149 (MQ=255)
                                                                                                                                                                                                                   atggtagtgtggggtctccccatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTTCGACTGAGCctttc                                                                                                                                                                                                                                                    <  5:860784‑M1/37‑2 (MQ=255)
                                                                                                                                                                                                                       ccGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGACAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGAc                                                                                                                                                                                                       <  8:913490/129‑1 (MQ=255)
                                                                                                                                                                                                                        cggcggATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGa                                                                                                                                                                                                                                         >  6:1378872/1‑94 (MQ=255)
                                                                                                                                                                                                                                    ccTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTc                                                                                                                                                                                                                                     >  5:1489314/1‑86 (MQ=255)
                                                                                                                                                                                                                                    ccTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTc                                                                                                                                                                                                                                     <  6:1489314/86‑1 (MQ=255)
                                                                                                                                                                                                                                      tACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTc                                                                                                                                                                    <  5:50921/149‑1 (MQ=255)
                                                                                                                                                                                                                                      tACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTc                                                                                                                                                                    >  7:873077/1‑149 (MQ=255)
                                                                                                                                                                                                                                      ccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTGa                                                                                                                                                                                                                                                           <  8:139277‑M1/30‑1 (MQ=255)
                                                                                                                                                                                                                                      ccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTGa                                                                                                                                                                                                                                                           >  7:139277‑M1/33‑62 (MQ=255)
                                                                                                                                                                                                                                       aCTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCt                                                                                                                                                                   <  7:1641454/149‑1 (MQ=255)
                                                                                                                                                                                                                                         tCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATgg                                                                                                                                                                                                             >  7:332781/1‑105 (MQ=255)
                                                                                                                                                                                                                                           gggAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGt                                                                                                                                                               >  6:887627/1‑149 (MQ=255)
                                                                                                                                                                                                                                             gagaGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACCTCTGAGTTc                                                                                                                                                             >  7:178486/1‑149 (MQ=255)
                                                                                                                                                                                                                                                 gtgtTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCa                                                                                                                                                         >  5:761735/1‑149 (MQ=40)
                                                                                                                                                                                                                                                  tgtTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATggg                                                                                                                                                     >  4:204664/1‑152 (MQ=39)
                                                                                                                                                                                                                                                    ggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTCCGACTGAGCCtt                                                                                                                                                                                                                                                      <  7:528990‑M1/35‑1 (MQ=255)
                                                                                                                                                                                                                                                     tCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATggg                                                                                                                                                     <  5:1560706/149‑1 (MQ=40)
                                                                                                                                                                                                                                                      cACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCa                                                                                                                                                 <  2:204667/152‑1 (MQ=38)
                                                                                                                                                                                                                                                       aCCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGt                                                                                                                                                   >  5:331203/1‑149 (MQ=255)
                                                                                                                                                                                                                                                       aCCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGt                                                                                                                                                   >  8:408248/1‑149 (MQ=40)
                                                                                                                                                                                                                                                           acaaacaaCAGATGAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATg                                                                                                                                                                                                                                       <  7:130396/61‑1 (MQ=255)
                                                                                                                                                                                                                                                           acaaacaaCAGATGAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATg                                                                                                                                                                                                                                       >  8:130396/1‑61 (MQ=255)
                                                                                                                                                                                                                                                               acaacaGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTggg                                                                                                                                           <  8:836722/149‑1 (MQ=40)
                                                                                                                                                                                                                                                                caacaGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGa                                                                                                                                          >  6:642390/1‑149 (MQ=39)
                                                                                                                                                                                                                                                                caacaGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGa                                                                                                                                          >  6:1319255/1‑149 (MQ=40)
                                                                                                                                                                                                                                                                 aacaGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGa                                                                                                                                                                 <  5:365881/125‑1 (MQ=255)
                                                                                                                                                                                                                                                                 aacaGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGa                                                                                                                                                                 >  6:365881/1‑125 (MQ=40)
                                                                                                                                                                                                                                                                 aacaGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGac                                                                                                                                         >  6:488833/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                  acaGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGacc                                                                                                                                        >  8:1202820/1‑149 (MQ=39)
                                                                                                                                                                                                                                                                     aATAAAACAAAAGGCTCAGTCGGAAGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGaccacc                                                                                                                                     <  8:585356‑M1/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                      caaaacCAAAAGGCCACGTTTTCCGTCTGAGCTNTTTTTTTAATTTGTTGGTGGGCGAANCCCGACTCTCGCATGGGGAGACCCCACACTACCAGCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGNGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAg  <  3:167811/277‑1 (MQ=25)
                                                                                                                                                                                                                                                                      ATAAAACAAAAGGCTCAGTCGGAAGACTGGGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACcg                                                                                                                                    <  7:795641‑M1/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                      ATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAgtggt                                                                         <  2:472590‑M1/208‑1 (MQ=255)
                                                                                                                                                                                                                                                                         aaaCAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCgcgc                                                                                                                                 >  7:1243117/1‑149 (MQ=40)
                                                                                                                                                                                                                                                                         aaaCAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCgcgc                                                                                                                                 <  8:1161953/149‑1 (MQ=40)
                                                                                                                                                                                                                                                                          aaCAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCt                                                                                                                                >  5:661791/1‑149 (MQ=39)
                                                                                                                                                                                                                                                                           aCAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTa                                                                                                                               >  5:99303/1‑149 (MQ=40)
                                                                                                                                                                                                                                                                           aCAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTCCCAAGGGTTAGCTCTACCAACTGAACCACATCAg  >  2:359146/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                            cAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTAc                                                                                                                              <  8:873077/149‑1 (MQ=39)
                                                                                                                                                                                                                                                                             aaaaGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACt                                                                                                                             >  8:643722/1‑149 (MQ=39)
                                                                                                                                                                                                                                                                             aaaaGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACt                                                                                                                             >  8:825367/1‑149 (MQ=38)
                                                                                                                                                                                                                                                                             aaaaGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACt                                                                                                                             >  5:636828/1‑149 (MQ=38)
                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                               
GAGCAAGCTGCACCTGATGAATTTTTTGCGCATCCTAAATCAGAGCGTACGAGGGCATTTTTATCGCAGGTAATCCATTAATTGAATGTTAGTTCGAAAAGCAAAAAGGCCATCCTTTCGGATGGCCTTTCGCTTGATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAATTCTTTTCTAATCTGCCGAACTTTAACCTAAAAAGTGGTGCTGATACCCAGAGTCGAACTGGGGACCTCACCCTTACCAAGGGTGCGCTCTACCAACTGAGCCATATCAG  >  NC_000913/3423114‑3423654

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: