BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,009,779	T→G	5.3%	L204V (TTG→GTG)	pldB →	lysophospholipase L2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,009,779	0	T	G	5.3%	304.7 / 5.1	105	L204V (TTG→GTG)	pldB	lysophospholipase L2
Reads supporting (aligned to +/- strand): ref base T (42/56); new base G (3/3); total (46/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 107 total aligned reads displayed.

GTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTCTGGCTGGCGAACAGGTGTTAG  >  NC_000913/4009630‑4009926

gtgtATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGt                                                                                                                                                     <  1:789103/150‑1 (MQ=255)

gtgtATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGt                                                                                                                                                     >  1:815799/1‑150 (MQ=255)

  gtATGTGACGCCATTGCTCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCCACTGGGCCGAAGCGCATCCACGTGTCCGTGATGGCTATGCAATAGGCACCTGGCGCTGGAGCGCGTTg                                                                                                                                                   >  1:274522/1‑150 (MQ=255)

  gtATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTg                                                                                                                                                   <  2:232753/150‑1 (MQ=255)

   tATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGc                                                                                                                                                  <  1:469471/150‑1 (MQ=255)

    aTGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGcc                                                                                                                                                 >  2:1078187/1‑150 (MQ=255)

         ccGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTg                                                                                                                                            <  2:590608/149‑1 (MQ=255)

          cGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGc                                                                                                                                           >  1:1246513/1‑150 (MQ=255)

          cGCCATTCCGCTAACTGCGCCACTGTTTGGGATCGTGATTAGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGCAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGc                                                                                                                                           <  2:851540/150‑1 (MQ=255)

           gCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGTTTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCAACGCGCGCCAGCGCCCGGTGCCTATTGca                                                                                                                                          <  1:1534525/150‑2 (MQ=255)

            ccATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCTATTGcat                                                                                                                                         <  2:825934/150‑3 (MQ=255)

                     tAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTAc                                                                                                                                <  2:517082/150‑1 (MQ=255)

                     tAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTAc                                                                                                                                <  2:873684/150‑1 (MQ=255)

                        cTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTgg                                                                                                                             >  1:1092490/1‑149 (MQ=255)

                        cTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGa                                                                                                                             >  1:1436634/1‑150 (MQ=255)

                            gcCAATGTTTGGGATCGAGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTAGTGAcccc                                                                                                                         >  1:350032/1‑149 (MQ=255)

                             ccAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCcac                                                                                                                        >  2:1003965/1‑150 (MQ=255)

                             ccAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATACTCAACTGGGCCGAAGCGCATACACGTTTCCGTGATGGCTATGCAATAGGCACCGTGCGCTTGCGCGCGTTGCCGGTTGCTATCAACGTACTGaacaac                                                                                                                        >  1:627615/1‑145 (MQ=255)

                              cAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCAGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCAGTTTGCTATCAACGTACTGACCcaca                                                                                                                       >  1:864211/1‑150 (MQ=255)

                               aaTGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCAcag                                                                                                                      >  2:432269/1‑150 (MQ=255)

                                  gTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCAcagcag                                                                                                                   >  1:1271078/1‑150 (MQ=255)

                                  gTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGAGTTGCCGTTTGCTATCAACGTCCTGCCCCAcagcag                                                                                                                   >  1:866656/1‑150 (MQ=255)

                                      gtgATCGTGATTCGTATGCCGTCATTTATGGCACGGAAGATCCTCATCGGGGCCGAAGCGCATCCACGTTTCCGTGATTGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAg                                                                                                               <  1:338407/148‑1 (MQ=255)

                                       ggATCGTGATTCGTATGCCGTCATTTATGGCCAGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGAACCGTGCGCTGGCTCGCGTTGCCGGTTGCTATCCACGTACTGACCCACAGCAGACAGc                                                                                                              >  1:458568/1‑150 (MQ=255)

                                        gATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCg                                                                                                             >  2:634361/1‑150 (MQ=255)

                                        gATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGACCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCg                                                                                                             <  2:890014/150‑1 (MQ=255)

                                            gTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGatat                                                                                                         <  1:463221/150‑1 (MQ=255)

                                               aTTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGa                                                                                                      <  1:1678440/150‑1 (MQ=255)

                                                   gTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTa                                                                                                  <  2:351967/150‑1 (MQ=255)

                                                   gTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTa                                                                                                  <  2:1052504/150‑1 (MQ=255)

                                                    tATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTaa                                                                                                 >  1:1661060/1‑150 (MQ=255)

                                                      tGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGCTTTCGACGTAACt                                                                                               >  2:1112038/1‑150 (MQ=255)

                                                         cGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTAc                                                                                            <  2:299126/150‑1 (MQ=255)

                                                          gTCATTTATGGCACGGCAGATCCTCATCTGGGCCTATGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACg                                                                                           <  1:67377/150‑1 (MQ=255)

                                                           tCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGc                                                                                          <  1:640858/150‑1 (MQ=255)

                                                            cATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGCCCCACAGCAGACAGCGATATCGACGTAACTTACGCt                                                                                         >  1:511096/1‑150 (MQ=255)

                                                            cATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCt                                                                                         <  1:975079/150‑1 (MQ=255)

                                                             aTTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCtt                                                                                        <  2:754461/150‑1 (MQ=255)

                                                              tttATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTc                                                                                       <  1:45242/150‑1 (MQ=255)

                                                               ttATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCt                                                                                      <  2:375603/150‑1 (MQ=255)

                                                                tATGGCATGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCCCCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGAGATCGACGTACCTTACGCTTCTa                                                                                     >  1:688257/1‑150 (MQ=255)

                                                                 aTGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTAt                                                                                    <  2:616542/150‑1 (MQ=255)

                                                                 aTGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTAt                                                                                    <  1:895626/150‑1 (MQ=255)

                                                                  tGGCACGGCAGATCCTCAACTGTGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACTTAACTTACGCTTCTATg                                                                                   <  2:404534/150‑1 (MQ=255)

                                                                  tGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGGTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATg                                                                                   >  1:765435/1‑150 (MQ=255)

                                                                   ggCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAAGGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGc                                                                                  <  2:1204465/150‑1 (MQ=255)

                                                                       cGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGTCTATGCAATAGGCACCTGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTTCTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCtgat                                                                              <  1:775524/150‑1 (MQ=255)

                                                                        ggCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCtgatg                                                                             <  2:1092490/150‑1 (MQ=255)

                                                                          cAGATCCTAAACTGGGCCGAAGCGCATCCACGTTTCCGTGAAGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGCTGAc                                                                           >  2:1593264/1‑150 (MQ=255)

                                                                           aGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGAcc                                                                          <  1:1095966/150‑1 (MQ=255)

                                                                           aGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGAcc                                                                          <  2:1472966/150‑1 (MQ=255)

                                                                           aGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGAcc                                                                          >  1:714968/1‑150 (MQ=255)

                                                                           aGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATAGACGTAACTTACGCTTCTATGCTGATGAcc                                                                          >  1:708753/1‑150 (MQ=255)

                                                                               ccTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAAc                                                                      <  1:11971/150‑1 (MQ=255)

                                                                                cTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGCCCCACAGCACACAGCGATTTCGACGTAACTTACGCTTCTATGCTGATGACCCAACg                                                                     >  2:900026/1‑150 (MQ=255)

                                                                                    aCTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTc                                                                 <  2:807481/150‑1 (MQ=255)

                                                                                     cTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGGTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTcg                                                                <  1:989527/150‑1 (MQ=255)

                                                                                       gggCCGAAGCGCATCCACGTTTCCGTGATGTCTATGCAATAGGCACCGGGCGCTGACGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTcgcg                                                              <  2:1603914/150‑1 (MQ=255)

                                                                                       gggCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTcgcg                                                              <  2:1083453/150‑1 (MQ=255)

                                                                                        ggCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGt                                                             <  2:1273652/150‑1 (MQ=255)

                                                                                         gCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCCCAGCAGACAGCGATATCGACGTAACTTTCGCTTCTATGCTGATGACCCAACGATTCGCGTc                                                            >  2:951268/1‑150 (MQ=255)

                                                                                         gCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTc                                                            <  2:472585/150‑1 (MQ=255)

                                                                                             aaGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTg                                                        <  2:1511193/150‑1 (MQ=255)

                                                                                              aGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCCGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTgg                                                       >  1:1547320/1‑150 (MQ=255)

                                                                                              aGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTgg                                                       >  2:1429359/1‑150 (MQ=255)

                                                                                              aGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTgg                                                       >  1:389771/1‑150 (MQ=255)

                                                                                               gcgcATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTggg                                                      <  2:815799/150‑1 (MQ=255)

                                                                                                 gcATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACTTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGcc                                                    >  2:1052710/1‑150 (MQ=255)

                                                                                                 gcATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGcc                                                    <  2:1204102/150‑1 (MQ=255)

                                                                                                  cATCCACGGTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGTCGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCg                                                   <  2:1135734/150‑1 (MQ=255)

                                                                                                   aTCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTAGCAACGTACTGACCCACAGCAGACAGCGATATCTACGTAACTTACGCTTCAATGCTGATGACCCAACGATTCGCGTCCGTGGGACGa                                                  >  2:882298/1‑150 (MQ=255)

                                                                                                    tCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGATGCCGTTTGCTATCAACGTACTGACCCACAGCCGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGAc                                                 >  1:1590386/1‑150 (MQ=255)

                                                                                                             ccGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTg                                        <  1:778849/150‑1 (MQ=255)

                                                                                                              cGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATAAACGTACTGACCCACAGCAGACAGCGATATCTACGTAAATTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTgg                                       <  2:627615/150‑1 (MQ=255)

                                                                                                               gTGATGGCTATGCAATCGGCACCGGGCGCGGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACCGCGATATCGACGAAACTTACGCTTCTATGCTGATGCCCCAACGCTTCGCGTCGGTGGGCCGACCTACCATTggg                                      <  2:334756/150‑1 (MQ=255)

                                                                                                                tGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCATCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGt                                     <  1:807075/150‑1 (MQ=255)

                                                                                                                   tGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTAcg                                  >  2:566726/1‑150 (MQ=255)

                                                                                                                     gCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTAcgcg                                >  2:1367957/1‑150 (MQ=255)

                                                                                                                     gCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCCATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGCGGGGCCGACCTACCATTGGGTAcgcg                                >  2:345747/1‑150 (MQ=255)

                                                                                                                         tGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAg                            <  2:1479340/150‑1 (MQ=255)

                                                                                                                         tGCAATAGGCACCGGGCGCTGGCGAGGGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGGCGGTGGGCCGACCTACCATTGGGTACGCGAAAg                            <  2:350032/150‑1 (MQ=255)

                                                                                                                         agcaaacggCAACGCGCGCCAGCGCCCGGTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAg                            >  1:901019/8‑150 (MQ=255)

                                                                                                                           caaacggCAACGCGCGCCAGCGCCCGGTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTa                          >  1:906802/6‑150 (MQ=255)

                                                                                                                           cAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTa                          >  1:1677929/1‑150 (MQ=255)

                                                                                                                              gcggCACCGCGCGCCAGCGCCCGGTGCCGGTTGCTATCACCGTACTGACCCACAGCAGACAGCGATATCGCCGTAACTTACGCTTCTATGCTGCTCCCCCCACGCTTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTc                       >  1:75321/3‑150 (MQ=255)

                                                                                                                               aGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTct                      <  2:864211/150‑1 (MQ=255)

                                                                                                                                ggCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTCTGCTGATGACCCATCGATTCGCGTCGGTGGGCCGACCTAACATTGGGTACGAGAAAGTATTctg                     >  1:227415/1‑150 (MQ=255)

                                                                                                                                  cACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTctggc                   <  1:634361/150‑1 (MQ=255)

                                                                                                                                    ccGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCTCAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTctggctg                 >  2:1380626/1‑150 (MQ=255)

                                                                                                                                    ccGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTctggctg                 <  2:1233642/150‑1 (MQ=255)

                                                                                                                                       ggCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGCCCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGCTGACCCAACGATTCGCGTCGGTGGGCCGACCTACCAGTGGGTACGCGAAAGTATTCTGGCTGGCg              >  1:1360832/1‑150 (MQ=255)

                                                                                                                                          gcTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTCTGGCTGGCGAAc           <  1:6765/150‑1 (MQ=255)

                                                                                                                                             ggCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCACCGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTCTGGCTGGCGAACAgg        >  1:233070/1‑150 (MQ=255)

                                                                                                                                             ggCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTCTGGCTGGCGAACAgg        >  1:96292/1‑150 (MQ=255)

                                                                                                                                               cgcgcgTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTCTGGCTGGCGAACAGgtg      <  2:783574/150‑1 (MQ=255)

                                                                                                                                               cgcgcgTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTCTGGCTGGCGAACAGgtg      >  2:57245/1‑150 (MQ=255)

                                                                                                                                                 cgcgTTGCCGTTTGCTATCAACGTACTGACCCTCAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTCTGGCTGGCGAACAGGTGtt    <  1:368611/150‑1 (MQ=255)

                                                                                                                                                 cgcgTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTCTGGCTGGCGAACAGGTGtt    <  2:1014513/150‑1 (MQ=255)

                                                                                                                                                   ggTTGCCGTTTGATATCAACGTACTGACCCAGAGCAGACAGGGATATCGACGTTATTTACGCTTCTTGGCTGATGCGCCATCGATTCGCGTCGGTGGGCGGACCTAGCATTGGGTACGCGAAAGTATTCTGGCTGGCGAACAGGTGTTAg  <  1:324929/149‑1 (MQ=255)

                                                                                                                                                   cgTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTCTGGCTGGCGAACAGGTGTTAg  >  2:1676672/1‑150 (MQ=255)

GTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCGTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTCTGGCTGGCGAACAGGTGTTAG  >  NC_000913/4009630‑4009926

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