BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP016816	203,912	0	G	T	39.2%	944.2 / 267.0	839	intergenic (+296/+162)	JCVSYN2_00845/JCVSYN2_00850	lipoprotein/inorganic diphosphatase
Reads supporting (aligned to +/- strand): ref base G (147/346); new base T (318/0); total (493/346)
Fisher's exact test for biased strand distribution p-value = 0.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 892 total aligned reads displayed.

AAAAAATAAAAATTTTTCATTTAAAGCCATTTGTAAAAAAGCCTATGTTCTAGGCTTTTTATTATACTTATCAACATTACAACGCTTTCAAAG‑TGTATTTTTTAG‑G‑‑‑C‑‑‑‑TTG‑‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑GTAAAAATAAAAGTAAGATTTGGATCAACCAAAACT‑ATAT‑CAAAAAAATTAAAAAAAATAGTTTAGTTGAGTTTCAAATCTGCATATATAGTACAAAGGAATTTATGATAGAAAAAAATGGATTAGTT  >  CP016816/203786‑204039

 aaaaataaaaatTTTTCATTTAAAGCCATTTGTAAAAAAGCCTATGTTCTAGGCTTTTTATTATACTTATCAACATTACAACGCTTTCAAAG‑TGTATTTTTTAG‑G‑‑‑A‑‑‑‑TTG‑‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑GTa                                                                                                                                 >  2:755604/1‑128 (MQ=255)

      taaaaatTTTTAATTTAAATCTATTTTTAAAATAGCCTATTTTTTATTCTTTTTATTATAATTTTAAAATTTAAAAATCTTTCAAAT‑TTTTTTTTTTTT‑G‑‑‑T‑‑‑‑TTG‑‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑GTAAAAAt                                                                                                                            <  1:755604/128‑1 (MQ=255)

       ataatttttttATTTTAATCCATTTTTAAAAAATCAAATTTTATATTGTTTTTTTTATAAATTTTAAAATTTAAAAGCATTAAAAT‑TTTTATTTTTAT‑G‑‑‑T‑‑‑‑TTT‑‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑GTAAAAATa                                                                                                                           <  2:3164253/123‑1 (MQ=255)

        aaaatTTTTTTTTTTAATTCTTTTTTTAAAAATCCTTTTTTTTAGGATTTTTTTTATAATAATAAACATTAAAACGCTTTAAAAGTTTTTTTTTTTGG‑G‑‑‑T‑‑‑‑TTT‑‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑GTAAAAATa                                                                                                                           <  1:289622/128‑1 (MQ=255)

           tttttttATTTAAATACATTTTTAATAAAGCTTTTTTTTTTTGTTTTTTTTTTTATTTTTTTAAATTTAATAGTTTTAAAAT‑TGTTTTTTTTTG‑G‑‑‑CT‑‑‑TTT‑‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑GTAAAAATaaaa                                                                                                                        <  2:1853218/127‑1 (MQ=255)

            tttttCATTTTAAGCCATTTGTAAAAAAGCCTATGTTCTAGGCTTTTTATTATACTTATCAACATTACAACGCTTTCAAAG‑TGTATTTTTTAG‑G‑‑‑C‑‑‑‑TTG‑‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑TTAAATATAAAAGt                                                                                                                      >  1:216696/1‑128 (MQ=255)

            tttttCATTTAAAGCCATTTGTAAAAAAGCCTATGTTCTAGGCTTTTTATTATACTTATCAACATTACAACGCTTTCAAAG‑TGTATTTTTTAG‑G‑‑‑C‑‑‑‑TTGT‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑TTATAAATAAAAg                                                                                                                       >  2:4275523/1‑128 (MQ=255)

             ttttTATTTAATATCATTTTTAAAAAATCATATTTTCTATGTTTTTTTTTAAAATTATAAAAAATAAAAATCTTTTAAAT‑TTTTTTTTTTAGTGT‑‑T‑‑‑‑TTT‑‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑GTAAAAATAAAAg                                                                                                                       <  1:2257799/128‑1 (MQ=255)

             ttttCATTTAAAGCCATTTGTAAAAAAGCCTATGTTCTAGGCTTTTTATTATACTTATCAACATTACAACGCTTTCAAAG‑TGTATTTTTTAG‑G‑‑‑C‑‑‑‑TTGTT‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑TTAAAAATaaaat                                                                                                                       >  1:4419628/1‑127 (MQ=255)

             ttttCATTTAAAGCCATTTGTAAAAAAGCCTATGTTCTAGGCTTTTTATTATACTTATCAACATTACAACGCTTTCAAAG‑TGTATTTTTTAG‑G‑‑‑C‑‑‑‑TTGTT‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑GTAAAAATAAAAg                                                                                                                       >  1:3748959/1‑128 (MQ=255)

             ttttCATTTAAAGCCATTTGTAAAAAAGCCTATGTTCTAGGCTTTTTATTATACTTATCAACATTACAACGCTTTCAAAG‑TGTAGTTTTTAG‑G‑‑‑C‑‑‑‑TTGT‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑TTAAACATaaactt                                                                                                                      >  1:2334480/1‑125 (MQ=255)

                taaTTTATAGCCTTTTTTAAAAAACCCTTTGTTTTAGTTTTTTTTTTATTATTATAAAAATTAAAACGCTTTAAATG‑TTTATTTTTTAG‑G‑‑‑G‑‑‑‑TTGT‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑GTAAAAATAAAAGTAAg                                                                                                                   <  1:3672239/126‑1 (MQ=255)

                      aaatcaatatGTAAAAAAGCCTTTTTTTAAGGATTTTTATTAAATTATTAAAAATTACAACGCTTACAAAT‑TTTATTTTTTAG‑C‑‑‑C‑‑‑‑TTT‑‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑GTAAAAATAAAAGTAAGATTTGGa                                                                                                            <  2:3713419/119‑1 (MQ=255)

                       aaGCCATTTGTAAAAAAGCCTATGTTGTAGGCTTTTTATTATACTTATCAACATTACAACGCTTTCAAAG‑TGTATTTTTTAG‑G‑‑‑C‑‑‑‑TTGT‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑TTAAAAATaaaatatatatttgta                                                                                                            >  2:2086404/1‑116 (MQ=255)

                       aaGCCATTTGTAAAAAAGCCTATGTTCTAGGCTTTTTATTATACTTATCAACATTACAACGCTTTCAAAG‑TGTATTTTTTAG‑G‑‑‑C‑‑‑‑TTGTT‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑GTAAAAATAAAAGTAATATTTgg                                                                                                             >  2:3530500/1‑128 (MQ=255)

                         accATTTGTAAAAAACCATATTTTCATATTATTTTATTAAACTTATAAAAATAAAAACTCTTTCTAAT‑TGTTTTTTTTTT‑T‑‑‑C‑‑‑‑TTT‑‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑GTAAAAATAAAAGTAAGATTTGGATCa                                                                                                         <  1:1101355/127‑1 (MQ=255)

                                 tAAAAAATCTTTTTTTCTATGTTTTTTATTATAATTATCAAAATTAAAAATCTTTAAATG‑TGTTTTTTTTAG‑G‑‑‑C‑‑‑‑TTG‑‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑GTAAAAATAAAAGTAAGATTTGGATCAACCAAAAc                                                                                                 <  1:1498264/128‑1 (MQ=255)

                                 tAAAAAATCATATTTTATATGCTTTTTTTTATACTTATTAAAATTAAAAATCTTTAAAAT‑TTTATTTTTTAG‑G‑‑‑C‑‑‑‑TTGT‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑GTAAAAATAAAAGTAAGATTTGGATCAACCaaaa                                                                                                  <  2:4577554/128‑1 (MQ=255)

                                 tAAAAAAGCCTATGTTCTAGGCTTTTTATTATACTTATCAACATTACAACGCTTTCAAAG‑TGTATTTTTTAG‑G‑‑‑C‑‑‑‑TTGTT‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑TTTAAAATAAAAGTAATATTTGTTTAAACCaaa                                                                                                   >  1:4037994/1‑128 (MQ=255)

                                 tAAAAAAGCCTATGTTCTAGGCTTTTTATTATACTTATCAACATTACAACGCTTTCAAAG‑TGTATTTTTTAG‑G‑‑‑C‑‑‑‑TTGT‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑GTAAAAATAAAATATATATTTGGATAAACaaata                                                                                                  >  2:3625887/1‑126 (MQ=255)

                                 tAAAAAAGCCTATGTTCTAGGCTTTTTATTATACTTATCAACATTACAACGCTTTCAAAG‑TGTATTTTTTAG‑G‑‑‑C‑‑‑‑TTG‑‑‑‑‑‑‑TTTTTTTTTTTTTATTT‑GTATAAATAAAATTAATATTTGGTTCAACCaaaa                                                                                                  >  2:642285/1‑128 (MQ=255)

                                  aaaaaataatatgttttatgtttTTTATTATTCTTATAAAAATTAAAAATTTTTAAAAT‑TTTATTTTTTAT‑T‑‑‑G‑‑‑‑TTGT‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑GTAAAAATAAAAGTAAGATTTGGATCAACCAAAAc                                                                                                 <  2:1272094/119‑1 (MQ=255)

                                    aaaaGCCTATGTTCTAGGCTTTTTATTATACTTATCAACATTACAACGCTTTCAAAG‑TGTATTTTTTAG‑G‑‑‑C‑‑‑‑TTG‑‑‑‑‑‑‑TTTTTTTGTTTTT‑TTT‑TTAAAATAAAAAGTAAGATTTGGAAAAAACAAAACt‑at                                                                                             >  1:3458591/1‑128 (MQ=255)

                                     aaaGCCTATGTTCTAGGCTTTTTATTATACTTATCAACATTACAACGCTTTCAAAG‑TGTATTTTTTAG‑G‑‑‑C‑‑‑‑TTGTT‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑GTAAAAAAAAAAGTAAAATTTGGATAAAACAAAAAt‑a                                                                                              >  2:3288348/1‑128 (MQ=255)

                                     aaaGCCTATGTTCTAGGCTTTTTATTATACTTATCAACATTACAACGCTTTCAAAG‑TGTATTTTTTAG‑G‑‑‑C‑‑‑‑TTGT‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑GTAAAAATAAAAGTAAGATTTTGAACAACAAAAACttt                                                                                              >  2:3436544/1‑126 (MQ=255)

                                          cTATGTTCTAGGCTTTTTATTATACTTATCAACATTACAACGCTTTCAAAG‑TGTATTTTTTAG‑G‑‑‑C‑‑‑‑TTGT‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑GTAAAAAATAAATTAAATTTTTGATAAAACAAAACT‑TTAT‑Taa                                                                                       >  2:2205396/1‑128 (MQ=255)

                                           tagaacatAGGCTTTTTATTATACTTATCAACATTACAACGCTTTCAAAG‑TGTATTTTTTAG‑G‑‑‑C‑‑‑‑TTGT‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑GTAAAAATAAAAGTAATATTTGGATCAACCAAAACT‑AAAT‑aaaa                                                                                      >  2:136414/8‑128 (MQ=255)

                                            aTGTTCTAGGCTTTTTATTATACTTATCAACATTACAACGCTTTCAAAG‑TGTATTTTTTAG‑G‑‑‑C‑‑‑‑TTGT‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑TTAAAAATAAAATTAAGATTTGGATAAACCaaaaaaattaaaaaa                                                                                       >  2:4021113/1‑117 (MQ=255)

                                            aTGTTCTAGGCTTTTTATTATACTTATCAACATTACAACGCTTTCAAAG‑TGTATTTTTTAG‑G‑‑‑C‑‑‑‑TTG‑‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑TTAAAAATAAAAGTAAGATTTGGATAAACCAAAACT‑ATAT‑Caaaaa                                                                                    >  2:1590291/1‑128 (MQ=255)

                                              aacatAGGCTTTTTATTATACTTATCAACATTACAACGCTTTCAAAG‑TGTATTTTTTAG‑G‑‑‑C‑‑‑‑TTG‑‑‑‑‑‑‑TTTTTTTTTTTTT‑TT‑‑‑TAAAAATAAAATTAAGATTTGTTTCAACCAAAACT‑AAAT‑aaaaaaaaat                                                                                >  1:1008672/5‑126 (MQ=255)

                                                  tAGGCTTTTTATTATACTTATCAACATTACAACGCTTTCAAAG‑TGTATTTTTTAG‑G‑‑‑CTTGTTTG‑‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑GTAAAAATAAAAGTTAATTTTGGGTCACCCAAAACT‑ATAT‑Ctaaaaaa                                                                                  >  1:3882158/1‑128 (MQ=255)

                                                   aGGCTTTTTATTATACTTATCAACATTACAACGCTTTCAAAG‑TGTATTTTTTAG‑G‑‑‑C‑‑‑‑TTGT‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑TTAAAAATAAAAGTATGATTTGTTTGAACCAAAAAT‑AATT‑CAAAAAAATaaa                                                                              >  2:4335134/1‑128 (MQ=255)

                                                    ggCTTTTTATTATACTTATCAACATTACAACGCTTTCAAAG‑TGTATTTTTTAG‑G‑‑‑C‑‑‑‑TTGT‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑TTTAAAATAAATTTAATATTTTGATAGACCAAAACT‑ATAT‑AAAAAAAATTaaa                                                                             >  1:4611803/1‑128 (MQ=255)

                                                         tttATTATACTTATCAACATTACAACGCTTTCAAAG‑TGTATTTTTTAG‑G‑‑‑C‑‑‑‑TTGTTT‑‑‑‑TTTTTTTTTTTTT‑TTT‑GTAAAAATAAAATTAATATTTGTATAAACCAAAAAT‑ATTT‑AAAAAAAATAaaaaaa                                                                          >  1:1151102/1‑128 (MQ=255)

                                                         tttATTATACTTATCAACATTACAACGCTTTCAAAG‑TGTATTTTTTAG‑G‑‑‑C‑‑‑‑TTGTT‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑TTAAAATTAAAAATAATATTTTGATAAACAAAAAATAATAT‑‑AAAAAAATTaaaaaaa                                                                         >  2:1440606/1‑128 (MQ=255)

                                                           tattatACTTATCAACATTACAACGCTTTCAAAG‑TGTATTTTTTAG‑G‑‑‑C‑‑‑‑TTG‑‑‑‑‑‑‑TTTTTTTTTTTTT‑TTTATTAAAAATAAAAGTAATATTTTTAAAAAAAAATAAT‑TTAT‑AAAAAATATTAAAAAAATta                                                                      >  2:1145250/1‑128 (MQ=255)

                                                           tattatACTTATCAACATTACAACGCTTTCAAAG‑TGTATTTTTTAG‑G‑‑‑C‑‑‑‑TTG‑‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑GTTAAAATAAAATTAAAATTTTGTTATCACAAAAAT‑ATAT‑AAAAAATAATAAAAAAAAtaa                                                                     >  1:3020952/1‑127 (MQ=255)

                                                               ataCTTATCAACATTACAACGCTTTCAAAG‑TGTATTTTTTAG‑G‑‑‑C‑‑‑‑TTGTTT‑‑‑‑TTTTTTTTTTTTT‑TTT‑TTAAAAAAAAAATTTAATTTTGTTTAAACCAAAAAA‑TTTt‑aaaaaaaattaaaaaaaattat                                                                    >  2:3793718/1‑125 (MQ=255)

                                                                taCTTATCAACATTACAACGCTTTCAAAG‑TGTATTTTTTAG‑G‑‑‑C‑‑‑‑TTG‑‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑TTAAAAATAAAATTAATATTTGGATCAACCAAAAAT‑ATAT‑AAAAAAAATTAAAAAAAAtattttat                                                                >  2:4639417/1‑127 (MQ=255)

                                                                 aCTTATCAACATTACAACGCTTTCAAAG‑TGTATTTTTTAG‑G‑‑‑C‑‑‑‑TTGTTT‑‑‑‑TTTTTTTTTTTTT‑TTT‑TTAAAAATAAAATTAATATTTTGATAAACCAAAAAT‑ATAT‑AAAAAAAATTAAAAAAAAtatttt                                                                  >  2:255474/1‑128 (MQ=255)

                                                                 aCTTATCAACATTACAACGCTTTCAAAG‑TGTATTTTTTAG‑G‑‑‑C‑‑‑‑TTGTT‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑TTAAAAAAAAAATTAAGATTTGGAAAAACCAAAAAT‑ATAA‑AAAAAAAATTAAAAAAAAaatttta                                                                 >  2:4477949/1‑128 (MQ=255)

                                                                  cTTATCAACATTACAACGCTTTCAAAG‑TGTATTTTTTAG‑G‑‑‑C‑‑‑‑TTG‑‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑TTTAAAAAAAAAATAAAATATTTATAAAAAAAAAAT‑AAAT‑AAAAAAAATAAAAAAAAAtaatttatta                                                              >  2:2714545/1‑125 (MQ=255)

                                                                  cTTATCAACATTACAACGCTTTCAAAG‑TGTATTTTTTAG‑G‑‑‑C‑‑‑‑TTG‑‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑‑TAAAAAAAAAAATAAATATTTGTTAAAACAAAAAT‑AAAT‑AAAAAAAATTAAAAAAAAtaattttatat                                                             >  2:4392298/1‑123 (MQ=255)

                                                                  attattaAAATTTCAAATCTTTTAAAT‑TTTTTTTTTTAG‑T‑‑‑T‑‑‑‑TTT‑‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑GTAAAAATAAAAGTAAGATTTGGATCAACCAAAACT‑ATAT‑CAAAAAAATTAAAAAAAAtagtttagtt                                                              <  1:1590291/127‑1 (MQ=255)

                                                                   ttATCAACATTACAACGCTTTCAAAG‑TGTATTTTTTAG‑G‑‑‑C‑‑‑‑TTG‑‑‑‑‑‑‑TTTATTTTTTTTT‑TTT‑TTAAAAAAAAAAATTAAAATTTGATTAACAAAAAaa‑ataa‑aaaaaaaaataaaaaaaaaaatattcttt                                                             >  2:929750/1‑117 (MQ=255)

                                                                   ttATAAAAATTACAAAGCTTTAAAAT‑TGTATTTTTTTT‑G‑‑‑C‑‑‑‑TTGT‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑GTAAAAATAAAAGTAAGATTTGGATCAACCAAAACT‑ATAT‑CAAAAAAATTAAAAAAAAtagtttagtt                                                              <  1:4518193/128‑1 (MQ=255)

                                                                    tatataaaTTACAAATCTTTTAAAT‑TGTATTTTTTAG‑G‑‑‑T‑‑‑‑TTGT‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑GTAAAAATAAAAGTAAGATTTGGATCAACCAAAACT‑ATAT‑CAAAAAAATTAAAAAAAATAGTTTAGTTg                                                             <  2:3585104/121‑1 (MQ=255)

                                                                    tATCAACATTACAACGCTTTCAAAG‑TGTATTTTTTAG‑G‑‑‑C‑‑‑‑TTGTTT‑‑‑‑TTTTTTTTTTTTT‑TTT‑TTAAAAATAAAATTAAAATTTTGATAAACCAAAAAA‑TTAT‑AAAAAAAATTAAAAAAAAtatttttgt                                                               >  2:2506077/1‑128 (MQ=255)

                                                                    tATCAACATTACAACGCTTTCAAAG‑TGTATTTTTTAG‑G‑‑‑C‑‑‑‑TTG‑‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑GTAAAAATAAAATTAAGATTTGGATCAACCAAAAAT‑ATAA‑AAAAAAAATTAAAAAAAAAAGTTTATTTGa                                                            >  2:4228465/1‑128 (MQ=255)

                                                                     attaaaaTTTCAACGATTTAAAAT‑TGTATTTTTTAG‑GTGTG‑‑‑‑TTT‑‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑GTAAAAATAAAAGTAAGATTTGGATCAACCAAAACT‑ATAT‑CAAAAAAATTAAAAAAAAtagtttagtt                                                              <  1:4223489/122‑1 (MQ=255)

                                                                     aTCAACATTACAACGCTTTCAAAG‑TGTATTTTTTAG‑G‑‑‑C‑‑‑‑TTGTTTTTTTTTTTTTTTTTTTT‑TTT‑TTAAAAATAAAATTAATATTTGGTAAAACCAAAAAT‑TTAA‑AAAAAAAATTAAAAAAAAaatttt                                                                  >  2:577902/1‑128 (MQ=255)

                                                                       cAACATTACAACGCTTTCAAAG‑TGTATTTTTTAG‑G‑‑‑C‑‑‑‑TTG‑‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑TTAAAAATAAAATTAATATTTGTATCACCAAAAAGT‑TTTA‑AAAAAAAATTAAAAAAATTTTTTTATTTgtttt                                                         >  2:2899380/1‑124 (MQ=255)

                                                                           aTTACAACGCTTTCAAAG‑TGTATTTTTTAG‑G‑‑‑C‑‑‑‑TTG‑‑‑‑‑‑‑‑TTTTTTTTTTTT‑TTT‑GTAAAAAAAAAATTAATATTTTGTTAAACCAAAACT‑ATTA‑AAAAAAAATTAAAAAAAATAATTTATTTGAGTTTCaaa                                                    >  1:1801111/1‑128 (MQ=255)

                                                                            ttACAACGCTTTCAAAG‑TGTATTTTTTAG‑G‑‑‑C‑‑‑‑TTGTT‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑TTAAAAATAAAATTAATATTTGTTAAAACCAAAACT‑ATTA‑AAAAAAAAATAAAAAAAAAATTTTTGTTGAGTTTCa                                                      >  1:512462/1‑128 (MQ=255)

                                                                              aCAACGCTTTCAAAG‑TGTATTTTTTAG‑G‑‑‑C‑‑‑‑TTGT‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑TTAAAAAAAAAAGTAAGATTTGTTCAAACAAAAACA‑ATTA‑CAAAAAAATAAATAAAAATTATTTAATTTATTTTaaaaa                                                   >  1:2063271/1‑127 (MQ=255)

                                                                              aCAACGCTTTCAAAG‑TGTATTTTTTAG‑G‑‑‑C‑‑‑‑TTGT‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑GTAAAAATAAAAGTAATTTTTGGATCAACCAAAAAA‑AAAA‑AAAAAAAAATAAAAAAAATATTTTATTTGATTTTAAAAt                                                   >  2:2331416/1‑128 (MQ=255)

                                                                              aCAACGCTTTCAAAG‑TGTATTTTTTAG‑G‑‑‑C‑‑‑‑TAGT‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑TTAAAAAAAAAATTAATAATTGGATAAAAAAAAAAA‑AATT‑AAAAAAAATTAAAAAAAAATATTTAATTAAGTTTAAAAt                                                   >  1:1693263/1‑128 (MQ=255)

                                                                               cAACGCTTTCAAAG‑TGTATTTTTTAG‑G‑‑‑C‑‑‑‑TTGTTT‑‑‑‑TTTTTTTTTTTTT‑TTT‑TTAAAAATAAAAGTAATATTTGGATAAACCAAACAA‑TTAT‑CAAAAAAATTAAAAAAAAAAGTTTATTTTAGTTTCaaa                                                    >  2:1031770/1‑128 (MQ=255)

                                                                                aaCGCTTTCAAAG‑TGTATTTTTTAG‑G‑‑‑C‑‑‑‑TTGT‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑TTAAAAATAAAAGTAAGATTTGGATAAACCAAAAAT‑ATAT‑CAAAAAAATTAAAAAAAATATTTTATTTGAGTTTAAAAtat                                                 >  1:1427285/1‑126 (MQ=255)

                                                                                 aCGCTTTCAAAG‑TGTATTTTTTAG‑G‑‑‑C‑‑‑‑TTGTT‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑TTAAAAATAAAATTAATATTTTGTTAAACAAAAAAT‑AATT‑AAAAAAAATTAAAAAAAATATTTTATTTTTGTTTAAAAtgt                                                 >  1:3983774/1‑126 (MQ=255)

                                                                                   gCTTTCAAAG‑TGTATTTTTTAG‑G‑‑‑C‑‑‑‑TTGTTT‑‑‑‑TTTTTTTTTTTTT‑TTT‑TTAAAAATAAAATTAATATTTGGTTCAACCAAAACT‑ATAT‑CAAAAAAATTAAAAAAAATATTTTATTTTAGattacaaaatg                                                >  2:1214524/1‑120 (MQ=255)

                                                                                   gCTTTCAAAG‑TGTATTTTTTAG‑G‑‑‑C‑‑‑‑TTGTT‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑TTAAAAAAAAAATTTATATTTTGATCAAACAAAAAT‑AAAT‑AAAAAAAATTAAAAAAAATAGTTTTTGTTAGTTTCAAAtaaga                                               >  1:4093302/1‑124 (MQ=255)

                                                                                   gCTTTCAAAG‑TGTATTTTTTAG‑G‑‑‑C‑‑‑‑TTGT‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑TTAAAAATAAAAGTAATATTTTGATCAACCAAAAAT‑ATAT‑AAAAAAAATTAAAAAAAATATTTTATTTTATTTTCAAAtattca                                              >  2:2789148/1‑123 (MQ=255)

                                                                                    cTTTCAAAG‑TGTATTTTTTAG‑G‑‑‑C‑‑‑‑TTGT‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑TTTAAAAAAAAATTAAATTTTTGTTAAACAAAAAAA‑ATAT‑AAAAAAAAAAAAAAAAAATTTTTTATTTTTTTTTTAAATCTGCat                                             >  1:3916761/1‑128 (MQ=255)

                                                                                      ttCAAAG‑TGTATTTTTTAG‑G‑‑‑C‑‑‑‑TTGT‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑TTAAAAAAAAAATTAAAATTTTGATAAAACAAAAAT‑ATAT‑AAAAAAAATTAAAAAAAATTTTTTAATTTATTTTATAATAAGCatat                                           >  1:1213321/1‑128 (MQ=255)

                                                                                       tCAAAG‑TGTATTTTTTAG‑G‑‑‑C‑‑‑‑TTG‑‑‑‑‑‑‑‑TTTTTTTTTTTT‑TTT‑GTAAAAATAAAAGTAAGATTTGGATAAACCAAAACA‑ATAT‑CAAAAAAATTAAAAAAAATATTTTAGTTTTGTTTCAAAGATTCAtatata                                        >  2:3842712/1‑128 (MQ=255)

                                                                                          aaatttttttttttAT‑G‑‑‑T‑‑‑‑TTGT‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑GTAAAAATAAAAGTAAGATTTGGATCAACCAAAACT‑ATAT‑CAAAAAAATTAAAAAAAATAGTTTAGTTGAGTTTCAAATCTGCAtatata                                        <  2:3492645/120‑1 (MQ=255)

                                                                                              gtatttttttTG‑T‑‑‑T‑‑‑‑TTT‑‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑GTAAAAATAAAAGTAAGATTTGGATCAACCAAAACT‑ATAT‑CAAAAAAATTAAAAAAAATAGTTTAGTTGAGTTTCAAATCTGCATATATAGTACa                                   <  1:4279693/124‑1 (MQ=255)

                                                                                               gtATTTTTTAG‑G‑‑‑C‑‑‑‑TTG‑‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑GTAAAAATAAAAGTAAGTTTTGGATCAACCAAAACT‑ATAT‑CAAAAAAATTAAAAAAAATATTTTATTTGAGTTTAAAATATGCAAATATATTACaa                                  >  2:2655756/1‑128 (MQ=255)

                                                                                                 aTTTTTTAG‑G‑‑‑C‑‑‑‑TTGT‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑TTAAAAAAAAAATTAAGATTTGGAAAAACCAAAACA‑ATAT‑CAAAAAAATTAAAAAAAAATGTTTATTTTAGTTATAAATGTGAATATATAGTAaaaa                                 >  1:556554/1‑128 (MQ=255)

                                                                                                  ttttttAG‑G‑‑‑C‑‑‑‑TTGT‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑TTTAAAATAAAATTAATATTTGGATAAACCAAAAAT‑ATTT‑AAAAAAAATTAAAAAAAATATTTTATTTGAGTTTAAAATATGCATTTATATTAaaaat                                >  2:412308/1‑127 (MQ=255)

                                                                                                    ttttAG‑G‑‑‑C‑‑‑‑TTG‑‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑TTAAAAATAAAAGTAAGATTTGGATCAACCAAAAAT‑ATAT‑AAAAAAAAATAAAAAAAATAGTTTTATTTAGTTTAAAATCAGCATATAAATTAAAAATGaa                             >  2:1746582/1‑128 (MQ=255)

                                                                                                    ttttAG‑G‑‑‑C‑‑‑‑TTG‑‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑GTAAAAATAAAAGTAAGATTTGGATCAACCAAAACT‑ATAT‑CAAAAAAATTAAAAAAAATAGTTTAGTTGAGTTTAAAATCTGCATATATAGTACAAAGGaa                             >  2:1651208/1‑128 (MQ=255)

                                                                                                    ttttAG‑G‑‑‑C‑‑‑‑TTG‑‑‑‑‑‑‑‑TTTTTTTTTTTT‑TTT‑GTAAAAATAAAAGTAAGATTTGGATCAACCAAAAAT‑ATAT‑AAAAAAAATTAAAAAAAATAGTTTATTTTAGTTTCAAATCTTCATATATAGTAAAAATGAAt                            >  2:84214/1‑128 (MQ=255)

                                                                                                     tttAG‑G‑‑‑C‑‑‑‑TTG‑‑‑‑‑‑‑‑TTTTTTTTTTTT‑TTT‑GTAAAAATAAAAGTAAGATTTGGATCAACCAAAACT‑ATAT‑CAAAAAAATTAAAAAAAATATTTTATTTGAGTTTCAAATCTGCATATATATTACAAAGGAttt                           >  1:3575491/1‑128 (MQ=255)

                                                                                                      ttAG‑G‑‑‑C‑‑‑‑TTGT‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑GTAAAAATAAAAGTAAGATTTGGATCAACCAAAAAT‑ATTT‑AAAAAAAATTAAAAAAAATATTTTATTTGATTTTCAAAACTGCATTTATATTTAAAAGGAAt                            >  2:2965402/1‑128 (MQ=255)

                                                                                                               c‑‑‑‑TTG‑‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑GTAAAAAAAAAATTAAGATTTTGATCAACAAAAAAT‑ATTT‑AAAAAAAATTAAAAAAAATATTTTATTTGAGTTTCAAAAATAAAAATATAGTAAAAAGCAAttattca                      >  2:1979349/1‑123 (MQ=255)

                                                                                                                      tattttttttttttttttttt‑ttt‑GTAAAAATAAAAGTAAGATTTGGATCAACCAAAACT‑ATAT‑CAAAAAAATTAAAAAAAATAGTTTAGTTGAGTTTCAAATCTGCATATATAGTACAAAGGAAttt                          <  1:265460/120‑1 (MQ=255)

                                                                                                                      ggttttttttttttttttttt‑ttt‑GTAAAAATAAAAGTAAGATTTGGATCAACCAAAACT‑ATAT‑CAAAAAAATTAAAAAAAATAGTTTAGTTGAGTTTCAAATCTGCATATATAGTACAAAGGAAttt                          <  1:3611662/120‑1 (MQ=255)

                                                                                                                      ggcttgttttttttttttttt‑ttt‑ttAAAAAAAAAATAAATATTTGGGTAAACAAAAACT‑ATATAAAAAAAAATTAAAAATATTATTTTATTTTATTTTAAAATATTCATATAAAGTAAAAAGAAAtt                           >  2:1490680/9‑128 (MQ=255)

                                                                                                                        gtttgtttttttttttttt‑ttt‑GTAAAAATAAAAGTAAGATTTGTATCAACCAAAACT‑ATAT‑CAAAAAAATTAAAAAAAATAGTTTAGTTGAGTTTCAAATCTGCATATATAGTACAAAGGAATTTAt                        <  1:4170555/122‑1 (MQ=255)

                                                                                                                         tttgtgtttttttttttt‑ttt‑GTAATAATAATAGTAAGATTTGGATCAACAAAAACT‑ATAT‑CAAAAAAATTATAAAAAATAGTTTAGTTGAGTTTCAAATCTGCATATATAGTACAAAGGAATTTATg                       <  2:536996/122‑1 (MQ=255)

                                                                                                                          ttttttttttttttttt‑ttt‑GTAATATTTAAAGTAAGTTTTGGTTAAACTAAAAAT‑ATTT‑AAAAAAAATTAAATAAAATAGTTTAGTTTAGTTTCAAATCTTAATATATAGTACAAAGGAATTTATGa                      <  1:3302209/124‑1 (MQ=255)

                                                                                                                          ttgtttttttttttttt‑ttt‑ttAAAAATAAAATTAAGATTTGGATAAACCAAAAAT‑ATAT‑AAAAAAAATTAAAAAAAATATTTTATTTGATTTTCAAATATGCATATATAGTACAAATGAATTTATGa                      >  2:3985029/5‑128 (MQ=255)

                                                                                                                           tttttttttttttttt‑ttt‑GTAATAATAATAGTAAGATTTGGATCAACCAAAACT‑ATAT‑CATAAAAATTAAAAAAAATAGTTTAGTTGAGTTTCAAATCTGCATATATAGTACAAAGGAATTTATGAt                     <  2:3606168/125‑1 (MQ=255)

                                                                                                                               tttttttttttt‑ttt‑GTAAAAATAAAAGTTAGATTTTTATCAACCAAAACT‑ATAT‑CAAAAAAATTAAAAAAAATAGTTTAGTTGAGTTTCAAATCTGCATATATAGTACAAAGGAATTTATGATAGaa                 <  1:3816645/128‑1 (MQ=255)

                                                                                                                               tttttttttttt‑ttt‑GTAAAAATAAAAGTAAGATTTGGATCAACCAAAACT‑ATAT‑CAAAAAAATTAAAAAAAATAGTTTAGTTGAGTTTCAAATCTGCATATATAGTACAAAGGAATTTATGATAGaa                 <  1:3183986/128‑1 (MQ=255)

                                                                                                                                tttttttttat‑ttt‑GTAAAAATAAAATTAAGATTTGGATCAAACAAAATT‑ATAT‑CAAAAAAATTAAAAAAAATAGTTTAGTTGAGTTTCAAATCTGCATATATAGTACAAAGGAATTTATGATAGaaa                <  2:408392/128‑1 (MQ=255)

                                                                                                                                    ttttttt‑ttt‑ttAAATTTATTAGTTAGATTTTGATAAATCAAAAAT‑ATAT‑CAAAAATATTAAAAAAAATAGTTTAGTTGAGTTTCAAATCTGCATATATAGTACAAAGGAATTTATGATAGaaaaaaa            <  1:628223/128‑1 (MQ=255)

                                                                                                                                    ttttttt‑ttt‑GTAAAAATAAAAGTAAGATTTGGATCAAACAAAACT‑ATAT‑CAAAAAAATTAAAAAAAATAGTTTAGTTGAGTTTCAAATCTGCATATATAGTACAAAGGAATTTATGATAGaaaaaaa            <  1:852196/128‑1 (MQ=255)

                                                                                                                                      ttttt‑ttt‑ttAAAAATAAAATTAATTTTTTGATAAACCAAAAAT‑ATAT‑CATAAAATTTAAAAAAATTAGTTTAGTGGAGTTTCAAATCTGCATATATAGTACAAAGGAATTTATGATAGAAAAAAATg          <  2:2404190/128‑1 (MQ=255)

                                                                                                                                      ttttt‑ttt‑GTAAAAATAAAATTAAGATTTGGGTCAACCAAAAAT‑ATTT‑AAAAAAAATTAAAAAAAATATTTTTATTTATTTTAAAATATGAATATATATTAAAAAGGAATTTATTATATAAAAAAAtt          >  2:394185/1‑127 (MQ=255)

                                                                                                                                       tttt‑ttt‑GTAAAAATAAAAGTAAGATTTGGATCAACCAAAACT‑ATAT‑CAAAAAAATTAAAAAAAATAGTTTAGTTGAGTTTCAAATCTGCATATATAGTACAAAGGAATTTATGATAGAAAAAAATgg         <  1:2582524/128‑1 (MQ=255)

                                                                                                                                        ttt‑ttt‑ttAAAAATAAAATTAAAATTTGGATAAACCAAAAAT‑AAAT‑AAAAAAAATTAAAAAAAATAGTTTATTTTATTTTAAAATATGGATATATAGTAAAAAGGAATTTATGATATAAAAAAATTGa        >  1:4339378/1‑128 (MQ=255)

                                                                                                                                        ttt‑ttt‑GTAAAAATAAAATTAAGATTTGGATCAACCAAAACT‑ATAT‑AAAAAAAATTAAAAAAAATAGTTTAGTTGAGTTTCAAATCTGCATATATAGTACAAAGGAATTTATGATAGAAAAAAATGGa        <  1:2467621/128‑1 (MQ=255)

                                                                                                                                            ttt‑tTAATAATAAAATTAATATTTGGATCAAACAAAAAT‑ATAT‑CAAAAAAATTAAAAAAAATAGTTTAGTTGAGTTTCAAATCTGCATATATAGTACAAAGGAATTTATGATAGAAAAAAATGGATTa     <  1:936684/128‑1 (MQ=255)

                                                                                                                                            ttt‑GTAATAATAAAAGTTAGATTTTGATCAACAAAAACT‑ATAT‑CAAAAAAATTAAAAAAAATAGTTTAGTTGAGTTTCAAATCTGCATATATAGTACAAAGGTATTTATGATAGAAAAAAATGGATTa     <  1:1346739/128‑1 (MQ=255)

                                                                                                                                            ttt‑GTAAAAATAAAAGTTAGTTTTGTATCAACCATAACT‑ATAT‑AAAAAAAATTAAAAAAAATAGTTTAGTTGAGTTTCAAATCTGCATATATAGTACAAAGGAATTTATGATAGAAAAAAATGGATTa     <  2:3044547/128‑1 (MQ=255)

                                                                                                                                            ttt‑GTAAAAATAAAAGTAAGATTTGGATCAAACAAAAAT‑ATAT‑CAAAAAAATTAAAAAAAATAGTTTAGTTGAGTTTCAAATCTGCATATATAGTACAAAGGAATTTATGATAGAAAAAAATGGATTa     <  2:3880356/128‑1 (MQ=255)

                                                                                                                                             tt‑GTAAATTTATAAGTAATATTTTGTTTATTAAATACT‑ATAT‑ATAAAATATTAAAAAAAATAGTTTAGTTGAGTTTCAAATATGCATATATAGTACAAAGGAATTTATGATAGAAAAAAATGGATTat    <  1:1979349/128‑2 (MQ=255)

AAAAAATAAAAATTTTTCATTTAAAGCCATTTGTAAAAAAGCCTATGTTCTAGGCTTTTTATTATACTTATCAACATTACAACGCTTTCAAAG‑TGTATTTTTTAG‑G‑‑‑C‑‑‑‑TTG‑‑‑‑‑‑‑TTTTTTTTTTTTT‑TTT‑GTAAAAATAAAAGTAAGATTTGGATCAACCAAAACT‑ATAT‑CAAAAAAATTAAAAAAAATAGTTTAGTTGAGTTTCAAATCTGCATATATAGTACAAAGGAATTTATGATAGAAAAAAATGGATTAGTT  >  CP016816/203786‑204039

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 15 ≤` `ATCG/ATCG` `< 21 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`