Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,332,791	0	A	C	42.9%	121.9 / 11.2	63	V161G (GTG→GGG)	basS	sensory histidine kinase in two‑component regulatory system with BasR
Reads supporting (aligned to +/- strand):  ref base A (31/5);  new base C (0/27);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 4.76e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATGTTGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGA‑CGCTCTCCATCATCTGATCAAGCCGTG‑CCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑A‑GATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCG  >  NC_000913/4332642‑4332941                                                                                                                                                           |                                                                                                                                                       cATGTTGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGA‑CGCTCTCCATCATCTGATCAAGCCGTG‑CCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑A‑GATGCAAACGCAc                                                                                                                                                       <  2:1050937/151‑1 (MQ=255) cATGTTGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGA‑CGCTCTCCATCATCTGATCAAGCCGTG‑CCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑A‑GATGCAAACGCAc                                                                                                                                                       <  1:625185/151‑1 (MQ=255)  aTGTTGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGA‑CGCTCTCCATCATCTGATCAAGCCGTG‑CCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑A‑GATGCAAACGCAcc                                                                                                                                                      >  2:190480/1‑151 (MQ=255)      tGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGA‑CGCTCTCCATCATCTGATCAAGCCGTG‑CCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑A‑GATGCAAACGCAccccc                                                                                                                                                   <  2:286629/150‑1 (MQ=255)        aTAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGA‑CGCTCTCCATCATCTGATCAAGCCGTG‑CCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑A‑GATGCAAACGCACCCCCGcc                                                                                                                                                >  1:804433/1‑151 (MQ=255)          aaTTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGA‑CGCTCTCCATCATCTGATCAAGCCGTG‑CCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑A‑GATGCAAACGCACCCCCGcc                                                                                                                                                >  1:190533/1‑149 (MQ=255)          aaTTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCGGGGAGA‑CGCTCTCCATCATCTGATCAAGCCGTG‑CCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑A‑GATGCAAACGCACCCCCGCCa                                                                                                                                               <  1:1091011/150‑1 (MQ=255)           aTTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGA‑CGCTCTCCATCATCTGATCAAGCCGTG‑CCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑A‑GATGCAAACGCACCCCCGCCAg                                                                                                                                              >  1:165124/1‑150 (MQ=255)              aCCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGA‑CGCTCTCCATCATCTGATCAAGCCGTG‑CCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑A‑GATGCAAACGCACCCCCGCCAGTGGc                                                                                                                                          >  2:1099164/1‑151 (MQ=255)                  gaagaaAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGA‑CGCTCTCCATCATCTGATCAAGCCGTG‑CCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCT‑A‑GATGCAAACGCACCCCCGCCAGTGGCGTTc                                                                                                                                      >  1:638487/1‑151 (MQ=255)                   aagaaAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGA‑CGCTCTCCATCATCTGATCAAGCCGTG‑CCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑A‑GATGCAAACGCACCCCCGCCAGTGGCGTTCg                                                                                                                                     >  2:155766/1‑151 (MQ=255)                   aaaaaaaTGACTGCCCGGCCCGCCCCAGTTGCAGCAGCGGGGAGC‑CCCCCTCCATCATTTGATCAAGCCGTC‑CCACAAACGGTG‑CTACATCAATGTGATGCTTTTCCGCCAGCAGTTCC‑A‑GATGAAACCGCCCCCCCGCCAGTGGCGTTc                                                                                                                                      <  1:27007/147‑1 (MQ=255)                    aaaaaaTGACTGCCGGGCCCGCGCCAGTGGCAGCGGTGGGGGGC‑CCCTCCCCATCTTCTGACCAGGCCGTC‑CCACAAACGGTG‑CTACATCAATGTGATGCTTTTCCCCCAGCATTTCC‑A‑GAGGAAAACCCCCCCCCGCCAGTGGCGTTCGc                                                                                                                                    <  2:221773/149‑1 (MQ=255)                     aaaaaGGACTTCCCGGCCCCCGCCGGTTGCAGCAGCTGGGAGA‑CCCTCCCCATCATCTGATAAAGCGGTC‑CCCCTAACGGGG‑CTACATCAATGTGATGCGTTTCCGCCAGCAGTTCC‑A‑GAGCCAAACGCCCCCCCGCCAGTGGCGTTCGCa                                                                                                                                   <  1:1132451/150‑1 (MQ=255)                         atgcttGTCCGGCACGCGCCTGTTGCAGCGGCTGGGAGA‑CCCTCTCCATCATCTGACCAAGCCGTG‑CCACTAACGGTG‑CTACATAAATGTGATGCTTTTTCCCCAGCAGTCCC‑A‑GATGCAACCCCCCCCCCGCCAGTGGCGTTCGCAGTTc                                                                                                                               <  1:968606/146‑1 (MQ=255)                              tGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGA‑CGCTCTCCATCATCTGATCAAGCCGTG‑CCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑A‑GATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTgc                                                                                                                           >  2:309418/1‑150 (MQ=255)                                   ggCACGCGCCAGTTGCAGCAGCTGGGAGA‑CGCTCTCCATCATCTGATCAAGCCGTG‑CCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑A‑GATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACg                                                                                                                     >  1:601741/1‑151 (MQ=255)                                     cccgccCCAGTGGCAGCAGCGGGGGGC‑CCCTCCCCCCCACTTAACCAACCCGCG‑CCAC‑AACGGGGCCTACATCAATGTGTTGCGTTTGCGCCAGAGGTTCC‑A‑GATAAAAACCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGcgcg                                                                                                                        <  2:276281/144‑1 (MQ=255)                                      aCGCGCCAGTTGCAGCAGCTGGGAGA‑CGCTCTCCATCATCTGATCAAGCCGTG‑CCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑A‑GATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTc                                                                                                                   >  2:203766/1‑150 (MQ=255)                                        gcgcCAGTTGCAGCAGCTGGGAGA‑CGCTCTCCATCATCTGATCAAGCCGTG‑CCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑A‑GATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAg                                                                                                                 >  1:333183/1‑150 (MQ=255)                                        gcgcCAGTTGCAGCAGCTGGGAGA‑CGCTCTCCATCATCTGATCAAGCCGTG‑CCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑A‑GATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAg                                                                                                                 >  2:808968/1‑150 (MQ=255)                                         cgcCAGTTGCAGCAGCTGGGGGC‑CCCTCTCCATCATCTGATCAAGCCGTG‑CCACTAACGGTG‑CTACATCAATGTGATGCGTTTCCGCCAGCAGTTCC‑A‑GATGCAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGc                                                                                                                <  1:18503/150‑1 (MQ=255)                                         cgcCAGTTGCAGCAGCTGGGAGA‑CGCTCTCCATCATCTGATCAAGCCGTG‑CCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑A‑GATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCg                                                                                                               >  2:465730/1‑151 (MQ=255)                                           ccAGTTGCAGCAGCTGGGAGA‑CGCTCCCCATCATCTAACCAAGCCGGC‑CCACAACGGGTG‑CTACATAATTGTGATGCTTTTCCGCCAGCAGTCCC‑A‑GATGAAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAg                                                                                                                 <  1:16864/147‑1 (MQ=255)                                               ttGCAGC‑GGTGGGGGACCGCTCCCCATCATCTGATAAGGCCGTC‑CCCCAAACGGTG‑CAACACAAATGTGATTTTTTTCCGCCAGCAGTCCC‑A‑GAGGAAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTaa                                                                                                           <  2:221057/149‑1 (MQ=255)                                                    gcagcTGGGAGA‑CGCTCTCCATCATCTGATCAAGCCGTG‑CCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑A‑GATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAaacaac                                                                                                      >  1:921488/1‑149 (MQ=255)                                                       ggtGGGAGC‑CGCTCCCCATCATCTGAACAGCCCGTC‑CCACTAGGGGTG‑CTACATCAATGGGATGCGTTTTCGCCAGC‑GTTCCAA‑GAGGAAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAaaca                                                                                                        <  2:1106355/142‑1 (MQ=255)                                                       gcggggAGC‑CGCTCCCCATCACCTGATCAGGCCGTG‑CCACTAACGGTG‑CTACATCATTTGGATGTGTTTCCGCCGGCATTTCC‑A‑GATGCAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTc                                                                                                 <  2:978392/148‑1 (MQ=255)                                                        cTGGGAGA‑CGCTCTCCATCATCTGATCAAGCCGTG‑CCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑A‑GATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCg                                                                                                >  1:957649/1‑151 (MQ=255)                                                         ggggAGC‑CGCTCCCCACCATCTGACCAGGCCGGC‑CCCCTACGGGTG‑TTAAAAAAGTGGGTTGCTTTTCCCCCAGCAGTCCA‑A‑GAGGCAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTaaa                                                                                                          <  2:47279/139‑1 (MQ=255)                                                          gggAGA‑CGCTCTCCATCATCTGATCAAGCCGTG‑CCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑A‑GATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGtt                                                                                              >  1:290634/1‑151 (MQ=255)                                                           ggggtccGCTCCCCATCATCGGAGCAAGCCTTC‑CCACAACGGGTG‑CACCAAAATTGTGATGCGTTTGCGCCAGCAGTCCC‑A‑GATCCAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTc                                                                                                 <  1:1115329/142‑1 (MQ=255)                                                              ga‑CGCTCTCCATCATCTGATCAAGCCGTG‑CCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑A‑GATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATcc                                                                                          >  1:990553/1‑151 (MQ=255)                                                                 cGCTCTCCATCATCTGATCAACCCGCG‑CCACTAAGGGTG‑CAACATAAATGGGATGCTTTTTCCCCGGCAGTTAC‑A‑GATGCAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCa                                                                                         <  1:1037727/150‑1 (MQ=255)                                                                 cGCTCTCAATCTTCTGA‑CCAGGCGTGCCCATTAGGGGTG‑CAACAAAATTGTGTTGGGTTTTCGCGA‑CAATTCC‑AGGATGAAAGCGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTAt                                                                                            <  2:1028489/147‑1 (MQ=255)                                                                   ctctCCATCATCTGATCAAGCCGTG‑CCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑A‑GATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCg                                                                                      >  1:974317/1‑151 (MQ=255)                                                                   ctctCCATCATCTGATCAAGCCGTG‑CCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑A‑GATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCg                                                                                      >  1:533675/1‑151 (MQ=255)                                                                     ccccATCATCTGACAAAGCCGTG‑CCACTACCGGTC‑CTACATAAATGTGATGCTTTTTCGCCGGCAGTTCC‑A‑GAGCCAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCa                                                                                         <  2:804422/144‑1 (MQ=255)                                                                          tctcttGATAAAGCCGTC‑CCCTTAACGGTG‑CAACATCAAGGGGATGCTTTTCCGCCAGCATTTCC‑A‑AAGGCAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTg                                                                                 <  2:560266/144‑1 (MQ=255)                                                                                gATCAAGCCGTG‑CCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑A‑GATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCg                                                                         >  2:318483/1‑151 (MQ=255)                                                                                  tCAAGCCGTG‑CCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑A‑GATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGAc                                                                       >  2:640758/1‑151 (MQ=255)                                                                                        cggtcccccTACCGGTG‑CTACATAATTGTGATGCGTTTCCGCCAGCAGTTCA‑A‑GATGCAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGa                                                                    <  2:51952/141‑1 (MQ=255)                                                                                        cgatccccccgagggTG‑ATCCATCAATGTGATGC‑TTTTCCCCCCGAGTTCCGA‑AATGCAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGAc                                                                   <  2:781982/137‑1 (MQ=255)                                                                                        cGTG‑CCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑A‑GATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGAcc                                                                  >  1:235423/1‑150 (MQ=255)                                                                                             ccacaacgggTG‑CTAAATCATTGGGATGCGTTTCCGCCAGCATTCCC‑A‑GATGAAAACCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGAt                                                              <  1:190504/142‑1 (MQ=255)                                                                                                     ggTG‑TTACATCATGGTGGTGCGTTTCCCCCACCAGTCCA‑A‑GTGGAAACCCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGttt                                                         <  2:891011/147‑1 (MQ=255)                                                                                                             cATCAATGTGATGCTTTTCCCCCAGCAGTTCC‑A‑GATGCAAACCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGaaa                                               <  2:94224/150‑1 (MQ=255)                                                                                                                cAATGTGATGCGTTTTCGCCAGCAGTTCC‑A‑GATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAaccac                                           >  1:111305/1‑151 (MQ=255)                                                                                                                cAATGTGATGCGTTTTCGCCAGCAGTTCC‑A‑GATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAaccac                                           >  2:799072/1‑151 (MQ=255)                                                                                                                            ttttCGCGAGCAGTCCA‑A‑GATGCAAACCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGAtt                                 <  2:732585/149‑1 (MQ=255)                                                                                                                            tttcccccAGCATTTCC‑A‑GATGAAAACCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGAtt                                 <  1:846226/143‑1 (MQ=255)                                                                                                                            tttcCGCCAGCAGTCCC‑A‑GATGCAACCCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGAttt                                <  2:886344/150‑1 (MQ=255)                                                                                                                             tttCGCCAGCAGTTCC‑A‑GATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCg                              >  2:117078/1‑151 (MQ=255)                                                                                                                                  ccAGCAGTTCC‑A‑GATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGt                         >  2:1068763/1‑151 (MQ=255)                                                                                                                                      gctttcaaa‑GATACAAACGCCCCGCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGa                             <  2:642505/136‑1 (MQ=255)                                                                                                                                      cagTTCC‑A‑GATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGgcg                     >  1:627156/1‑151 (MQ=255)                                                                                                                                         ttCC‑A‑GATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCtg                  >  2:872122/1‑151 (MQ=255)                                                                                                                                         tccc‑A‑GATGCAACCCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGgcgc                    <  2:949686/147‑1 (MQ=255)                                                                                                                                          tCC‑A‑GATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCtgt                 >  1:220036/1‑151 (MQ=255)                                                                                                                                            c‑A‑GATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCtgtg                >  1:537991/1‑150 (MQ=255)                                                                                                                                                gATGCAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGaa              <  2:447179/150‑1 (MQ=255)                                                                                                                                                 agacAAACGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATg            <  1:155783/148‑1 (MQ=255)                                                                                                                                                     aaaCGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCa         >  1:467634/1‑150 (MQ=255)                                                                                                                                                           aCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCg  >  2:25350/1‑151 (MQ=255)                                                                                                                                                           |                                                                                                                                                       CATGTTGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGA‑CGCTCTCCATCATCTGATCAAGCCGTG‑CCACTAACGGTG‑CTACATCAATGTGATGCGTTTTCGCCAGCAGTTCC‑A‑GATGCAAACGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCG  >  NC_000913/4332642‑4332941

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