BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,009,765	G→C	5.5%	G199A (GGG→GCG)	pldB →	lysophospholipase L2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,009,765	0	G	C	5.5%	927.4 / 30.0	269	G199A (GGG→GCG)	pldB	lysophospholipase L2
Reads supporting (aligned to +/- strand): ref base G (129/125); new base C (5/10); total (134/135)
Fisher's exact test for biased strand distribution p-value = 2.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

Only 100 of 274 total aligned reads displayed.

GCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTCTGGCTGGCG  >  NC_000913/4009616‑4009914

        aTCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGgc                                                                                                                                                <  1:921218/150‑1 (MQ=255)

         tCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGgcg                                                                                                                                               >  1:1174358/1‑150 (MQ=255)

           cAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGgcgcg                                                                                                                                             <  2:1162052/150‑1 (MQ=255)

           cAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGCGCGCCAgcgcc                                                                                                                                             >  1:1993216/1‑149 (MQ=255)

            aGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGgcgcgc                                                                                                                                            <  2:1214478/150‑1 (MQ=255)

             ggTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGcgcgcg                                                                                                                                           >  2:1852976/1‑150 (MQ=255)

                gtATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTg                                                                                                                                        >  1:904246/1‑150 (MQ=255)

                gtATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTg                                                                                                                                        >  2:350616/1‑150 (MQ=255)

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

                 tATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCAACGCGCGCCAGCGCCCGGTGc                                                                                                                                       <  2:548393/150‑1 (MQ=255)

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

                   tgtgACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGGTTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCg                                                                                                                                     <  2:582269/150‑1 (MQ=255)

                    gtgACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑t                                                                                                                                   >  2:241943/1‑150 (MQ=255)

                    gtgACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑t                                                                                                                                   >  2:492692/1‑150 (MQ=255)

                        cGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCTGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGc                                                                                                                               <  2:571387/150‑1 (MQ=255)

                        cGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGc                                                                                                                               <  2:1049451/150‑1 (MQ=255)

                        cGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGc                                                                                                                               >  2:1552948/1‑150 (MQ=255)

                          ccATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAAACGGCAACGCGCGCCAGCGCCCGGTGCCT‑ATTGcat                                                                                                                             <  2:911093/150‑3 (MQ=37)

                            aTTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTTCCG‑TTTGCTATc                                                                                                                           <  2:1485894/150‑1 (MQ=255)

                            aTTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATc                                                                                                                           <  1:628565/150‑1 (MQ=255)

                             ttGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAAACGGCAACGCGCTCCAGCGCCCGGTGCCT‑ATTGcatagc                                                                                                                          <  2:1459673/150‑6 (MQ=37)

                               gcgcTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAAc                                                                                                                        >  2:508685/1‑150 (MQ=255)

                                cgcTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACg                                                                                                                       >  1:1437119/1‑150 (MQ=255)

                                  cTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTa                                                                                                                     >  1:1291083/1‑150 (MQ=255)

                                    aaCTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACt                                                                                                                   >  1:1916686/1‑150 (MQ=255)

                                        gcgcCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGAcc                                                                                                               <  2:608720/150‑1 (MQ=255)

                                         cgcCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGAccc                                                                                                              <  2:746926/150‑1 (MQ=255)

                                          gcCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCca                                                                                                             >  2:1588252/1‑150 (MQ=255)

                                            cAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCACGCGTTGCCG‑TTTGCTATCAACGTACTGACCcaca                                                                                                           <  1:1706889/150‑1 (MQ=255)

                                              aTGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCAcagc                                                                                                         <  1:1869114/150‑1 (MQ=255)

                                              aTGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCAcagc                                                                                                         <  2:1380039/150‑1 (MQ=255)

                                              aTGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGA‑CCAcagca                                                                                                        <  2:922917/150‑1 (MQ=255)

                                               tGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCAcagca                                                                                                        <  1:69192/150‑1 (MQ=255)

                                                gTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCcac                                                                                                            <  2:170794/145‑1 (MQ=255)

                                                gTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCAcagcag                                                                                                       <  1:1742939/150‑1 (MQ=255)

                                                     ggATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCcaca                                                                                                           <  2:197413/141‑1 (MQ=255)

                                                      gATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCg                                                                                                 >  2:1387818/1‑150 (MQ=255)

                                                      gATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCg                                                                                                 <  1:625596/150‑1 (MQ=255)

                                                      gATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCg                                                                                                 >  1:270603/1‑150 (MQ=255)

                                                      gATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCg                                                                                                 >  1:1470000/1‑150 (MQ=255)

                                                         cGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGata                                                                                              >  2:1944841/1‑150 (MQ=255)

                                                         cGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGata                                                                                              >  2:1678060/1‑150 (MQ=255)

                                                         cGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGata                                                                                              <  1:934939/150‑1 (MQ=255)

                                                          gTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGatat                                                                                             >  1:1260912/1‑150 (MQ=255)

                                                           tGATTCGTATGCCGTCATTTATGGCACTCCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTACGTCATGGCTATGCAATAGGCACCGGGCACTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATc                                                                                            >  2:110320/1‑150 (MQ=255)

                                                            gATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCg                                                                                           >  2:1731633/1‑150 (MQ=255)

                                                            gATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCg                                                                                           >  1:964986/1‑150 (MQ=255)

                                                             cttCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGGTTCCGTGATGGCTATGTAATAGGCACCGGGCGCTGGCGCGCGGTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGa                                                                                          <  2:1655316/149‑1 (MQ=255)

                                                             aTTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGa                                                                                          <  2:1412818/150‑1 (MQ=255)

                                                               tCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACg                                                                                        <  1:1157638/150‑1 (MQ=255)

                                                                 gTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTa                                                                                      >  2:1524406/1‑150 (MQ=255)

                                                                   aTGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAAc                                                                                    <  1:592913/150‑1 (MQ=255)

                                                                      ccGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTa                                                                                 <  1:74896/150‑1 (MQ=255)

                                                                        gTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACg                                                                               >  1:320284/1‑150 (MQ=255)

                                                                          cATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCt                                                                             >  2:853817/1‑150 (MQ=255)

                                                                            tttATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCGACGTTTCCGTTATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTc                                                                           <  1:1661503/150‑1 (MQ=255)

                                                                              tATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACt                                                                                                                   >  1:1398604/1‑108 (MQ=255)

                                                                                tGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATg                                                                       >  1:1034059/1‑150 (MQ=255)

                                                                                tGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATg                                                                       <  1:1552948/150‑1 (MQ=255)

                                                                                  gCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCt                                                                     >  1:1490562/1‑150 (MQ=255)

                                                                                  gCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCt                                                                     >  1:206147/1‑150 (MQ=255)

                                                                                    aCGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCtg                                                                    >  1:719808/1‑149 (MQ=255)

                                                                                     cGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCtgat                                                                  <  2:1296882/150‑1 (MQ=255)

                                                                                     cGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCtgat                                                                  >  1:1787349/1‑150 (MQ=255)

                                                                                          gATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGAccc                                                             <  2:1260912/150‑1 (MQ=255)

                                                                                          gATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGAccc                                                             <  1:1408476/150‑1 (MQ=255)

                                                                                              cTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACg                                                         >  2:298725/1‑150 (MQ=255)

                                                                                                cAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGAt                                                       >  2:1112674/1‑150 (MQ=255)

                                                                                                cAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGAt                                                       >  2:801304/1‑150 (MQ=255)

                                                                                                 aaCTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGAtt                                                      >  1:1358897/1‑150 (MQ=255)

                                                                                                                cATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCg                                       >  2:938204/1‑150 (MQ=255)

                                                                                                                   ccACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGAcc                                    >  2:1002451/1‑150 (MQ=255)

                                                                                                                   ccACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGAcc                                    <  1:485549/150‑1 (MQ=255)

                                                                                                                     aCGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTa                                  >  1:179704/1‑150 (MQ=255)

                                                                                                                      cGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTAc                                 >  1:548582/1‑150 (MQ=255)

                                                                                                                             gTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTggg                          <  1:1588252/150‑1 (MQ=255)

                                                                                                                             gTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTggg                          <  1:51169/150‑1 (MQ=255)

                                                                                                                              tgttGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGt                         >  1:1104159/4‑150 (MQ=255)

                                                                                                                               gATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCg                                       >  1:475720/1‑135 (MQ=255)

                                                                                                                                   gCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTAcgcg                    >  1:436778/1‑150 (MQ=255)

                                                                                                                                     tATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGaa                  <  1:1712213/150‑1 (MQ=255)

                                                                                                                                      aTGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTAcgc                     <  1:1955296/146‑1 (MQ=255)

                                                                                                                                       tGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATTTCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAg                <  2:1285802/150‑1 (MQ=255)

                                                                                                                                       tGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAg                <  2:758606/150‑1 (MQ=255)

                                                                                                                                       tGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAg                <  2:1787110/150‑1 (MQ=255)

                                                                                                                                       tGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAg                <  1:494984/150‑1 (MQ=255)

                                                                                                                                         cAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTa                        >  2:725647/1‑140 (MQ=255)

                                                                                                                                           aacggCAACGCGCGCCAGCGCCCGGTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTAtt            >  1:1528354/4‑150 (MQ=255)

                                                                                                                                           aacggCAACGCGCGCCAGCGCCCGGTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTAtt            >  1:970119/4‑150 (MQ=255)

                                                                                                                                            tAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTc           >  2:332890/1‑150 (MQ=255)

                                                                                                                                             aGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTct          <  2:541216/150‑1 (MQ=255)

                                                                                                                                              ggCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTAc                       >  1:1065231/1‑136 (MQ=255)

                                                                                                                                              ggCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTctg         >  2:1904826/1‑150 (MQ=255)

                                                                                                                                              ggCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTctg         >  2:977460/1‑150 (MQ=255)

                                                                                                                                               gCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTctgg        >  1:83473/1‑150 (MQ=255)

                                                                                                                                               gCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTctgg        <  1:1678060/150‑1 (MQ=255)

                                                                                                                                               gCAACGCGCGCCAGCGCCCGGTGCCGTTTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTctg         >  2:179272/1‑150 (MQ=255)

                                                                                                                                                cACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTctggc       >  1:1445276/1‑150 (MQ=255)

                                                                                                                                                  ccGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTctggctg     <  2:699668/150‑1 (MQ=255)

                                                                                                                                                     ggCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTCTGGCTGGCg  <  1:249725/150‑1 (MQ=255)

GCAACGCCATCCAGGTGTATGTGACGCCATTGCGCTAACTGCGCCAATGTTTGGGATCGTGATTCGTATGCCGTCATTTATGGCACGGCAGATCCTCAACTGGGCCGAAGCGCATCCACGTTTCCGTGATGGCTATGCAATAGGCACCGGGCGCTGGCGCGCGTTGCCG‑TTTGCTATCAACGTACTGACCCACAGCAGACAGCGATATCGACGTAACTTACGCTTCTATGCTGATGACCCAACGATTCGCGTCGGTGGGCCGACCTACCATTGGGTACGCGAAAGTATTCTGGCTGGCG  >  NC_000913/4009616‑4009914

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