Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	15,382	0	A	G	41.8%	86.0 / 16.1	55	intergenic (+84/‑63)	dnaJ/insL1	chaperone Hsp40, DnaK co‑chaperone/IS186 transposase
Reads supporting (aligned to +/- strand):  ref base A (4/28);  new base G (23/0);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 9.40e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACT‑‑CCCATAAGCGCTAACTTAAGGG‑TTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCA‑GCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGG  >  NC_000913/15246‑15532                                                                                                                                         |                                                                                                                                                           cAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACC‑‑CCCAt                                                                                                                                                >  2:29782/1‑147 (MQ=255) cAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACT‑‑CCCATaa                                                                                                                                              <  2:778650/149‑1 (MQ=255)  aaaGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCACT‑‑CCCACaa                                                                                                                                              >  1:947915/1‑148 (MQ=255)  aaaGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACT‑‑CCCATAAgc                                                                                                                                            <  2:508549/150‑1 (MQ=255)   aaGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGGACT‑‑CCCCTAAacg                                                                                                                                           >  1:957883/1‑150 (MQ=255)        cttcttTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCAAT‑‑CCCCTAAtcgccca                                                                                                                                       >  1:93149/1‑146 (MQ=255)          tcttTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACT‑‑CCCATAAGCGCTAACtt                                                                                                                                    <  2:713840/150‑1 (MQ=255)                aTGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCACT‑‑CCCCTAAACGGCTACTTaaag                                                                                                                                >  1:260876/1‑146 (MQ=255)                aTGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCAAT‑‑CCCAACAAAGCTAACTTAAagg                                                                                                                               >  2:361294/1‑149 (MQ=255)                  ggTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTAACG‑‑ACCAAAAAAGCCAACTTAAagg                                                                                                                               >  1:675468/1‑147 (MQ=255)                  ggTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACT‑‑CCCATAAGCGCTAACTTAAGGG‑Ttg                                                                                                                           <  1:573113/150‑1 (MQ=255)                     gtgAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACT‑‑CCCATAAGCGCTAACTTAAGGG‑TTGTgg                                                                                                                        <  2:458381/150‑1 (MQ=255)                     gtgAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACT‑‑CCCATAAGCGCTAACTTAAGGG‑TTGTGGt                                                                                                                       <  2:334650/151‑1 (MQ=255)                      tgAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTAACT‑‑CCCCTAAACGCCAAATTAAAGG‑Gtg                                                                                                                           >  1:922494/1‑146 (MQ=255)                        aagaagTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACT‑‑CCCATAAGCGCTAACTTAAGGG‑TTGTGGTAtt                                                                                                                    <  1:1151739/151‑1 (MQ=255)                        aagaagTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACT‑‑CCCATAAGCGCTAACTTAAGGG‑TTGTGGTAtt                                                                                                                    <  1:610628/151‑1 (MQ=255)                         agaagTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACT‑‑CCCATAAGCGCTAACTTAAGGG‑TTGTGGTATTa                                                                                                                   <  1:804984/151‑1 (MQ=255)                          gaagTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACT‑‑CCCATAAGCGCTAACTTAAGGG‑TTGTGGTATTAc                                                                                                                  <  2:976733/151‑1 (MQ=255)                                    gacgacCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACT‑‑CCCATAAGCGCTAACTTAAGGG‑TTGTGGTATTACGCCTGatat                                                                                                         <  2:716667/150‑1 (MQ=255)                                    gacgacCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACT‑‑CCCATAAGCGCTAACTTAAGGG‑TTGTGGTATTACGCCTGATATg                                                                                                        <  1:103884/151‑1 (MQ=255)                                      cgacCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACT‑‑CCCATAAGCGCTAACTTAAGGG‑GTGTGGTGTTACGCCTGatat                                                                                                         >  2:850436/1‑148 (MQ=255)                                      cgacCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGGACC‑‑CCCACAAAAGCCAACTTTAGGG‑GTGTggggt                                                                                                                     >  1:974609/1‑133 (MQ=255)                                       gacCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACT‑‑CCCCTAAACGCTAACTTTAAGG‑TTGTGGGATTACGCCTGATATGAtt                                                                                                     >  2:353995/1‑151 (MQ=255)                                         ccTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACT‑‑CCCATAAGCGCTCAATTAAGGG‑GTGTGGGATTTAGCCTGATATTAttt                                                                                                    >  2:185008/1‑150 (MQ=255)                                                cGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTGACT‑‑CCCACAAACGCCAAATTAAAGG‑GTTTGGGATTTCGCCTGATATTATTTTAAGTGcc                                                                                            >  1:614733/1‑151 (MQ=255)                                                 gCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACT‑‑CCCATAAGCGCTAACTTAAGGG‑TTGTGGTATTACGCCTGATATGATTTAACGTGCCg                                                                                           <  2:370721/151‑1 (MQ=255)                                                   tAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCACT‑‑TCCATAAAACCTCACCTTAAGG‑TTGTGGGAGTTTGACTTGTATTATTTAACg                                                                                                >  1:1103539/1‑144 (MQ=255)                                                    aaCCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACT‑‑CCCATAAGCGCTAACTTAAGGG‑TTGTGGTATTACGCCTGATATGATTTAACGTGCCGAt                                                                                         <  1:645778/150‑1 (MQ=255)                                                    aaCCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAAGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACT‑‑CCCATAAGCGCTAACTTAAGGG‑TTGTGGTATTACGCCTGATATGATTTAACGTGCCGAt                                                                                         <  1:392164/150‑1 (MQ=255)                                                       cTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACT‑‑CCCACAAACGCTAACTTAAAGG‑GTGTGGGATTTCGCCCGGTATGATTTTAAGTGCCGGTGaa                                                                                      >  2:960459/1‑150 (MQ=255)                                                         ccccAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACT‑‑CCCATAAGCGCTAACTTAAGGG‑TTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAAt                                                                                     <  2:1065860/149‑1 (MQ=255)                                                            cAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCACC‑‑CCCCTAAACGCTAACTTAAGGG‑GTGTGGGATTTCGCCCGGTATTATTTAAAGTGCCGATGaa                                                                                      >  2:1133163/1‑145 (MQ=255)                                                               aaGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACT‑‑CCCATAAGCGCTAACTTAAGGG‑TTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCa                                                                             <  1:353360/151‑1 (MQ=255)                                                                 gCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACT‑‑CCCATAAGCGCTAACTTAAGGG‑TTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACg                                                                           <  2:426305/151‑1 (MQ=255)                                                                     gCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACT‑‑CCCATAAGCGCTAACTTAAGGG‑TTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGAt                                                                         <  1:772537/149‑1 (MQ=255)                                                                      cccGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCACC‑‑CCCCTAAGCGCCAAATTAAGGG‑GTGTGGGATTAAGCCTGCTATGATTTTAACTGCCGATGGATTATTCTCTCGa                                                                          >  2:441759/1‑147 (MQ=255)                                                                             gCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACT‑‑CCCATAAGCGCTAACTTAAGGG‑TTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCA‑g                                                              <  2:455194/151‑1 (MQ=255)                                                                                ggCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACT‑‑CCCATAAGCGCTAACTTAAGGG‑TTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCA‑GCaa                                                           <  1:950820/151‑1 (MQ=255)                                                                                         aaaaaTAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGTACC‑‑CCCCCAAACGGTAACTTTAGGG‑GTGTGGGATTACGCCTGGTATTATTTTAAGAGCCGGTGAATTAATCTCACGATAACTAGGTA‑GGAAATCTTGccc                                                  >  1:901103/1‑151 (MQ=255)                                                                                         aaaaaTAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACT‑‑CCCATAAGCGCTAACTTAAGGG‑TTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCA‑GCAATTCTGGccc                                                  <  1:500153/151‑1 (MQ=255)                                                                                                gggTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACT‑‑CCCATAAGCGCTAACTTAAGGG‑TTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCA‑GCAATTCTGGCCCATATTgg                                           <  1:404161/151‑1 (MQ=255)                                                                                                     cGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACT‑‑CCCCTAAGCGCTAACTTAAGGG‑GTGTGGGATTAAGCCTGATATGATTTAAAGTGCCGATGAATTAATCTCACGATAAATGGTCA‑GCAAATCTGGCCCCTATTTGTaa                                        >  2:735729/1‑149 (MQ=255)                                                                                                      gTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGAAAA‑‑CCCCTAAAGGCCCACTTAAGGG‑GTGTGGGATTACGGCCGGTATGAATTTACGAGGCGATGAATTACTTTCACCATAAATGGGCA‑Gcgaattttgccc                                                  >  1:913395/1‑127 (MQ=255)                                                                                                       ttGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGAACTTCCCCAAAAGC‑CTAA‑TTAAGGGTTTTTGGTGTT‑TGGCTGCTATGATTTTAAGGGTCGACGGATTAATCTCCCCCTAACTGGTCAGGCATTTTTGGcccttctggtaaaccc                                     >  1:702764/1‑138 (MQ=255)                                                                                                         gAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACT‑‑CCCATAAGCGCTAACTTAAGGG‑TTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCA‑GCAATTCTGGCCCATATTGGTAAGCCCga                                  <  1:576514/151‑1 (MQ=255)                                                                                                          aaGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACT‑‑CCCCTAAACGCCAAATTAAAGG‑GTGTGGGATTTCGGCTTATATTATTTTAAGTGGCGAAGGATTTATCTCCCGAAAAATGGTCA‑GCAATTTTTGCCCCTTTTTGTaaaac                                     >  1:853376/1‑144 (MQ=255)                                                                                                          aaGATATGCGAGCACCTGTAAAGTGGCGGGGATCACT‑‑CCCATAAGCGCTAACTTAAGGG‑TTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCA‑GCAATTCTGGCCCATATTGGTAAGCCCgaa                                 <  1:29781/151‑1 (MQ=255)                                                                                                             atatGCGAGCACCTGTAAAGTGGCGGGGGACACT‑‑CCCCTAAACGCTAACTTAAGGG‑GTGTGGGGTTATTCCTGGTATGATTTTTAGTGCCGATGAATTAATCCCCCGCTAAATCGTCA‑GCAAATCTTGCCCATATTTGTAAGacc                                    >  1:68442/1‑145 (MQ=255)                                                                                                                   gAGCACCTGTAAAGTGGCGGGGAGTACT‑‑CCCCTAAGAGCCTACTTTAGGG‑GTGGGGGGTTTCGCCCGATATGATTTTAAGTGCCGATGAAATACTCTCACGATAAATGGTCA‑GCAATTCTTGGCCCTATTTGTAAACCCCAAGAACTgg                          >  1:50244/1‑149 (MQ=255)                                                                                                                        ccTGTAAAGTGGCGGGGATCACT‑‑CCCATAAGCGCTAACTTAAGGG‑TTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCA‑GCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTc                   <  2:922517/151‑1 (MQ=255)                                                                                                                           gTAAAGTGGCGGGGATCACT‑‑CCCATAAGCGCTAACTTAAGGG‑TTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCA‑GCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCg                  <  2:110215/149‑1 (MQ=40)                                                                                                                                   gCGGGGATCACT‑‑CCCATAAGCGCTAACTTAAGGG‑TTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCA‑GCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATg        <  2:469893/151‑1 (MQ=34)                                                                                                                                    cGGGGATCACT‑‑CCCATAAGCGCTAACTTAAGGG‑TTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCA‑GCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGc       >  2:228162/1‑151 (MQ=33)                                                                                                                                     ggggATCACT‑‑CCCATAAGCGCTAACTTAAGGG‑TTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCA‑GCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGcc      >  2:1146305/1‑151 (MQ=33)                                                                                                                                     ggggATCACT‑‑CCCATAAGCGCTAACTTAAGGG‑TTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCA‑GCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGc       >  2:60363/1‑150 (MQ=33)                                                                                                                                         aTCACT‑‑CCCATAAGCGCTAACTTAAGGG‑TTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCA‑GCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCgggg  >  1:428780/1‑151 (MQ=33)                                                                                                                                         |                                                                                                                                                           CAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACT‑‑CCCATAAGCGCTAACTTAAGGG‑TTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCA‑GCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGG  >  NC_000913/15246‑15532

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