Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,378,802	0	A	C	54.4%	57.8 / 89.3	58	V48G (GTA→GGA)	ampC	penicillin‑binding protein; beta‑lactamase, intrinsically weak
Reads supporting (aligned to +/- strand):  ref base A (25/1);  new base C (0/31);  total (25/33)
Fisher's exact test for biased strand distribution p-value = 2.62e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.79e-10

CCCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGT‑CCGCATAG‑CCCCAGGT‑AAAGTAATAAGGTTTACCC‑TGATAAA‑TTACCG‑C‑CACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGG  >  NC_000913/4378654‑4378948                                                                                                                                                          |                                                                                                                                                     cccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGT‑CCGCATAG‑CCCCAGGT‑AAAGTAATAAGGTTTACCC‑TGATAAA‑TTaca                                                                                                                                                    <  2:587299/151‑2 (MQ=255)    tCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGT‑CCGCATAG‑CCCCAGGT‑AAAGTAATAAGGTTTACCC‑TGATAAA‑TTaccg‑c‑c                                                                                                                                               >  1:777981/1‑151 (MQ=255)       aGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGT‑CCGCATAG‑CCCCAGGT‑AAAGTAATAAGGTTTACCC‑TGATAAA‑TTaccg‑c‑cacc                                                                                                                                            >  2:839758/1‑151 (MQ=255)                gTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGT‑CCGCATAG‑CCCCAGGT‑AAAGTAATAAGGTTTACCC‑TGATAAA‑TTACCG‑C‑CACCGCCATAccc                                                                                                                                   >  2:887436/1‑151 (MQ=255)                gTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGT‑CCGCATAG‑CCCCAGGT‑AAAGTAATAAGGTTTACCC‑TGATAAA‑TTACCG‑C‑CACCGCCATAccc                                                                                                                                   >  2:667555/1‑151 (MQ=255)                  cGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGT‑CCGCATAG‑CCCCAGGT‑AAAGTAATAAGGTTTACCC‑TGATAAA‑TTACCG‑C‑CACCGCCATACCCgg                                                                                                                                 >  2:246918/1‑151 (MQ=255)                     caccaAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGT‑CCGCATAG‑CCCCAGGT‑AAAGTAATAAGGTTTACCC‑TGATAAA‑TTACCG‑C‑CACCGCCATACCCGGGAt                                                                                                                              >  1:460943/1‑151 (MQ=255)                       caacgCACGCAAGTAATTTTTTTGCTCCCCAACCCAAACAAAAACAACTTTTGCTGTGGGGGGGGCTGTTTTTGGGCGATGC‑CCGCATCC‑CCCGAGGA‑AAAGAAAT‑AGGTTTTCCC‑TTAAAAATTTCCCC‑C‑CACCGCCATACCCGGGATc                                                                                                                             <  1:582821/146‑1 (MQ=255)                             aCGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGT‑CCGCATAG‑CCCCAGGT‑AAAGTAATAAGGTTTACCC‑TGATAAA‑TTACCG‑C‑CACCGCCATACCCGGGATCTTttgtt                                                                                                                      >  1:500752/1‑151 (MQ=255)                                 cAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGT‑CCGCATAG‑CCCCAGGT‑AAAGTAATAAGGTTTACCC‑TGATAAA‑TTACCG‑C‑CACCGCCATACCCGGGATCTTTTGTTGct                                                                                                                   >  2:395051/1‑150 (MQ=255)                                 aataaatttttttGCGGACCAAACCTAACAAAAACAACTTTTGCTGGGGGCGGGGTTGTTTTTGGGGGATGC‑CCGCATAC‑CCCAAGTA‑AAAGAAATAGGTTTTCCCT‑TGAAAAA‑TTCCCG‑C‑CACCGCCATACCCGGGATCTTTTGTTGc                                                                                                                    <  2:92965/140‑1 (MQ=255)                                     aattGTTTTGCTGACCGACCCAAACTAAAAAATGTTTTGGGGGGTGAGGGGCTGTTTTTGGGCGTGCC‑CCCCCTAC‑CCCCGGAA‑AAAGTAAGAAGTTTTCCTC‑TAAAAAT‑TTCCCA‑C‑CACCGCCATACCCGGGATCTTTTGTTGctct                                                                                                                 <  1:859226/145‑1 (MQ=255)                                     aattGTTTTGCGAACCAACCAAAACACAAACAAGTTTTGCTGTGTGGGGGGTTTTTTTTTGGTATTGC‑CCCCAAAC‑CCCAAGGA‑AAAGTAATAATTTTTCCCC‑TGATAAA‑TCCCCC‑C‑CACCGCCATACCCGGGATCTTTTGTTGctct                                                                                                                 <  1:519817/145‑1 (MQ=255)                                      aaTTTTTTGCTGCCCAAACCAACCAAAAACAACGTTTGGTGTGGGAGGGGCTGTTTTTTGGC‑TTGTCCCCCATAC‑CCCCGGGA‑AAAGTAAAAATTTTTCCCC‑TGATAAT‑TTCCCG‑C‑CACCGCCATACCCGGGATCTTTTGTTGCTCTATaa                                                                                                             <  1:528072/151‑1 (MQ=255)                                         gTTTTGCTGACCGAACCTAACTCACATAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGT‑CCGCATCG‑TCCCAGGT‑AAAGTAGTAAGGTTTACCC‑TGATAAA‑TTACCG‑C‑CATCGCCATACCCGGGCTCTTTTGTTGctct                                                                                                                 >  1:721517/1‑144 (MQ=255)                                         gTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGT‑CCGCATAG‑CCCCAGGT‑AAAGTAATAAGGTTTACCC‑TGATAAA‑TTACCG‑C‑CACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCg                                                                                                          >  1:645203/1‑151 (MQ=255)                                                tGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGT‑CCGCATAG‑CCCCAGGT‑AAAGTAATAAGGTTTACCC‑TGATAAA‑TTACCG‑C‑CACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAAt                                                                                                   >  1:49976/1‑151 (MQ=255)                                                      aaCCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGT‑CCGCATAG‑CCCCAGGT‑AAAGTAATAAGGTTTACCC‑TGATAAA‑TTACCG‑C‑CACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCg                                                                                             >  1:674490/1‑151 (MQ=255)                                                      aaCCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGT‑CCGCATAG‑CCCCAGGT‑AAAGTAATAAGGTTTACCC‑TGATAAA‑TTACCG‑C‑CACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCg                                                                                             >  2:842760/1‑151 (MQ=255)                                                        ccaaccTCAACCACCGTTTGCTGTGTGAGGGGCTGTTTTTGGGCGATGC‑CCGCATGC‑CCCCAGGA‑AAAGTAATAAGGTTTCCCC‑TGATAAT‑TTCCCG‑C‑CACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGAt                                                                                           <  2:85991/146‑1 (MQ=255)                                                              tCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGT‑CCGCATAG‑CCCCAGGT‑AAAGTAATAAGGTTTACCC‑TGATAAA‑TTACCG‑C‑CACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGcaca                                                                                      >  1:1378/1‑150 (MQ=255)                                                               cAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGT‑CCGCATAG‑CCCCAGGT‑AAAGTAATAAGGTTTACCC‑TGATAAA‑TTACCG‑C‑CACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACaa                                                                                     >  1:758444/1‑150 (MQ=255)                                                                aaaaaaCGTTTGGTGTGTGGGGGTCTGCTTTTTGGCGAGCC‑CCGCAAAG‑CCCAAGGT‑AAAGAAATAAGTTTTACCT‑AGATAAT‑TCCCCG‑C‑CACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAata                                                                                   <  1:790773/151‑1 (MQ=255)                                                                aaaCAGCTTTGGCTGTGTGACGGGTTTTTTTTTGGCGATGT‑CCGCATAG‑CCCCGGAA‑AAAGAAAAAGGGTTTCCCT‑TGAAAAA‑TTCCCC‑C‑CACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAat                                                                                    <  2:632935/150‑1 (MQ=255)                                                                aaaCAACTTTTGCTTTGAGCCGGCGTGTTTTTTGGAGTTTT‑CC‑CAAAGCCCCCAGGA‑AAAAAAAAAAGGTTTACGC‑TAAAAAT‑TGCCCG‑C‑CGCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAata                                                                                   <  2:518203/151‑1 (MQ=255)                                                                       gTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGT‑CCGCATAG‑CCCCAGGT‑AAAGTAATAAGGTTTACCC‑TGATAAA‑TTACCG‑C‑CACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTg                                                                             >  1:573983/1‑150 (MQ=255)                                                                         ttgttgTGTGAGGGGCTGTTTTTTGGGGATGT‑CCGCATAC‑CCCCGGGA‑AAAGTAATAGGTTTTCCCC‑TAAAAAT‑TTCCCG‑C‑CACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAt                                                                           <  2:710539/150‑1 (MQ=255)                                                                           gttgtgtgAGGGGCTTGTTTTGGGCGATGC‑CCGCATAC‑CCGCAGAA‑AAAAAAAAAAGGTTTCCCC‑TGAAAAA‑TTCCCC‑C‑CACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttg                                                                        <  1:646897/149‑1 (MQ=255)                                                                           gCTGTGTGCGGGGCTGTTTTTTGGCGATGC‑CCGCATAG‑CCCCAGGT‑AAAGTAATAAGTTTTCCCC‑TGATAAA‑TTCCCG‑C‑CACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAtt                                                                          <  2:760696/149‑1 (MQ=255)                                                                              gtgtgAGGGGTTGTTTTTTGGCGTTGC‑CCCCATAC‑CCCCAGAA‑AAAAAAAAAAGTTTTACCC‑TAATAAT‑TTCCCC‑C‑CACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgttg                                                                     <  2:263869/151‑1 (MQ=255)                                                                               tgtgCCGGTCTTTTTTTTGGCGATGT‑CCGCATCC‑CCCCGGGA‑AAAG‑AAAAAGGTTTTCCCTTGATAAA‑TTCCCG‑C‑CACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgttg                                                                     <  2:774030/150‑1 (MQ=255)                                                                                gggaggggCTTTTTTTGGGCGATGT‑CCGCATCC‑CCCCAGGA‑AAAGTAAAAAGTTTTACCT‑TAAAAAT‑TTCCCG‑C‑CACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAt                                                                           <  2:603367/138‑1 (MQ=255)                                                                                  ggcGGGCTGTTTTTTGGCGATGT‑CCGCATAG‑CCCCAGGT‑AAAGTAATAAGTTTTCCCT‑TAAAAAA‑TACCCG‑C‑CACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAggg                                                                 <  1:749789/149‑1 (MQ=255)                                                                                       gctgctTTTTGGCGATGT‑CCGCATAG‑CCCCAGGT‑AAAGTAATAAGGTTTACCC‑TGATAAA‑TTACCG‑C‑CACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgcagc                                                            >  2:414884/1‑151 (MQ=255)                                                                                         tgttttttGGCGTTGC‑CCGCATAC‑CCCCGGGA‑AAAGTAATAAGTTTTCCCC‑TAAAAAA‑TTCCCG‑C‑CACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgcag                                                             <  2:190895/145‑1 (MQ=255)                                                                                            tttttGGCGATGT‑CCGCATAG‑CCCCAGGT‑AAAGTAATAAGGTTTACCC‑TGATAAA‑TTACCG‑C‑CACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAAtg                                                       >  1:284311/1‑151 (MQ=255)                                                                                             ttttGGCGAGTT‑CCACTTAG‑CCCATGTA‑AAAGTAATAA‑GTTTTCCC‑CGTTAAA‑TTTCCC‑CAAACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAAt                                                        <  2:396949/149‑1 (MQ=255)                                                                                               ttGGCGATGT‑CCGCATAG‑CCCCAGGT‑AAAGTAATAAGGTTTACCC‑TGATAAA‑TTACCG‑C‑CACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTgg                                                    >  1:347033/1‑151 (MQ=255)                                                                                                 ggggATGC‑CCGAATGC‑CCCCAGGA‑AAAGAAAAAGGTTTTCCCC‑TAATAAT‑TCCCCC‑C‑CACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGa                                                   <  1:803404/147‑1 (MQ=255)                                                                                                     aTGT‑CCGCATAG‑CCCCAGGT‑AAAGTAATAAGGTTTACCC‑TGATAAA‑TTACCG‑C‑CACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCa                                                >  2:146970/1‑149 (MQ=255)                                                                                                      tGT‑CCGCATCC‑CCCCAGGT‑AAAGTAATAAGGTTTCCCC‑TGATAAA‑TTCCCG‑C‑CACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAag                                              <  2:53719/150‑1 (MQ=255)                                                                                                       tgcccGCATAC‑CCCCAGGT‑AAAGTAATAAGTTTTCCCC‑TGATAAA‑TTCCCG‑C‑CACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCaa                                               <  2:179292/146‑1 (MQ=255)                                                                                                       gtgccccgcaacccccAGGAAAAAGTAA‑AATGTTTTCCC‑TTAAAAA‑TTCCCGCC‑CACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCaa                                               <  1:296986/139‑1 (MQ=255)                                                                                                        t‑CCGCATAG‑CCCCAGGT‑AAAGTAATAAGGTTTACCC‑TGATAAA‑TTACCG‑C‑CACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAgg                                           >  2:55085/1‑151 (MQ=255)                                                                                                         cccccATAC‑CCCCGGGA‑AAAGAAATAAGTTTTCCCC‑TAAAAAT‑TTCCCC‑C‑CACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAgg                                           <  1:617711/147‑1 (MQ=255)                                                                                                          ccGCATAC‑CCCAAGGA‑AAATAAATAAGTTTTCCCC‑TGAAAAA‑TTCCCC‑C‑CACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAaga                                             <  2:689948/148‑1 (MQ=255)                                                                                                           cGCACAC‑CCCCAGGA‑AAATTAATATGGTTTCCCT‑TTATAAT‑TTCCCG‑C‑CACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAag                                              <  1:273391/146‑1 (MQ=255)                                                                                                                   ccccAGTT‑AAAGTATTAGGTTTTCCCT‑TGATAAA‑TTACCG‑C‑CACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaatt                                   <  2:344355/150‑1 (MQ=255)                                                                                                                            aaataaataaGGTTTCCCC‑TGATAAA‑TCCCCG‑C‑CACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATtaata                               <  1:149757/141‑1 (MQ=255)                                                                                                                            aaaGTAATAAGTTTTCCCT‑TAATAAA‑TTCCCG‑C‑CACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAg                             <  2:754132/148‑1 (MQ=255)                                                                                                                             aaGTAATAAGGTTTACCC‑TGATAAA‑TTACCG‑C‑CACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcg                          >  1:601027/1‑150 (MQ=255)                                                                                                                                 aataaGTTTTCCCC‑TAAAAAA‑TTCCCC‑C‑CACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCa                        <  2:389155/148‑1 (MQ=255)                                                                                                                                 aataaGTTTTACCT‑GGATAAA‑TTCCCG‑C‑CACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCa                        <  2:49978/148‑1 (MQ=255)                                                                                                                                 aataaGGTTTACCC‑TGATAAA‑TTACCG‑C‑CACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCag                       >  2:58648/1‑149 (MQ=255)                                                                                                                                      ggTTTACCC‑TGATAAA‑TTACCG‑C‑CACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGcgtcg                >  1:362416/1‑151 (MQ=255)                                                                                                                                        tttCCCC‑TGAAAAT‑TTCCCG‑C‑CACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGcgtc                 <  2:173738/148‑1 (MQ=255)                                                                                                                                           aCCT‑TGTTAAA‑TTCCCC‑C‑CACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGtttt            <  2:30518/150‑1 (MQ=255)                                                                                                                                           aCCC‑TGATAAA‑TTACCG‑C‑CACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGCTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGCTTTg           >  1:617600/1‑151 (MQ=255)                                                                                                                                                    aaa‑TTCCCG‑C‑CACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAgg   <  2:403475/151‑1 (MQ=255)                                                                                                                                                     aa‑TTACCG‑C‑CACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAggg  >  1:415425/1‑151 (MQ=255)                                                                                                                                                      a‑TTACCG‑C‑CACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAggg  >  2:899858/1‑150 (MQ=255)                                                                                                                                                          |                                                                                                                                                     CCCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGT‑CCGCATAG‑CCCCAGGT‑AAAGTAATAAGGTTTACCC‑TGATAAA‑TTACCG‑C‑CACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGG  >  NC_000913/4378654‑4378948

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