Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,378,802	0	A	C	40.0%	89.8 / 21.1	55	V48G (GTA→GGA)	ampC	penicillin‑binding protein; beta‑lactamase, intrinsically weak
Reads supporting (aligned to +/- strand):  ref base A (31/2);  new base C (0/22);  total (31/24)
Fisher's exact test for biased strand distribution p-value = 2.12e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCT‑GTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAATT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGG  >  NC_000913/4378653‑4378948                                                                                                                                                           |                                                                                                                                                    ccccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCT‑GTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAATT‑ac                                                                                                                                                    <  2:1071906/151‑1 (MQ=255)   ccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGGTTGCT‑GTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGGAAAGTAATAAGGATT‑ACCC‑TG‑ATAAATT‑acc‑g                                                                                                                                                 >  1:194898/1‑151 (MQ=255)          cAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCT‑GTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAATT‑acc‑gccaccg                                                                                                                                           <  2:806276/150‑1 (MQ=255)           aaTAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCT‑GTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAATT‑acc‑gccaccgcc                                                                                                                                         >  2:1020274/1‑151 (MQ=255)           aaTAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCT‑GTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAATT‑acc‑gccaccgc                                                                                                                                          <  1:1059550/150‑1 (MQ=255)              aGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCT‑GTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAATT‑ACC‑GCCACCGCCAt                                                                                                                                       >  2:632562/1‑150 (MQ=255)                     ccaccaAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCT‑GTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAATT‑ACC‑GCCACCGCCATACCCggg                                                                                                                                >  2:330743/1‑150 (MQ=255)                         caAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCT‑GTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAATT‑ACC‑GCCACCGCCATACCCGGGATCtt                                                                                                                           >  1:138221/1‑151 (MQ=255)                             cACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCT‑GTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAATT‑ACC‑GCCACCGCCATACCCGGGATCTTttgt                                                                                                                       >  2:659104/1‑151 (MQ=255)                              aCGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCT‑GTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAATT‑ACC‑GCCACCGCCATACCCGGGATCTTttgtt                                                                                                                      >  2:386262/1‑151 (MQ=255)                               cgccataaatttttttGCGGCCCAACCCTAACTAAAACACCTTTTGTT‑GGGTGGGGGGCTTTTTTTGGGCGATGCCCGCATCCCCCAAGGAAAAGTAAAAGGTTTT‑CCCC‑TA‑AAAAATC‑CCC‑GCCACCGCCATACCCGGGATCTTttgtt                                                                                                                      <  2:253710/138‑1 (MQ=255)                                  cAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCT‑GTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAATT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGctc                                                                                                                  >  2:713271/1‑151 (MQ=255)                                  cAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCT‑GTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGGTT‑ACCC‑TG‑ATAAATT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGctc                                                                                                                  >  1:1105001/1‑151 (MQ=255)                                     tAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCT‑GTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAATT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCtat                                                                                                               >  1:17492/1‑151 (MQ=255)                                     tAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCT‑GTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAATT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCtat                                                                                                               >  1:58593/1‑151 (MQ=255)                                      aatttttttGCGGCCCAACCCTAAC‑AAAAAAAAGTTTGGTGGTGTGCGGGGCTGTTTTTGGGGGATGCCCGCATACCCCAAGGAAAAGAAATAAGTTTT‑CCCC‑TA‑AAAAATC‑CCC‑CCCACCGCCATACCCGGGATCTTTTGTTGctc                                                                                                                  <  2:784550/142‑1 (MQ=255)                                         ttttttGCTGCCCGAACCTAACTCAACCACTGTTTGCT‑GGGGGGGGGTCTTTTTTTTGGGGGTGCCCGAAAAGCCCGAGGAAAAATAAT‑AGGTTT‑CCCCTTA‑AAAAATT‑CCC‑GCTACCGCCATACCCGGGATCTTTTGTTGCTCTATAAg                                                                                                            <  2:1060794/148‑1 (MQ=255)                                          gTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCT‑GTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAATT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCg                                                                                                          >  1:1166789/1‑151 (MQ=255)                                                cTGACCGAACCTAACTCAAACAACGTTTGCT‑GTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAATT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTa                                                                                                     >  1:169259/1‑150 (MQ=255)                                                  gACCGAACCTAACTCAAACAACGTTTGCT‑GTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAATT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAAtt                                                                                                  >  1:410258/1‑151 (MQ=255)                                                   aCCGAACCTAACTCAAACAACGTTTGCT‑GTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAATT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATtg                                                                                                 >  2:539314/1‑151 (MQ=255)                                                    caaacccaaccaaaaaCACTGTTTGTT‑TTGTGAGGGTTTGTTTTTGGGCGATGCCCGCA‑GCCCCCAGGGAAAAAAAATATGTTT‑TCCC‑TT‑AAAAAAT‑CCCGGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAAt                                                                                                   <  1:415458/135‑1 (MQ=255)                                                        aCCTAACTCAAACAACGTTTGCT‑GTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAATT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCg                                                                                             >  1:922986/1‑150 (MQ=255)                                                            accaaaaaCAACGTTGGTA‑GGGTGACGGGTTGTTTTTTGGCGATGCCCGAATACCCCTAAGAAAAATAA‑AAGGTTT‑TCCC‑CGAAAAAATC‑ACC‑CCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGca                                                                                        <  2:428818/146‑1 (MQ=255)                                                            aaCTCAAACAACGTTTGCT‑GTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAATT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGc                                                                                         >  2:333432/1‑150 (MQ=255)                                                            aaCTCAAACAACGTTTGCT‑GTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAATT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGc                                                                                         >  2:1130065/1‑150 (MQ=255)                                                                cAAACAACGTTTGCT‑GTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAATT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAat                                                                                    >  1:144463/1‑151 (MQ=255)                                                                           ttgttttgtgTCGTGCTGTTTTTTGCCGATTTCCGCATGGCCCCAGGAAAAGTAATAAGTTTT‑ACCC‑TA‑ATAAA‑T‑CCC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAtt                                                                          <  1:550677/144‑1 (MQ=255)                                                                             cT‑GTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAATT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgt                                                                       >  2:528008/1‑151 (MQ=255)                                                                                 tgggcggggtttttttttGGCGTTGTCCACACACCCCCAGAAAAATAAAAAAGTTTT‑CCCT‑TA‑AAAAATT‑CCC‑CCCTCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgt                                                                       <  1:371485/135‑1 (MQ=255)                                                                                    gACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAATT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAggg                                                                 >  1:620731/1‑151 (MQ=255)                                                                                      gggggtttttttgGGCGATTCCCGCACA‑CCCCAGGTAAAAAAAAAAGGTTTCTCCT‑TT‑AAAAATC‑ACC‑GCCGCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgc                                                               <  1:258589/143‑1 (MQ=255)                                                                                         gctgctTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAATT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgcagc                                                            >  2:484817/1‑151 (MQ=255)                                                                                           tttttttgGGCGTTGCCCGCACAGCCCGAGGAAAAGAAATAAGTTTT‑CCCT‑TA‑AAAAATT‑CCC‑GCTTCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgcagc                                                            <  1:805402/146‑1 (MQ=255)                                                                                            ttttttgGGCATTGTCCAAATACCCCCAGAAAAAATAAAAGGTTTT‑CCCC‑AA‑ATAATTT‑CCC‑CCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgcagca                                                           <  2:804794/147‑1 (MQ=255)                                                                                               ttttGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACTC‑TG‑ATAAATT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAAtgt                                                      >  2:1112710/1‑151 (MQ=255)                                                                                                       aTGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAATT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCa                                                >  1:695055/1‑149 (MQ=255)                                                                                                         gTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAATT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGTTGCACAATATCGTTGATTTGGTGAGGGGCAGCAAATGTGGAGCAag                                              >  2:964814/1‑149 (MQ=255)                                                                                                          tCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAATT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAgg                                           >  2:1226174/1‑151 (MQ=255)                                                                                                          tCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAATT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAgg                                           >  2:678074/1‑151 (MQ=255)                                                                                                           ccGCATAGCCCCAGGAAAAGAAATAAGTTTT‑CCCT‑TA‑AAAAATT‑CCC‑CCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAagag                                            <  2:882487/149‑1 (MQ=255)                                                                                                             gCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAATT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCgg                                        >  2:489259/1‑151 (MQ=255)                                                                                                                   ccccgagaaaaaaaataaGTTTT‑TCCC‑TG‑AAATATA‑CCC‑CCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaat                                    <  1:591335/136‑1 (MQ=255)                                                                                                                   cccccagaaaaaTAAAAAGGTTT‑TCCC‑TT‑AAAAATTCCCC‑CCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaatt                                   <  1:1109612/143‑1 (MQ=255)                                                                                                                   ccccAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAATT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaatt                                   >  2:851596/1‑150 (MQ=255)                                                                                                                   ccccAGGAAAGAAAATAGGTTTT‑ACTA‑TA‑ATAATTT‑CCC‑GCCCCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaatta                                  <  2:1046976/151‑1 (MQ=255)                                                                                                                    cccaagaaaaaataaaaGGTTT‑TCCCTTG‑AAAATTC‑CCC‑CCCTCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaatta                                  <  1:141990/136‑1 (MQ=255)                                                                                                                    ccaagtaaaaGTAATAAGTTTT‑CCCT‑TA‑AAAATTT‑CCC‑CCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaattaa                                 <  2:1152689/144‑1 (MQ=255)                                                                                                                           aaattaatagttttt‑CCCC‑TA‑AAAAATT‑CCC‑CCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaattaat                                <  1:459338/133‑1 (MQ=255)                                                                                                                           aaaGAAATAAGTTTT‑CCCC‑TG‑ATAATTT‑CCC‑CCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATtaata                               <  2:50330/146‑1 (MQ=255)                                                                                                                                  ataGGTTT‑ACCCTTG‑AAAATTT‑CCC‑CCCTCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCa                        <  2:410270/145‑1 (MQ=255)                                                                                                                                        tt‑ACCT‑TA‑ATAAATC‑CCC‑CCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGcgtcgt               <  1:7654/149‑1 (MQ=255)                                                                                                                                         ttttcccttaaTAATTT‑CCC‑CCCTCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCa                        <  1:812059/132‑1 (MQ=255)                                                                                                                                         ttcccc‑TA‑ATAATTA‑CCC‑GCACCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGc                    <  1:1235158/141‑1 (MQ=255)                                                                                                                                           aCCC‑TG‑ATAAATT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGtttt            >  1:702592/1‑150 (MQ=255)                                                                                                                                              ccttaaTAAATC‑CCC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGtt              <  2:1088267/142‑1 (MQ=255)                                                                                                                                              c‑TG‑ATAAATT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGaa         >  1:917428/1‑150 (MQ=255)                                                                                                                                                     aaaTT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACa       >  1:189834/1‑147 (MQ=255)                                                                                                                                                      aaTT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAggg  >  2:24900/1‑151 (MQ=255)                                                                                                                                                       aTT‑CCC‑CCCCCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAgg   <  2:750167/149‑1 (MQ=255)                                                                                                                                                           |                                                                                                                                                    CCCCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCT‑GTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAATT‑ACC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGG  >  NC_000913/4378653‑4378948

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: ‑