Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,423,398 C→T 10.8% intergenic (+204/+25) 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,3980CT10.8% 216.0 / 15.4 74intergenic (+204/+25)yhdZ/rrfFputative amino acid ABC transporter ATPase/5S ribosomal RNA of rrnD operon
Reads supporting (aligned to +/- strand):  ref base C (37/29);  new base T (5/3);  total (42/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.32e-01

ATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAA  >  NC_000913/3423250‑3423544
                                                                                                                                                    |                                                                                                                                                  
atttgtctggcggcagtagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcagATAAAACAAAAGGCCCAGTCTTc                                                                                                                                                    <  4:457692‑M1/23‑1 (MQ=255)
 tttgtctggcggcagtagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTcc                                                                                                                                                   <  1:697432‑M1/24‑1 (MQ=255)
      ctggcggcagtagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctcctcatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTCCGACTg                                                                                                                                              <  3:359411‑M1/29‑1 (MQ=255)
        ggcggcagtagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcagATAAAACAAAAGGCCCAGTCTTCCGACTGAg                                                                                                                                            >  3:457692‑M1/119‑149 (MQ=255)
         gcggcagtagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcagATAAAACAAAAGGCCCAGTCTTCCGACTGAGc                                                                                                                                           <  4:641411‑M1/32‑1 (MQ=255)
           ggcagtagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctcctcatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTCCGACTGAGCCt                                                                                                                                         >  4:359411‑M1/116‑149 (MQ=255)
           ggcagtagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCt                                                                                                                                         >  3:1331794‑M1/116‑149 (MQ=255)
           ggcagtagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCt                                                                                                                                         >  1:399315‑M1/116‑149 (MQ=255)
            gcagtagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTTCGACTGAGCCtt                                                                                                                                        >  1:456454‑M1/115‑149 (MQ=255)
            gcagtagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCtt                                                                                                                                        >  2:697432‑M1/115‑149 (MQ=255)
                tagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcagATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGt                                                                                                                                    >  3:641411‑M1/111‑149 (MQ=255)
                tagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTGAGc                                                                                                                                           <  4:1378499‑M1/32‑1 (MQ=255)
                tagcgcggtggtcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTGAGc                                                                                                                                           >  3:1378499‑M1/111‑142 (MQ=255)
                           tcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCtttt                                                                                                                                      >  1:752766‑M1/100‑136 (MQ=255)
                           tcccacctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCtttt                                                                                                                                      <  2:752766‑M1/37‑1 (MQ=255)
                                 ctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctcctcatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTCCGACTGAGCCt                                                                                                                                         <  3:327972‑M1/34‑1 (MQ=255)
                                 ctgaccccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctcctcatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTCCGACTGAGCCt                                                                                                                                         >  4:327972‑M1/94‑127 (MQ=255)
                                      cccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCt                                                                                                                                         <  3:1428409‑M1/34‑1 (MQ=255)
                                      cccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCt                                                                                                                                         >  4:1428409‑M1/89‑122 (MQ=255)
                                       ccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTTCGACTGAGc                                                                                                                                           >  4:429407‑M1/88‑119 (MQ=255)
                                       ccatgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTTCGACTGAGc                                                                                                                                           <  3:429407‑M1/32‑1 (MQ=255)
                                        catgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCt                                                                                                                                         <  2:234009‑M1/34‑1 (MQ=255)
                                        catgccgaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCt                                                                                                                                         >  1:234009‑M1/87‑120 (MQ=255)
                                           cTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTc                                                                                                                                                    <  3:1029105/106‑1 (MQ=255)
                                           cTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTc                                                                                                                                                    >  4:1029105/1‑106 (MQ=40)
                                              gaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctcctcatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTCCGACTGAGcc                                                                                                                                          <  3:168154‑M1/33‑1 (MQ=255)
                                              gaactcagaagtgaaacgccgtagcgccgatggtagtgtggggtctcctcatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTCCGACTGAGcc                                                                                                                                          >  4:168154‑M1/81‑113 (MQ=255)
                                                  tcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccaggcatcaaATAAAACGAGAGGCCCAGTCTTTCGACTGAGc                                                                                                                                           >  3:178072‑M1/77‑108 (MQ=255)
                                                  tcagaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccaggcatcaaATAAAACGAGAGGCCCAGTCTTTCGACTGAGc                                                                                                                                           <  4:178072‑M1/32‑1 (MQ=255)
                                                    agaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCtttt                                                                                                                                      >  1:747899‑M1/75‑111 (MQ=255)
                                                    agaagtgaaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCtttt                                                                                                                                      <  2:747899‑M1/37‑1 (MQ=255)
                                                     ttGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGtt                                                                                                                                   <  2:1431621/113‑1 (MQ=255)
                                                     ttGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGtt                                                                                                                                   >  1:1431621/1‑113 (MQ=255)
                                                         ttCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTc                                                                                                                                                    <  3:977342/92‑1 (MQ=255)
                                                         ttCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTc                                                                                                                                                    >  4:977342/1‑92 (MQ=255)
                                                           aaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTTCGACTGAGc                                                                                                                                           <  1:1116491‑M1/32‑1 (MQ=255)
                                                           aaacgccgtagcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTTCGACTGAGc                                                                                                                                           >  2:1116491‑M1/68‑99 (MQ=255)
                                                             acgccgtagcgccgatggtagtgtggggtctcctcatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTCCGACTGAGCCtt                                                                                                                                        <  4:1410685‑M1/35‑1 (MQ=255)
                                                             acgccgtagcgccgatggtagtgtggggtctcctcatgcgagagtagggaactgccaggcatcaaATAAAACGAAAGGCCCAGTCTTCCGACTGAGCCtt                                                                                                                                        >  3:1410685‑M1/66‑100 (MQ=255)
                                                                 gTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCttt                                                                                                                                       <  3:791340/97‑1 (MQ=255)
                                                                 gTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCttt                                                                                                                                       >  4:791340/1‑97 (MQ=255)
                                                                     gcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcagATAAAACAAAAGGCCCAGTCTTCCGACTGAGcc                                                                                                                                          <  2:865828‑M1/33‑1 (MQ=255)
                                                                     gcgccgatggtagtgtggggtctccccatgcgagagtagggaactgccagacatcagATAAAACAAAAGGCCCAGTCTTCCGACTGAGcc                                                                                                                                          >  1:865828‑M1/58‑90 (MQ=255)
                                                                      aaTCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGAc                                                                                                                                                >  4:1036186/1‑83 (MQ=255)
                                                                      aaTCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGAc                                                                                                                                                <  3:1036186/83‑1 (MQ=255)
                                                                       aTCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGAc                                                                                                                                                <  3:260999/82‑1 (MQ=255)
                                                                       aTCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGAc                                                                                                                                                >  4:260999/1‑82 (MQ=255)
                                                                                 ggcggATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTAt                                                                                                                               <  1:992421/89‑1 (MQ=255)
                                                                                 ggcggATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTAt                                                                                                                               >  2:992421/1‑89 (MQ=255)
                                                                                 ggcggATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTcc                                                                                                             <  2:253651/107‑1 (MQ=255)
                                                                                 ggcggATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTcc                                                                                                             >  1:253651/1‑107 (MQ=255)
                                                                                  gcggATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACg                                                                  >  1:1442283/1‑149 (MQ=255)
                                                                                             cccatgcgagagtagggaactgccagacatcagATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCtt                                                                                                                                        >  2:157606‑M1/34‑68 (MQ=255)
                                                                                             cccatgcgagagtagggaactgccagacatcagATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCtt                                                                                                                                        <  1:157606‑M1/35‑1 (MQ=255)
                                                                                                 tCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGa                                                   <  2:1442283/149‑1 (MQ=255)
                                                                                                     agagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCtttt                                                                                                                                      <  2:1372437‑M1/37‑1 (MQ=255)
                                                                                                     agagtagggaactgccagacatcaaATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCtttt                                                                                                                                      >  1:1372437‑M1/26‑62 (MQ=255)
                                                                                                             tCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGAc                                                                                                                                                <  1:734708/44‑1 (MQ=40)
                                                                                                             tCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGAc                                                                                                                                                >  2:734708/1‑44 (MQ=40)
                                                                                                                       acaacaGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAg                                  <  3:1346152/144‑1 (MQ=40)
                                                                                                                       acaacaGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAg                                  >  4:1346152/1‑144 (MQ=40)
                                                                                                                               tAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATggg                                                                                              >  1:539010/1‑76 (MQ=255)
                                                                                                                               tAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATggg                                                                                              <  2:539010/76‑1 (MQ=255)
                                                                                                                                aaaaCAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATgg                                        >  3:1303368/1‑129 (MQ=255)
                                                                                                                                aaaaCAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATgg                                        <  4:1303368/129‑1 (MQ=255)
                                                                                                                                     aaaaGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACt               >  1:239662/1‑149 (MQ=38)
                                                                                                                                            ccAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCGTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCt                                                     >  1:948455/1‑104 (MQ=255)
                                                                                                                                            ccAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGga                                                                                            <  1:1162211/65‑1 (MQ=255)
                                                                                                                                            ccAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGga                                                                                            >  2:1162211/1‑65 (MQ=255)
                                                                                                                                            ccAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTc                                                           <  3:871082/98‑1 (MQ=255)
                                                                                                                                            ccAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCt                                                     <  2:948455/104‑1 (MQ=40)
                                                                                                                                            ccAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCa        >  1:74604/1‑149 (MQ=38)
                                                                                                                                            ccAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCa        >  3:854011/1‑149 (MQ=38)
                                                                                                                                            ccAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCAGCGGCGCTACGGCGTTTc                                                           >  4:871082/1‑98 (MQ=255)
                                                                                                                                              aGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCca                                                                                     >  3:655388/1‑70 (MQ=255)
                                                                                                                                              aGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCca                                                                                     <  4:655388/70‑1 (MQ=255)
                                                                                                                                                  ttCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACaaa  >  3:370527/1‑149 (MQ=40)
                                                                                                                                                  ttCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACaaa  >  2:112605/1‑149 (MQ=40)
                                                                                                                                                  ttCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACaaa  >  4:699254/1‑149 (MQ=40)
                                                                                                                                                  ttCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACaaa  <  4:854011/149‑1 (MQ=39)
                                                                                                                                                    ccGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTa                                                                    <  2:190515/81‑1 (MQ=40)
                                                                                                                                                    ccGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTa                                                                    >  1:190515/1‑81 (MQ=40)
                                                                                                                                                    |                                                                                                                                                  
ATTTGATGTCTGGCAGTTTATGGCGGGCGTCCTGCCCGCCACCCTCCGGGCCGTTGCTTCGCAACGTTCAAATCCGCTCCCGGCGGATTTGTCCTACTCGGGAGAGTGTTCACCGACAAACAACAGATAAAACAAAAGGCCCAGTCTTCCGACTGAGCCTTTTGTTTTATTTGATGTCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACTGCCGCCAGACAAA  >  NC_000913/3423250‑3423544

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