Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_000913	3,608,682	0	T→C	34.2%	13.5 ‑7.9 0.0	43	intergenic (+33/+69)	zntA/tusA	zinc, cobalt and lead efflux system/mnm(5)‑s(2)U34‑tRNA 2‑thiolation sulfurtransferase
Rejected: Polymorphism prediction not supported by minimum number of reads on both strands.
Rejected: Prediction has biased read strand distribution supporting polymorphism.
Reads supporting (aligned to +/- strand):  new base (0/14):  ref base (15/12):  total (16/27)
Fisher's exact test strand distribution p-value = 3.79e-04
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

ACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTT‑GCG‑CAGGAG‑ATAAGGCAAACC‑GATCGC‑AACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGG‑TTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACG  >  NC_000913/3608475‑3608884                                                                                                                                                                                                                     |                                                                                                                                                                                                            aCGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTT‑GCG‑CAGGAG‑ATAAGGCAAACC‑GATCGC‑AACATTGAGCGCGATCGGTcc                                                                                                                                                                                                           >  2:574802/1‑210 (MQ=255)       cACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTT‑GCG‑CAGGAG‑ATAAGGCAAACC‑GATCGC‑AACATTGAGCGCGATCGGTCCCCTCGc                                                                                                                                                                                                     <  2:526370/210‑1 (MQ=255)       cACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTT‑GCG‑CAGGAG‑ATAAGGCAAACC‑GATCGC‑AACATTGAGCGCGATCGGGCCCCTCGc                                                                                                                                                                                                     >  1:526324/1‑210 (MQ=255)         cTCACGCCAATATCCGCCAGAACATCAATAGTGAGCGGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGTTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTT‑GCG‑CAGGAG‑ATAAGGCAAACC‑GATCGC‑AACATTGAGCGCGATCGGTCCCCTCGccc                                                                                                                                                                                                   <  1:551159/210‑1 (MQ=255)         cTCACGCCAATAGCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGGTCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTT‑GCG‑CAGGAG‑ATAAGGCAAACC‑GATCGC‑AACATTGAGCGCGATCGGTCCCCTCGccc                                                                                                                                                                                                   <  1:461958/210‑1 (MQ=255)            aCGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTT‑GCG‑CAGGAG‑ATAAGGCAAACC‑GATCGC‑AACATTGAGCGCGATCGGTCCCCTCGCCCctc                                                                                                                                                                                                <  2:277824/210‑1 (MQ=255)                 aaTATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTT‑GCG‑CAGGAG‑ATAAGGCAAACC‑GATCGC‑AACATTGAGCGCGATCGGTCCCCTCGCCCCTCTgggg                                                                                                                                                                                           <  2:245105/210‑1 (MQ=255)                   tatCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTT‑GCG‑CAGGAG‑ATAAGGCAAACC‑GATCGC‑AACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGgag                                                                                                                                                                                         <  2:262372/210‑1 (MQ=255)                       cGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGGTAAGATTGTT‑GCG‑CAGGGG‑ATAAGGCAAACC‑GAACGC‑AACATTGAGCGCGATCCGTCCCCTCGCCCCTCCGGGGAGAggg                                                                                                                                                                                     >  1:118372/1‑210 (MQ=255)                                  cACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCCGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATGGTT‑GCG‑CAGGAG‑ATAAGGCAAACC‑GATCGA‑AACATTGAGCGCGATCGGTCCCCTCGCCCCCCTGGGGAGAGGG‑CTAGGGCGAgg                                                                                                                                                                         >  1:183367/1‑210 (MQ=255)                                    cTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGACTGCGTTAAGATTGTT‑GCG‑CAGGAG‑ATAAGGCAAACC‑GCTCGC‑AACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGG‑TTCGGGTGAgggg                                                                                                                                                                       >  1:457488/1‑210 (MQ=255)                                      aTTGCGCTGGGGCTGAAAGCGATCTTCCTCGTCACCACGCTGTTAGGTATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACCGGGGCGACGGTGCTGGTGACACCGAATGCGTTAAGATTGTT‑GCG‑CAGGAG‑ATAAGTCGAACC‑GATCGT‑AATATTGACCGCGCTCACTCCCCTCTCCCCCCTGTGCAGCgggcccggtatcgggtat                                                                                                                                                                      >  1:586447/1‑195 (MQ=255)                                                ggCTAAAAGGGATCTTCCTGGTCACAACGTTGTTGGGGTGGACCGGGTTGGGGATGGAAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTT‑GCG‑CAGGAG‑ATAAGGCAAACC‑GATCGC‑AACATTGGGCGCGATCGTTCCCCTCGCCCCTCTGGGGAGAGGG‑TTAGGGTGAGGGGAAAAGGCGGCAt                                                                                                                                                           <  2:277427/210‑1 (MQ=255)                                                                         cacGCTGTTGGGCATCACGGGTTTGTGGCTGGCAGTGTTGGCAGATACGGGGGCCACGGTCCTGGTGAGAGCGAATGCCTTAAGATTGTT‑GCG‑CAGGAG‑ATAGGCCAAACC‑CATCGC‑AACATTGAGCGCGTTCGGTCCCCTCGCCCCTCTGGGGAGAGGG‑TTGGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAAcc                                                                                                                                  <  2:423168/210‑1 (MQ=255)                                                                                  aggggggccggggTTGTGGCTGGAGGGGCGGGCGGAAGGGGGGGGGACGGTGCTGGGAACGGCAATGGCTTTAAGTTTGTT‑GCG‑AAGGAG‑AAAGGGAAACCC‑GATCGA‑AACATGGGGCGGGATCGGCCCCCCCCCCCCTCGGGGGAGAGGG‑TTAGGGTGGGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACg                                                                                                                         <  2:502531/200‑1 (MQ=255)                                                                                   gggATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTT‑GCG‑CAGGAG‑ATAAGGCAAACC‑GATCGC‑AACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGG‑TTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGa                                                                                                                        >  1:81611/1‑210 (MQ=255)                                                                                   gggAGGCCCGGTTTGGGCCTGGGAGTGGTGGCAGATGGGGGGGCACCGGTGCTGGGGCAAGCGAAAGGTTTAAGATTGTG‑GCG‑CAGGGA‑ATAAGCCAACCC‑AACGGA‑AACATTGGCCGGGATGGTCCCCCTCCCCCCTCTGGGGAGGGGG‑TTAGGGTGAGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGa                                                                                                                        <  2:358303/210‑1 (MQ=255)                                                                                     gATGACGGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGTTTGTG‑GCG‑CAGGAG‑ATAAGGCAAACC‑GATCGC‑ACCATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGG‑TTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAAt                                                                                                                      <  2:51419/210‑1 (MQ=255)                                                                                       tGACGGGGTTGGGGCTGGCGGGGCGGGCAGTTACGGGGGCGACGGTGCTGGGAACAGCGAATGCGTTAAGATTGTG‑GCG‑CAGGAG‑ATAAGGCAACCC‑GATCGC‑AACATTGAGCGCGATCGGCCCCCCCGCCCCTCTGGGGAGAGGG‑TTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCa                                                                                                                    <  2:542257/210‑1 (MQ=255)                                                                                                 gtgGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTT‑GCG‑CAGGAG‑ATAAGGCAAACC‑GATAGC‑AACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGG‑TTAGGGTGAGGGGAAAAGACGGCATCGAAGCCAATCAGCCCGTATAAACCGCCTTTATGAAGCAAATAACGaga                                                                                                          >  1:151696/1‑208 (MQ=255)                                                                                                       ggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTT‑GCG‑CAGGAG‑ATAAGGCAAACC‑GATCGC‑AACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGG‑TTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGAGAAGGCAg                                                                                                    >  1:561223/1‑210 (MQ=255)                                                                                                               gggcggatgggggggggcgggTGATGGTGACAGGCAATCGGTTAAAGTGGTG‑GCA‑TGGGAG‑ATAAGGAAAACC‑GATAGA‑AAATTTGGGCGCAATGGGCCCCCTCGCCCCTCGGGGGAGAGGG‑TTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCg                                                                                            <  2:15744/191‑1 (MQ=255)                                                                                                                ggCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTT‑GCG‑CAGGAG‑ATAAGGCAAACC‑GATCGC‑AACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGG‑TTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCAATCCGt                                                                                           >  1:316251/1‑210 (MQ=255)                                                                                                                      ggcGGGGGGGACGGGGCGGGTGACAGCGAAGGCGTTAAGATTGTT‑GCG‑AAGAGG‑ATAGGGAAAACC‑GTGCGC‑ACAGTGGGGCGCGGTGGGTCCCCTCGCCCCTCTGGGGAGGGGG‑TGAGGGGGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACAAATCAAATAACGATAAGGCAGTCCATCCGTCTCttt                                                                                     <  1:213978/208‑1 (MQ=255)                                                                                                                        cGGGGGCGACGGTGCTGTTGACAGCGAATGCGTTAAGATTGTT‑GCG‑CAGGAG‑ATAAGGCAAACC‑GATCGC‑AACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGG‑TTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCGTTACGAATCAAATAACGATAAGGCAGTCCAGCCTTCTCTTTAg                                                                                   >  1:64234/1‑210 (MQ=255)                                                                                                                                 gggTGCGGGTGAGCGAGAAGGCGTTAAGAGTGTT‑GCG‑CAGGGA‑GTAAGGAAACCC‑GATGGG‑AACATGGGGCGCGGTCGGCCCCCTCGCCCCTCTGGGGAGAGGG‑TTAGGGTGAGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGtt                                                                          <  1:29160/209‑1 (MQ=255)                                                                                                                                     actgtgaaCAGCAAGTGCTTCAAGTTTG‑T‑GCGCCAGGAG‑AAACGGCAACCC‑GA‑CGCAAACTTTGAGCACGATCGGTCCCCTCGCCGTCTTGGGTAGGGGT‑CTAGGGTAGGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgtg                                                                      <  2:151709/203‑1 (MQ=255)                                                                                                                                                    ggCGTTAAGATGGTG‑GCG‑CAGGAG‑ATAAGGCAAACC‑GATCGA‑AACATTGGGGGCGATCGGCCCCCTCGCCCCTCGGGGGAGAGGG‑TTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACa                                                       <  1:48980/209‑1 (MQ=255)                                                                                                                                                            gagtgtg‑gCG‑AAGGCG‑AAAAGGCAACCC‑GATCGC‑AACATGGGCGGCGATGGGTCCCCCCGCCCCTGTGGGGAGGGGT‑TTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAg                                               <  1:574752/207‑1 (MQ=255)                                                                                                                                                            gATTGTT‑GCG‑CAGGAG‑ATAAGGCAAACC‑GATCGC‑AACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGG‑TTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAg                                               >  2:38094/1‑210 (MQ=255)                                                                                                                                                            agttgttgGCG‑CAGGAG‑AACAGGCAACCA‑GCGCGA‑AACAGTGAGAGCGATCGGCCCCCCCGCCCCTCTGGGGGGGGGG‑TAAGGGGGGGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACGAAACCCa                                                <  1:15924/208‑1 (MQ=255)                                                                                                                                                             agtggt‑gCG‑CAGGGGAAAAAGGCAAACG‑ATTCGA‑AACTTGGGGGGCGGTGGGCCCCCTCGCCCTGCTGGGGGGAGGG‑TAGGGGTGGGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAg                                               <  1:27897/205‑1 (MQ=255)                                                                                                                                                              ttgtt‑gCG‑CAGGAG‑ATAAGGCAAACC‑GATCGC‑AACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGG‑TTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGa                                             >  2:473432/1‑210 (MQ=255)                                                                                                                                                               tgtg‑gCG‑AAGGAG‑ATAGGGCAAACG‑TGTCGA‑AACATGGGGCGGGTCCGGCCCCCCCGCCCGTCGGGGGGGGGGG‑TTAGGGGGGGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGaa                                            <  1:580098/210‑1 (MQ=255)                                                                                                                                                                      g‑cAGGAG‑ATAAGGCAAACC‑GATCGC‑AACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGG‑TTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATcgc                                      >  1:137489/1‑210 (MQ=255)                                                                                                                                                                      ccaagaag‑aTAAGGAAAACC‑GCTCGA‑AAGATTGAGCGGGAGCGGCCCCTGCCCCCCTCGGGGGAGAGGG‑TTAGGGAGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATGAAGGTACAAAACCCAGGAATATcg                                       <  2:316281/204‑1 (MQ=255)                                                                                                                                                                         aGGAG‑A‑AAGGCCAACCAGGTCGC‑AACATGGGGCGCGATCGCCCCCCCCGCCCCTCGGGGGAGAGGG‑TTAGGGTGAGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCgg                                    <  2:507874/210‑1 (MQ=255)                                                                                                                                                                                 aaGGCAAACC‑GATCGC‑AACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGG‑TTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATACAGCGGGGAGTgg                             >  1:56154/1‑210 (MQ=255)                                                                                                                                                                                       aacacccgacgaaagagtGAGGGCGCTGGGCCCCCCCGCCCGTATGGGGGGGGGG‑GTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCgg                         <  1:58437/193‑1 (MQ=255)                                                                                                                                                                                         cc‑GTTCGA‑AACAGTGAGAGCGAGGGGCACGCCCGCCCCTCTGGGCGGAGGG‑TTAGGGTGAGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCACCCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGAtcg                     <  2:118382/210‑1 (MQ=255)                                                                                                                                                                                             acgcc‑AACATGGGGCGCGTTCGGCCCCCTCGCCCCTGGGGGGAGGGGG‑TTAGGGTGAGGGAAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGAtcgtcg                  <  2:258335/206‑1 (MQ=255)                                                                                                                                                                                                       ttGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGG‑TTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATc         >  2:564979/1‑210 (MQ=255)                                                                                                                                                                                                              cgATTGGGCCCCTCG‑CCCTTTGGGGGGTGGGTTTAGGGTGAGGGGCAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACg  <  2:64242/210‑1 (MQ=255)                                                                                                                                                                                                                     |                                                                                                                                                                                                            ACGCGCCACTCACGCCAATATCCGCCAGAACATCACTATTGCGCTGGGGCTGAAAGGGATCTTCCTCGTCACCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTT‑GCG‑CAGGAG‑ATAAGGCAAACC‑GATCGC‑AACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGG‑TTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGCGATAATCAGCAACG  >  NC_000913/3608475‑3608884

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.