BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,052,201	C→A	5.5%	F145L (TTC→TTA)	typA →	GTP‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,052,201	0	C	A	5.5%	1308.8 / 32.6	328	F145L (TTC→TTA)	typA	GTP‑binding protein
Reads supporting (aligned to +/- strand): ref base C (155/155); new base A (8/10); total (163/165)
Fisher's exact test for biased strand distribution p-value = 8.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 329 total aligned reads displayed.

TGCTGGTGGTTGACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATCGCGGGTCTGGACCACGAAGATATGGCGGAAGACATGACCCCGCTGTACCAGGCGATTGTTGACCACGTTCC  >  CP009273/4052053‑4052350

tgctgGTGGTTGACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCg                                                                                                                                                      <  2:4926855/150‑1 (MQ=255)

   gggtggtTGACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATc                                                                                                                                                   >  2:1790079/2‑150 (MQ=255)

      tggtTGACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTAACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTAGATCTGt                                                                                                                                                >  1:2544861/1‑150 (MQ=255)

       ggtTGACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGtt                                                                                                                                               <  2:3229104/150‑1 (MQ=255)

          tGACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGt                                                                                                                                            <  1:1676817/150‑1 (MQ=255)

                aTTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCt                                                                                                                                      >  2:2370924/1‑150 (MQ=255)

                aTTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCt                                                                                                                                      <  2:1653214/150‑1 (MQ=255)

                aTTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCt                                                                                                                                      <  2:2547289/150‑1 (MQ=255)

                 tttGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTc                                                                                                                                     >  2:1282508/1‑150 (MQ=255)

                  ttGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCg                                                                                                                                    >  1:2768289/1‑150 (MQ=255)

                   tGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTCTTCGATCTGTTCGTTAACCTCGa                                                                                                                                   >  2:3918231/1‑150 (MQ=255)

                       ggCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGAcgcg                                                                                                                               >  2:65055/1‑150 (MQ=255)

                         cccGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGcgac                                                                                                                             <  1:2401157/150‑1 (MQ=255)

                         cccGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGGTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGcgac                                                                                                                             >  2:2122581/1‑150 (MQ=255)

                          ccGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGcgacc                                                                                                                            >  2:2207601/1‑150 (MQ=255)

                            gATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGcgaccga                                                                                                                          >  1:93745/1‑150 (MQ=255)

                            gATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTAGATCTGTTCGTTAACCTCGACGcgaccga                                                                                                                          >  1:448286/1‑150 (MQ=255)

                            gATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGAGTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGcgaccga                                                                                                                          >  1:1351408/1‑150 (MQ=255)

                             aTGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGcgaccgac                                                                                                                         <  1:2615220/150‑1 (MQ=255)

                                    aaaCGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGagcagc                                                                                                                  >  2:1732588/1‑150 (MQ=255)

                                         cgcTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGAc                                                                                                             <  1:3179053/150‑1 (MQ=255)

                                          gcTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACt                                                                                                            <  1:2985544/150‑1 (MQ=255)

                                           cTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACtt                                                                                                           <  2:1070867/150‑1 (MQ=255)

                                            ttCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTc                                                                                                          >  2:2941480/1‑150 (MQ=255)

                                              cGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTccc                                                                                                        <  1:3359027/150‑1 (MQ=255)

                                                tAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGa                                                                                                      >  2:4612276/1‑150 (MQ=255)

                                                tAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGa                                                                                                      >  1:405153/1‑150 (MQ=255)

                                                tAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGa                                                                                                      <  1:1818889/150‑1 (MQ=255)

                                                 aaacAAAAAATCGTTTGCTTACGTCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGTTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGAt                                                                                                     <  1:519479/147‑1 (MQ=255)

                                                 aaCCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGAt                                                                                                     >  2:3026568/1‑150 (MQ=255)

                                                  aCCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTAGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATc                                                                                                    >  2:596975/1‑150 (MQ=255)

                                                    cAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGt                                                                                                  >  2:3290675/1‑150 (MQ=255)

                                                     aaaaaaGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGtt                                                                                                 >  1:5099165/1‑150 (MQ=255)

                                                     aaaaaaGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGtt                                                                                                 >  1:1916075/1‑150 (MQ=255)

                                                     aaaaaaGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGtt                                                                                                 >  2:4852480/1‑150 (MQ=255)

                                                      aaaaaGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGttt                                                                                                <  1:1776644/150‑1 (MQ=255)

                                                      aaaaaGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGttt                                                                                                >  2:4664457/1‑150 (MQ=255)

                                                      aaaaaGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGttt                                                                                                >  2:799971/1‑150 (MQ=255)

                                                       aaaaGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTa                                                                                               >  2:2529208/1‑150 (MQ=255)

                                                         aaGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACg                                                                                             >  2:2169353/1‑150 (MQ=255)

                                                          aGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGc                                                                                            >  1:10167/1‑150 (MQ=255)

                                                              tttGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCt                                                                                        >  1:3267339/1‑150 (MQ=255)

                                                               ttGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTg                                                                                       <  1:147127/150‑1 (MQ=255)

                                                                tGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTgc                                                                                      <  2:4449013/150‑1 (MQ=255)

                                                                 gCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTgcg                                                                                     >  2:748507/1‑150 (MQ=255)

                                                                  cTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTgcgc                                                                                    <  2:2620969/150‑1 (MQ=255)

                                                                  cTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTgcgc                                                                                    >  1:3825206/1‑150 (MQ=255)

                                                                   ttACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCt                                                                                   >  2:2467890/1‑150 (MQ=255)

                                                                     aCGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGa                                                                                 >  2:3777047/1‑150 (MQ=255)

                                                                      cGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGaa                                                                                <  2:2306300/150‑1 (MQ=255)

                                                                      cGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGaa                                                                                >  2:3478901/1‑150 (MQ=255)

                                                                       ggCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAAc                                                                               <  2:2897129/150‑1 (MQ=255)

                                                                         ccTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACgg                                                                             >  1:489875/1‑150 (MQ=255)

                                                                          cTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGt                                                                            >  2:3355749/1‑150 (MQ=255)

                                                                            gAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTAGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTAt                                                                          >  2:2782557/1‑150 (MQ=255)

                                                                               gcagATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATAGTTTACGCTTCTGCGCTGAACGGTATcgc                                                                       <  1:4147684/147‑1 (MQ=255)

                                                                                ccGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTCTcgcg                                                                      >  2:3597664/1‑150 (MQ=255)

                                                                                ccGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATcgcg                                                                      >  1:4317577/1‑150 (MQ=255)

                                                                                ccGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATcgcg                                                                      <  2:1474203/150‑1 (MQ=255)

                                                                                 cGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATCGCgg                                                                     >  2:1342381/1‑150 (MQ=255)

                                                                                 cGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATCGCgg                                                                     <  2:308889/150‑1 (MQ=255)

                                                                                  gATTGTTGTTATCAACAAAGTTGACCGCCCGGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATCTCggg                                                                    >  1:4406200/1‑150 (MQ=255)

                                                                                   aTTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATCGCGGGt                                                                   <  2:1863320/150‑1 (MQ=255)

                                                                                    ttgttgTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATCGCGGGTc                                                                  <  2:4202919/150‑1 (MQ=255)

                                                                                     tgttgtTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATCGCGGGTCt                                                                 >  2:3103682/1‑150 (MQ=255)

                                                                                     tgttgtTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATCGCGGGTCt                                                                 <  2:4932615/150‑1 (MQ=255)

                                                                                       ttgttATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATCGCGGGTCTgg                                                               <  2:841908/150‑1 (MQ=255)

                                                                                          ttATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATCGCGGGTCTGGAcc                                                            >  1:108773/1‑150 (MQ=255)

                                                                                          ttATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATCGCGGGTCTGGAcc                                                            >  2:573146/1‑150 (MQ=255)

                                                                                          ttATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATCGCGGGTCTGGAcc                                                            >  1:4828665/1‑150 (MQ=255)

                                                                                               aacaaAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATCGCGGGTCTGGACCACGaa                                                       >  2:4554930/1‑150 (MQ=255)

                                                                                                  aaaGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCTGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGTCGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATCGCGGGTCTGGACCACGAAGat                                                    <  1:3377133/150‑1 (MQ=255)

                                                                                                  aaaGTCGACCGCGCTGGAGCGCGTGCTGTTTGGGTTGTGGATCAGGTATTCGATCTGTTATTTTACCTCGTCGCGACCGACGAGCAGCTGGTCGTCCCTATCGTTTTCGCTTCTGCGCTTAACTGTATCTCGGGTATGGACCTCgatgat                                                    <  2:3505935/150‑1 (MQ=255)

                                                                                                    aGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATCGCGGGTCTGGACCACGAAGatat                                                  <  1:2847650/150‑1 (MQ=255)

                                                                                                    aGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATCGCGGGTCTGGACCACGAAGatat                                                  <  1:650717/150‑1 (MQ=255)

                                                                                                          ccGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATCGCGGGTCTGGACCACGAAGATATGGCGGa                                            <  2:2400310/150‑1 (MQ=255)

                                                                                                             cccTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATCGCGGGTCTGGACCACGAAGATATGGCGGAAGa                                         >  2:100403/1‑150 (MQ=255)

                                                                                                             cccTGGCGCGCGTCCTGAGTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATCGCGGGTCTGGACCACGAAGATATGGCGGAAGa                                         <  1:2774293/150‑1 (MQ=255)

                                                                                                                tGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATCGCGGGTCTGGACCACGAAGATATGGCGGAAGACAt                                      <  2:1499852/150‑1 (MQ=255)

                                                                                                                tGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATCGCGGGTCTGGACCACGAAGATATGGCGGAAGACAt                                      >  2:5029555/1‑150 (MQ=255)

                                                                                                                tGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATCGCGGGTCTGGACCACGAAGATATGGCGGAAGACAt                                      >  1:1928966/1‑150 (MQ=255)

                                                                                                                    gcgcgTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATCGCGGGTCTGGACCACGAAGATATGGCGGAAGACATGAcc                                  <  1:3096060/150‑1 (MQ=255)

                                                                                                                      gcgTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATCGCGGGTCTGGACCACGAAGATATGGCGGAAGACATGAcccc                                <  1:4302455/150‑1 (MQ=255)

                                                                                                                        gTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATCGCGGGTCTGGACCACGAAGATATGGCGGAAGACATGACCCCGc                              >  1:55845/1‑150 (MQ=255)

                                                                                                                         tCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATCGCGGGTCTGGACCACGAAGATATGGCGGAAGACATGACCCCGCt                             <  1:2207601/150‑1 (MQ=255)

                                                                                                                         tCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATCGCGGGTCTGGACCACGAAGATATGGCGGAAGACATGACCCCGCt                             <  2:2095508/150‑1 (MQ=255)

                                                                                                                         tCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATCGCGGGTCTGGACCACGAAGATATGGCGGAAGACATGACCCCGCt                             >  1:4997848/1‑150 (MQ=255)

                                                                                                                           cTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATCGCGGGTCTGGACCACGAAGATATGGCGGAAGACATGACCCCGCTGt                           >  2:4260428/1‑150 (MQ=255)

                                                                                                                                tGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATCGCGGGTCTGGACCACGAAGATATGGCGGAAGACATGACCCCGCTGTACCAg                      >  1:1722101/1‑150 (MQ=255)

                                                                                                                                   gTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATCGCGGGTCTGGACCACGAAGATATGGCGGAAGACATGACCCCGCTGTACCAGGCg                   <  1:511430/150‑1 (MQ=255)

                                                                                                                                    ttGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATCGCGGGTCTGGACCACGAAGATATGGCGGAAGACATGACCCCGCTGTACCAGGCGa                  >  1:1990666/1‑150 (MQ=255)

                                                                                                                                     tgtgGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATCGCGGGTCTGGACCACGAAGATATGGCGGAAGACATGACCCCGCTGTACCAGGCGAt                 >  1:3817141/1‑150 (MQ=255)

                                                                                                                                       tgGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATCGCGGGTCTGGACCACGAAGATATGGCGGAAGACATGACCCCGCTGTACCAGGCGAttg               >  1:1125536/1‑150 (MQ=255)

                                                                                                                                         gATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATCGCGGGTCTGGACCACGAAGATATGGCGGAAGACATGACCCCGCTGTACCAGGCGAttgtt             <  1:4627709/150‑1 (MQ=255)

                                                                                                                                         gATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATCGCGGGTCTGGACCACGAAGATATGGCGGAAGACATGACCCCGCTGTACCAGGCGAttgtt             <  2:3821064/150‑1 (MQ=255)

                                                                                                                                            cAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATCGCGGGTCTGGACCACGAAGATATGGCGGAAGACATGACCCCGCTGTACCAGGCGATTGTTGAc          <  2:4776580/150‑1 (MQ=255)

                                                                                                                                            cAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATCGCGGGTCTGGACCACGAAGATATGGCGGAAGACATGACCCCGCTGTACCAGGCGATTGTTGAc          >  1:1254029/1‑150 (MQ=255)

                                                                                                                                             aGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATCGCGGGTCTGGACCACGAAGATATGGCGGAAGACATGACCCCGCTGTACCAGGCGATTGTTGAcc         >  1:3654221/1‑150 (MQ=255)

                                                                                                                                              ggTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATCGCGGGTCTGGACCACGAAGATATGGCGGAAGACATGACCCCGCTGTACCAGGCGATTGTTGACCa        >  1:1402748/1‑150 (MQ=255)

                                                                                                                                               gTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATCGCGGGTCTGGACCACGAAGATATGGCGGAAGACATGACCCCGCTGTACCAGGCGATTGTTGACCAc       <  2:2730745/150‑1 (MQ=255)

TGCTGGTGGTTGACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATCGCGGGTCTGGACCACGAAGATATGGCGGAAGACATGACCCCGCTGTACCAGGCGATTGTTGACCACGTTCC  >  CP009273/4052053‑4052350

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