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%	82.3 / 18.0	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

CAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CT‑AACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGG  >  NC_000913/15232‑15532                                                                                                                                                       |                                                                                                                                                          caacaGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGa                                                                                                                                                           <  2:988813/151‑1 (MQ=255)   acaGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATc                                                                                                                                                         <  1:760311/151‑1 (MQ=255)      gCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCAc                                                                                                                                                       >  2:679900/1‑150 (MQ=255)               cAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAg                                                                                                                                             <  1:855979/150‑1 (MQ=255)                 aaGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCCT‑aa                                                                                                                                              >  1:237411/1‑147 (MQ=255)                     gCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCACTCCCCC‑AAGCG‑GT‑aa                                                                                                                                     >  1:997668/1‑150 (MQ=255)                        tcttTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CT‑AACtt                                                                                                                                  <  1:532405/150‑1 (MQ=255)                        tcttTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CT‑AACtt                                                                                                                                  <  2:773374/150‑1 (MQ=255)                        tcttTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CT‑AACt                                                                                                                                   <  1:298794/149‑1 (MQ=255)                              aTGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CT‑AACTTAAGGGtt                                                                                                                           <  1:450090/151‑1 (MQ=255)                                   gtgAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CT‑AACTTAAGGGTTGTGGt                                                                                                                      <  2:334212/151‑1 (MQ=255)                                     gaagaaGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CT‑AACTTAAGGGTTGTGGTa                                                                                                                     <  2:940305/150‑1 (MQ=255)                                     gaagaaGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CT‑AACTTAAGGGTTGTGGTAt                                                                                                                    <  2:974152/151‑1 (MQ=255)                                       agaagTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCAACCCCCT‑AAACG‑CTAAACTAAAGGGtggggggtt                                                                                                                    >  2:928789/1‑142 (MQ=255)                                        gaagTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGACAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCCATTCCCT‑AAACGCCT‑AATTTAAGGGtgtgg                                                                                                                        >  1:854883/1‑141 (MQ=255)                                             tttttGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CT‑AACTTAAGGGTTGTGGTATTGCGCCt                                                                                                             <  2:589327/150‑1 (MQ=255)                                             tttttGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGAACACTCCCCC‑AAAAG‑CC‑CACTTAAGGGGTGTGGGATTATGCCt                                                                                                             >  1:617834/1‑150 (MQ=255)                                               tttGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACCCCCCT‑AAACG‑CT‑TAATTAAGGGGGGTggggtttagcccgtt                                                                                                          >  1:895848/1‑138 (MQ=255)                                                     gacCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACCCCCCT‑CAACG‑CC‑ACATTTAGGGGTGTGGGATTAAGCCTGTTATGa                                                                                                      >  2:712948/1‑149 (MQ=255)                                                        cTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCCT‑AAACG‑GT‑AACTTAAGGGGTGGGGTGTTAAGCCTGattttat                                                                                                     >  2:839839/1‑142 (MQ=255)                                                          gACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CT‑AACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACg                                                                                               <  2:309300/151‑1 (MQ=255)                                                              cGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGACACCCCCCC‑AAGCG‑CT‑AACTTTAGGGGTGTGGGGTTATGCCCGATATGAttt                                                                                                   >  1:9738/1‑143 (MQ=255)                                                               gCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CT‑AACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCg                                                                                          <  1:993872/151‑1 (MQ=255)                                                                cTAACCTCCCCAAAAGCCTTCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CT‑AACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCg                                                                                          <  2:721125/150‑1 (MQ=255)                                                                  aaCCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CT‑AACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGAt                                                                                        <  1:834486/150‑1 (MQ=255)                                                                   aCCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGATACTCCCAT‑AAACG‑CT‑AAATTTAAGGGTGTGGGGTTACGCCTGATATGATTTTAAGTGCCGAt                                                                                        >  1:680961/1‑149 (MQ=255)                                                                   aCCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CT‑AACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGa                                                                                      <  2:1006380/151‑1 (MQ=255)                                                                       ccccAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCACTCCCAT‑AAACG‑CG‑AACTTAAGGGGTGTGGGATTTCGCCTGGTATGATTTTAAGTGCCGGTg                                                                                       >  2:167739/1‑146 (MQ=255)                                                                       ccccAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CT‑AACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTa                                                                                  <  1:723505/151‑1 (MQ=255)                                                                        cccAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CT‑AACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCTATGAATTAc                                                                                 <  2:305630/151‑1 (MQ=255)                                                                         ccAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CT‑AACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAct                                                                                <  1:318493/151‑1 (MQ=255)                                                                         ccAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CT‑AACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAct                                                                                <  2:948297/151‑1 (MQ=255)                                                                                 cTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGGACTCCCCTCAAACC‑CA‑AACTT‑TGGGGTGTGGGGTTAAGCCCGGTATGATTTAAAAGGCCGATGAATTAATCTCCCggt                                                                        >  2:658576/1‑149 (MQ=255)                                                                                 cTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCCC‑AAACG‑CC‑AAATTTAGGGGTGTGGGATTACGCCTGGTATGATTTTAAGTGCCGATGAAATACTCTCACg                                                                          >  1:171782/1‑149 (MQ=255)                                                                                      cGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCAT‑AAGAG‑CT‑CACTTAAGGGGTGTGGTATTAATCCTGAGTTTATTTTAATTGCCGATGAATTTcttaccagcaaactg                                                                   >  2:908104/1‑138 (MQ=255)                                                                                      cGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CT‑AACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTg                                                                   <  1:886747/151‑1 (MQ=255)                                                                                         gggCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACCCCCCT‑AAACG‑CC‑AACTTAAGGGGTGTGGGATTACGCCCGATATGATTTTACGTGCCGATGAATTACTCTCACGCTAAATgg                                                                  >  1:578755/1‑149 (MQ=25)                                                                                            cAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCACTCCCCC‑AAACG‑CC‑AACTTTAGGGGTGTGGGATTTTGCCCGATATGATTTAAAGTGCCGCTGAAATAATCTCACGATAAATGGTCAGc                                                             >  1:893874/1‑151 (MQ=255)                                                                                              ggCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CT‑AACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCaa                                                           <  2:771742/151‑1 (MQ=255)                                                                                                ccTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAACG‑CT‑AAATTTAGGGGTGTGGGATTAAGCCCGATATGATTTTACGTGCCGGTGAATTAATCTCACCATAACTGGGCAGCAAtt                                                         >  1:173395/1‑151 (MQ=255)                                                                                                 cTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGCACTCCCCT‑AAGCG‑CC‑AACTTAAGGGGTGTGGGATTACGCCTGATATGATTTTAAGTGCCGATGAATTACTCTCACGATAAATGGTCAGCaaat                                                         >  1:882135/1‑148 (MQ=21)                                                                                                       aaaaaTAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCCT‑AAAAG‑CC‑AACTTAAAGGGTGTGGGATTTTGCCCGATTTGATTTTAAGTGCCGGTGAATTATTCTCACCATAAATTGTCAGCCATTTTTGccc                                                  >  2:358896/1‑151 (MQ=255)                                                                                                            tAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCCT‑AAACG‑GC‑AACTTAAGGGTTGTGGTGTTTTGCCTGGTTTGATTTTACGTGCCGGTGAATTAATCTCACCATAAATGGTCAGGAATTCTTGCCCCTAtt                                             >  2:465186/1‑151 (MQ=255)                                                                                                              gggTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCCT‑AAACG‑CC‑AACTTTAAGGGTGTGGTATTACGCCTGATATGATTTAAAGTGCCGATGAATTAATCTCACGATAACTTGGCAGCAATTCTTGCCCATAtt                                             >  1:883846/1‑149 (MQ=18)                                                                                                                  gCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CT‑AACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTaa                                        <  1:920626/150‑1 (MQ=255)                                                                                                                     ttGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGAACTCCCAT‑AAACG‑CT‑AAATTAAGGGGTGTGGGATTACGCCCGATATGATTTAAAGTGCCGGTGAATTAATCTCAAGATAAATGGGCAGCAAATCTGGCCCCTATTTGTAAGacc                                    >  1:606002/1‑151 (MQ=255)                                                                                                                       gAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCCT‑AAAAG‑CC‑AAATTAAGGGGTGTGGGGTTATGCCTGATATGTTTTAAAGTGCCGAGGAAATACTCTCACGACAACAGGTCAGCAATTCTGGCCCCTATTGGTAAGccc                                    >  2:504434/1‑149 (MQ=255)                                                                                                                        aaGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CT‑AACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCgaa                                 <  1:560668/151‑1 (MQ=255)                                                                                                                          gATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CT‑AACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCgaaga                               <  1:637032/151‑1 (MQ=255)                                                                                                                            tatGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCCT‑AAGCG‑CT‑AACTTTAGGGGTGTGGGGTTACGCCTGGTATGATTTTAAGTGCCGATGAATTACTCTCACGATAAATGGGCAGCAAATCTTGGCCCTATTGGTAAGCCCgaagaa                              >  1:59678/1‑150 (MQ=17)                                                                                                                             atGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCAT‑AAACG‑CT‑AACTTTAAGGGTGGGGTATTTAGCCCGGTATGATTTTAAGTGCCGATGAATTATTCTCCCGATAAATGGTCAGCAATTCTTGCCCCTATTTGTAAGCCCgaagaa                              >  1:700356/1‑149 (MQ=255)                                                                                                                                      ccTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CT‑AACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTc                   <  2:171777/151‑1 (MQ=255)                                                                                                                                           aaaGTGGCGGGGATCACTCCCAT‑AAGCG‑CT‑AACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCa               <  2:237405/150‑1 (MQ=39)                                                                                                                                            aaGTGGCGGGGATCACTCCCAT‑AAGCG‑CT‑AACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACg             <  2:179550/151‑1 (MQ=38)                                                                                                                                               tGGCGGGGATCACTCCCAT‑AAGCG‑CT‑AACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTa           >  2:85608/1‑150 (MQ=35)                                                                                                                                                    gggATCACTCCCAT‑AAGCG‑CT‑AACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCg     <  2:866597/151‑1 (MQ=33)                                                                                                                                                       aTCACTCCCAT‑AAGCG‑CT‑AACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCgggg  >  2:123172/1‑151 (MQ=33)                                                                                                                                                       |                                                                                                                                                          CAACAGCCCGCGCTCAAAGAGCTTCTTTGATGGTGTGAAGAAGTTTTTTGACGACCTGACCCGCTAACCTCCCCAAAAGCCTGCCCGTGGGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCAT‑AAGCG‑CT‑AACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGG  >  NC_000913/15232‑15532

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